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Abstract 


PRAIRIE,  MICHAEL  WILLIAM.  Optical  Properties  of  Thermally  Intermixed 
Quantum  Wells  and  Superlattices.  (Under  the  direction  of  Robert  M.  Kolbas.) 

The  electronic  structure  and  optical  properties  of  thermally  intermixed 
GaAs/AlGaAs  quantum  wells  and  superlattices  are  examined  using  thermal  annealing  in 
a  controlled  environment,  time-integrated  and  time-resolved  photoluminescence  (PL 
and  TPL,  respectively),  and  theoretical  models  based  on  the  envelope-function 
approximation.  A  computer  model  to  calculate  the  energy  structure  of  intermixed 
quantum  wells  is  refined  and  used  with  PL  and  TPL  techniques  to  characterize  the 
luminescence  properties  of  thin  (1  monolayer  (ML)  to  8  ML)  annealed  quantum  wells. 
A  model  for  intermixed  superlattices  is  developed  next  which  predicts  the  changes  in 
energy  dispersion  and  density  of  states  with  annealing  time.  Both  PL  and  TPL  are  used 
to  verify  the  validity  of  the  model. 

Selective  intermixing  techniques  using  a  gallium  focused  ion  beam  (FIB),  SiN, 
encapsulation,  and  patterned  GaAs  grown  at  low  substrate  temperature  (low- 
temperature  GaAs,  or  LT-GaAs)  are  then  studied.  The  FIB  technique  tends  to  suppress 
layer  interdiffusion  in  and  around  the  implant  region.  The  SiN,  encapsulation  is  most 
promising  with  a  selectivity  ratio  of  more  than  4:1.  Initial  LT-GaAs  experiments  show 
a  small  selectivity  between  regions. 

Finally,  a  model  of  the  dielectric  constant  for  partially  intermixed  quantum  wells 
is  developed  using  a  modified  optical  matrix  element  which  takes  into  account  the 
electron-hole  wavefunction  overlap,  as  well  as  the  penetration  of  the  wavefunctions  into 
the  barriers  of  thin  and  intermixed  wells.  The  radiative  lifetime  is  calculated  from  the 
matrix  element,  and  compared  with  TPL  data.  The  results  suggest  the  transport  of 


crystal  defects  from  the  surface  during  the  annealing  process.  Next,  the  index  of 
refraction  of  partially  intermixed  wells  is  calculated,  and  a  3-dB  power  splitter  is 
designed  as  an  example  of  using  selective  intermixing  in  integrated  optics. 
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1  Introduction 


This  chapter  introduces  previous  work  in  thermally  intermixed  quantum  wells  and 
superlattices  and  provides  an  overview  of  the  work  done  in  this  thesis.  The  evolution 
of  semiconductor  quantum  structures  for  integrated  optics  is  presented  first,  and  serves 
as  the  motivation  for  studying  the  optical  properties  of  partially  intermixed  quantum 
wells.  Then,  an  overview  of  the  goals  and  an  outline  of  this  thesis  are  given. 

1.1  Motivation 

Ever  since  1974  when  Esaki  and  Chang  demonstrated  resonant  tunneling  across 
potential  barriers,i‘i  and  Dingle  et  al.  measured  the  quantized  energy  levels  in  a 
man-made  quantum  well,Pi  the  study  of  ultrathin  semiconductor  layers  has  exploded. 

As  crystal  growth  improved  and  new  materials  were  used  to  make  these  layers,  many 
new  electronic  structures  were  possible,  some  of  which  have  found  their  way  into 
electronic  and  optical  devices,  pi  The  underlying  advantage  of  these  structures  is  that 
the  optical  and  electronic  properties  can  be  tailored  by  choosing  appropriate  materials 
and  layer  thicknesses.  One  example  of  this  ability  to  engineer  the  properties  of  a 
semiconductor  device  is  the  quantum  well  heterostructure  laser  in  which  the  emission 
wavelength  can  be  selected  by  varying  the  thickness  of  the  well.w 

The  dielectric  response,  which  determines  the  index  of  refraction  and  the 
absorption  properties  of  semiconductors  is  also  dependent  on  the  electronic  structure. 
These  characteristics  can  be  adjusted  through  choices  of  materials  and  layer  thicknesses 
in  quantum  wells  and  superlattices.  Early  work  by  Streifer  et  al.  estimated  the  index  of 
refraction  of  a  multiple  quantum  well  laser  diode  (MQW-LD)  to  be  that  of  the  average 
AlGaAs  alloy  over  the  entire  superlattice  or  MQW-LD  structure,  pi  Ladig  et  a/,  pi  and 
Holonyak  et  a/., pi  however,  showed  that  the  index  of  a  superlattice  which  had  been 
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thermally  intermixed  into  an  average  alloy  was  lower  than  adjacent  areas  of  the 
superlattice  which  had  been  preserved.  This  effect  had  obvious  applications  to 
integrated  optics.  The  selective  intermixing  was  achieved  through  impurity-induced 
layer  disordering  (IILD),  which  was  used  by  Deppe  et  and  Thornton  et  a/.PJ  to 
fabricate  index-guided  MQW-LD's  which  exhibited  some  of  the  lowest  threshold 
currents  of  any  semiconductor  laser.  The  IILD  technique  was  later  used  by  Julien  et  al. 
to  fabricate  impurity-induced  disorder-delineated  passive  optical  waveguides. 

While  IILD  demonstrates  the  ability  of  selective  disordering  to  delineate  optical 
waveguides,  it  also  introduces  impurities  into  the  waveguide  structure.  The  optical 
field  is  then  attenuated  through  scattering  losses  from  three  factors:  (1)  the  increased 
number  of  defects;  (2)  increased  interband  absorption  due  to  the  shift  in  the  absorption 
edge  by  the  impurities,  and  (3)  increased  free  carrier  absorption  since  the  impurities 
generate  more  free  carriers.  In  order  to  make  low-loss  waveguides,  a  method  of 
selective  disordering  must  be  used  which  avoids  the  use  of  impurities.  One  way  to  do 
this  is  to  control  the  point  defects  generated  during  high-temperature  annealing  stq)s 
using  overpressures  of  the  same  elements  that  constitute  the  semiconductor  structure.  1*^1 
The  delineation  of  structures  using  this  impurity-free  method  requires  selective 
exposure  of  areas  of  the  semiconductor  to  these  overpressures.  This  can  be 
accomplished  by  using  a  dielectric  to  mask  selected  regions  for  thermal  annealing  in  an 
overpressure. Another  effective  method  is  to  implant  constituent  atomic  species  in 
or  near  the  heterostructure  to  mechanically  intermix  the  wells  or  to  create  a  local 
internal  "overpressure,"  respectively. 

The  study  of  selective  intermixing  for  use  in  integrated  optical  devices  is  in  its 
infancy,  and  a  great  deal  of  work  is  needed  before  the  technique  can  be  used  to 
fabricate  sophisticated  devices.  For  example,  the  index  of  refraction  and  electro-optic 
coefficients  of  intermixed  quantum  structures  must  be  known  precisely  in  order  to 
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determine  the  dimensions  and  process  conditions  to  build  a  directional-coupler  switch. 
While  data  and  models  exist  (and  continue  to  be  refined)  for  the  index  of  quantum  wells 
and  superlatticesf’-*®!  as  well  as  alloys,  we  would  be  naive  to  think  that  a 
selectively  intermixed  structure  would  be  made  up  of  regions  of  a  perfectly  preserved 
superlattice  which  abruptly  joins  to  regions  of  completely  intermixed  alloys.  In  fact  we 
would  expect  both  the  "preserved"  and  "intermixed"  regions  to  be  partially  intermixed 
as  shown  in  Figure  1.1,  the  differences  being  the  degree  to  which  they  are  intermixed. 
To  date,  no  attempt  to  quantify  the  dielectric  properties  of  partially  intermixed  quantum 
structures  has  been  published.  With  that  in  mind,  the  purpose  of  this  thesis  is  to  extend 
the  theories  mentioned  above  to  include  partially  intermixed  quantum  structures. 

1.2  Overview  of  Thesis 

In  this  thesis,  the  optical  properties  of  thermally  intermixed  GaAs/AlGaAs 
quantum  wells  and  superlattices  are  investigated.  Specifically,  the  effects  of  atomic 
interdiffusion  at  the  heterostructure  interfaces  on  the  optical  matrix  element  are 
examined,  and  a  theory  is  developed  which  describes  the  matrix  element,  and  hence, 
the  dielectric  properties  for  partially  intermixed  wells.  Several  interdiffusion 
experiments  were  performed  to  characterize  the  intermixing  process,  and  several 
existing  theories  were  applied  to  the  results  to  develop  a  new  model.  The  emphasis 
throughout  the  thesis  is  on  the  application  of  intermixed  quantum  wells  in  optical 
waveguides  for  integrated  optoelectronics. 

A  review  of  recent  research  which  is  relevant  to  this  work  is  presented  in 
Chapter  2.  A  record  of  the  work  performed  in  this  thesis  begins  in  Chapter  3,  where 
the  experimental  methods  are  described.  Chapters  4  and  5  contain  the  studies  of  the 
intermixing  process.  The  former  describes  the  work  performed  to  characterize  the 
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(b) 


Figure  1.1  Schematic  of  a  selectively  intermixed  quantum  well.  The  diffusion 
length  L  =  -[Wt  is  (a)  infinity  in  the  intermixed  region  and  zero  in  the  preserved 
region  for  the  ideal  case,  and  (b)  finite  for  the  realistic  case.  Here,  a  10;  1  ratio  of 
the  interdiffusion  coefficient  is  assumed  for  the  intermixed  and  preserved 
regions,  respectively,  (5  x  10->*:5  x  lO-*’)  for  a  one-hour  annealing  process. 
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interdiffusion  coefficient  between  GaAs/AlGaAs  heterostructure  interfaces,  and  the 
latter  examines  several  techniques  to  selectively  intermix  quantum  wells  and 
superlattices. 

The  optical  properties  of  intermixed  wells  are  addressed  in  Chapter  6,  where  an 
original  model  for  the  optical  matrix  element  in  intermixed  quantum  wells  is  presented. 
This  new  matrix  element  is  used  to  calculate  the  dielectric  constant,  which  then  allows 
for  the  first-ever  calculation  of  the  refractive  index  for  partially  intermixed  quantum 
wells. 

The  work  is  summarized  in  Chapter  7,  and  conclusions  are  drawn  regarding  the 
applicability  of  using  selectively  intermixed  quantum  wells  and  superlattices  in 
integrated  optoelectronic  devices.  Finally,  the  reference  list  appears  in  Chapter  8, 
followed  by  a  collection  of  appendices  describing  in  depth  some  of  the  basic  theories 
used  in  this  thesis. 
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2  Review  of  Thermal  Intermixing  and  Optical  Properties  of 
Quantum  Wells  and  Superlattices 


In  this  chapter,  recent  work  in  characterizing  the  mechanisms  of  thermal 
intermixing  of  quantum  wells  and  superlattices  is  presented  first.  The  development  of 
theories  to  describe  the  dielectric  constant,  or  optical  properties  of  quantum  structures 
(and  hence,  their  refractive  index  and  absorption  properties)  is  presented  separately. 
Although  brief  reference  has  been  made  to  selective  intermixing  to  produce  local 
changes  in  the  refractive  index,  these  two  methods  have  not  been  quantitatively 
combined  until  this  study  and  will  be  reported  in  subsequent  chapters  where  they  are 
integrated  for  the  first  time. 

2.1  Thermally  Intermixed  Quantum  Wells  and  Superlattices 

When  first  grown,  quantum  wells  and  superlattices  are  composed  of  layers  of 
different  materials  with  sharp  interfaces.  If  these  heterostructures  are  exposed  to  high 
enough  temperatures,  defects  are  generated  in  the  crystal  lattice  which  allow  the 
materials  at  the  heterojunctions  to  intermix.  This  results  in  a  graded  interface  between 
the  materials,  which  changes  the  shape  of  the  original  potential.  The  first  task  in 
controlling  the  intermixing  of  quantum  structures,  then,  is  to  characterize  the 
mechanisms  of  diffusion  and  find  a  way  to  control  them. 

2.1.1  Heterojunction  Intermixing 

The  AlGaAs  material  system  has  the  zinc  blende  crystral  structure  which  is 
composed  of  two  inter-penetrating  face-centered  cubic  sublattices.  The  Group-Ill 
sublattice  contains  the  gallium  and  aluminum  atoms,  while  the  Group-V  sublattice  is 
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made  up  entirely  of  arsenic  atoms.  Since  atomic  diffusion  on  the  Group- V  sublattice  is 
inconsequential,  the  intermixing  of  heterojunctions  depends  only  on  the  diffusion  of 
gallium  and  aluminum  atoms  on  the  Group-Ill  sublattice. 

Diffusion  of  the  atoms  on  the  Group-III  sublattice  occurs  via  defects  in  the  lattice. 
Some  important  defects  are  shown  in  Figure  2. 1  (a).  The  vacancy  (Vn,)  is  an  empty 
lattice  site,  the  interstitial  defect  (Im)  is  an  atom  belonging  to  the  Group-III  sublattice 
that  is  not  at  a  lattice  site,  and  the  substitutional  defect  can  be  either  a  Group-III  atom  in 
the  Group- V  sublattice  (niy)  or  vice-versa  (ASm).  These  basic  defects  are  combined  to 
make  other,  more  complex  defects.  For  example,  a  Frenkel  Pair  is  an  interstitial  atom 
which  is  associated  with  a  nearby  vacancy  where  that  atom  should  be.  The  diffusion 
process  via  Group-III  vacancies  occurs  when  a  Vm  is  filled  with  an  aluminum  or 
gallium  atom  from  the  nearest  Group-III  lattice  point,  or  it  may  be  filled  from  a  passing 
interstitial  atom.  Group-m  interstitials,  Im,  on  the  other  hand,  move  about  the  lattice 
until  switching  places  with  an  atom  in  a  Group-III  lattice  site,  thereby  creating  another 
In,,  which  may  be  an  Iq,  or  an  I^,.  The  substitutional,  or  antisite,  defect  process 
involves,  for  instance,  As^  which  migrates  through  the  crystal,  being  replaced  by  either 
Gam  o*"  ^  it  leaves  each  lattice  site.  Some  interdiffusion  mechanisms  are  shown  in 

Figure  2.1  (b).  The  rate  at  which  the  atoms  diffuse  depends  on  the  equilibrium 
concentration  of  defects.  Assuming  that  interstitial  atoms  and  vacancies  are  the 
dominating  defects,  their  concentrations  are  denoted  by  [  /  m  ]  and  [Km],  respectively. 
The  diffusion  coefficient,  D,,,,  can  then  be  described  by  the  simplified  modelP^i 

(2-1) 

where  /  i  and  /  2  contain  information  about  the  crystal  structure,  is  the  Group-III 
vacancy  diffusion  rate,  zxidD  is  the  Group-III  interstitial  diffusion  rate.  By 
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Figure  2.1  Schematic  of  the  interdiffusion  process,  (a)  The  major  point  defects 
are  shown,  including  the  Group-in  vacancy  (Vqj)  and  interstitial  (Iqj),  and  the 
antisite  defects  Asq,  (shaded  circle),  (b)  Schematic  of  (i)  vacancy,  (ii)  interstitial, 
and  (iii)  antisite  interdiffusion  processes.  The  filled  circles  can  be  considered  to 
be  gallium  atoms,  the  open  circles  aluminum,  and  shaded  circles  arsenic.  Only 
the  Group-m  sublattice  is  depicted. 
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controlling  [  K  ,„  ]  or  [  /  ],  the  diffusion  coefficient  can  be  controlled.  This  can  be  done 
by  adjusting  the  crystal  Fermi  levelP3.Mi  ^ith  impurities,  or  by  controlling  the 
stoichiometry  of  the  crystal  near  the  surface  with  various  overpressures. 

Because  V„j's  are  negatively  charged,  the  concentration  of  vacancies  tends  to 
follow  the  concentration  of  donor  impurities,  or  more  precisely,  the  crystal  Fermi  level. 
Since  ionized  donors  are  positively  charged,  more  vacancies  will  be  created  by  the 
crystal  in  order  to  preserve  charge  neutrality.  The  argument  is  similar  for  negatively- 
charged  ionized  acceptors  and  positively-charged  interstitial  defects.  Tan  et  ah  recently 
published  a  review  in  which  they  claim  that  Al-Ga  interdiffusion  is  controlled  by  triply 
negative  Group-m  vacancies,  V  f,;" ,  under  intrinsic  and  n-type  conditions,  while  the 
doubly-positive  Group-III  interstitial,  Ifu  ,  dominates  for  heavy  p-type  doping.P^i 
Baroff  and  Schluter,  however,  predicted  that  the  arsenic  antisite  defect  As ,  a  stable 
donor  in  p-type  GaAs,  was  the  dominant  defect.  In  the  same  review  by  Tan  et  ah , 
data  from  several  authors  were  compiled  which  resulted  in  the  intrinsic  gallium 
self-diffusion  coefficient  of 

(2.2) 

where  D„  =  2.9y.  10®  ernes  '  and  the  activation  energy  Q  =  6  eV.  Earlier  work  by 

Goldstein  determined  that  the  coefficients  for  self-diffusion  =10^  cm^s  *,  and 

Q  =  5.6  eV  for  temperatures  between  1100  and  1200°C.i^J  The  Ga-Al  interdiffusion 
coefficient  is  a  function  of  aluminum  composition,  as  was  shown  in  the  early  work  of 
Chang  and  Koma,(^^]  but  most  authors,  including  Tan  et  ah  take  the  gallium 
self-diffusion  coefficient  as  the  Ga-Al  interdiffusion  coefficient. 

The  dominant  diffusion  mechanism  near  the  surface  (within  a  few  microns)  can  be 
controlled  to  some  extent  by  controlling  the  conditions  at  the  surface  of  the 
semiconductor.  This  is  done  by  annealing  the  semiconductor  in  an  environment  with  an 
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overpressure  of  gallium  or  arsenic.  Hsieh  et  al.  showed  that  the  interdiffusion 
coefficient  was  enhanced  when  annealed  in  an  arsenic  overpressure,  while  interdiffusion 
was  suppressed  when  annealed  in  a  gallium  overpressure.  021  For  the  case  of  arsenic 
overpressure,  where  the  surface  stoichiometry  is  preserved  during  the  annealing 
process,  the  parameters  that  fit  Equation  (2.2)  were  £>0  =  4.4  x  10’  cm^-s  *,  and 
Q  =  6.2  eV.  This  data  is  shown  in  Figure  2.2. 

A  unique  method  of  enhancing  the  interdiffusion  coefficient  was  reported  by 
Hwang  et  c/.P’i  in  which  they  grew  a  GaAs  buffer  layer  by  molecular  beam  epitaxy  at 
low  substrate  temperature  (200 "C).  This  layer  is  highly  resistive  and  optically  inactive 
due  to  the  incorporation  of  excess  arsenic,  poj  The  arsenic  then  acts  as  an  effective 
arsenic  overpressure  which  enhances  the  interdiffusion  process. 

In  addition  to  the  presence  of  the  overpressure,  the  amount  of  overpressure  is  also 
important.  Guido  et  al.  measured  the  interdiffusion  coefficient  for  As4-overpressures 
ranging  from  0.1  atm  to  10  atm.n^J  Hsieh  performed  a  similar  experiment  in  which  the 
gallium  source  weight  was  varied  (the  actual  pressure  is  difficult  to  determine  since  part 
of  the  source  remains  liquid),  and  he  showol  that  a  minimum  occurs  in  the 
interdiffusion  coefficient  for  a  10  mg  source  in  a  3-  to  4-cm^  volume  at  914‘’C.P‘J  Both 
curves  are  reproduced  in  Figure  2.3. 

The  rate  of  interdiffusion  of  heterostructures  can  also  be  controlled  to  some 
degree  by  passivation  of  the  surface  with  a  dielectric  such  as  SiN,  or  SiOj.  Hsieh 
performed  a  diffusion  study  which  indicated  that  the  interdiffusion  for  temperatures 
below  850®C  was  effectively  stopped  when  the  wafer  was  encapsulated  with  a  SiN, 
layer  1000  A  thick.  P‘i  This  method  preserves  the  stoichiometry  of  the  surface,  which 
limits  the  generation  of  excess  point  defects  (gallium  vacancies)  in  the  crystal. 
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Figure  2.2  Al-Ga  interdiffusion  coefficient  as  a  function  of  annealing 
environment.  The  data  is  keyed  as  follows:  square,  arsenic  overpressure;  circle, 
no  overpressure;  diamond,  gallium  overpressure.Pn 
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Figure  2.3  Al-Ga  interdiffusion  coefficient  as  a  function  of  overpressure  source 
quantity.  For  arsenic  overpressure  (a),  the  interdiffusion  coefficient  is  minimized 
near  1  atm  (solid  circles), while  a  minimum  for  gallium  overpressure  (b)  occurs 
for  approximately  10  mg  of  gallium  in  a  sealed  quartz  ampoule  with  a  volume  of 
3-4  cm3.  Gallium  weight  is  used  in  (b)  since  all  the  gallium  does  not  vaporize,  so 
the  exact  overpressure  is  unknown. P'J 
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Deppe  et  al.  examined  both  SiN,  and  SiOj  as  encapsulants  and  found  that  while 
the  SiN,  was  a  good  mask  for  the  introduction  of  defects  at  the  surface,  SiOj  was  a 
source  of  defects  (VqJ  since  the  SiOj  absorbed  gallium  from  the  crystal.P^i  Select  areas 
of  a  wafer  containing  a  superlattice  were  covered  with  SiN,,  then  SiOj  was  deposited  on 
the  entire  surface.  After  annealing,  they  found  that  the  SiNj,-covered  superlattices  were 
preserved,  while  the  SiOj-covered  areas  were  intermixed.  Guido  et  performed  a 
similar  study  with  SiOj  as  an  encapsulant.  They  found  that  the  dielectric  acted  as  a 
source  for  silicon  impurities  if  the  heterostructure  was  too  close  to  the  surface  of  the 
wafer,  thereby  increasing  the  interdiffusion  rate.  Hsieh  suspected  the  same  problem 
from  SiN,  at  temperatures  above  850° C,  but  didn't  rule  out  the  possibility  of  strain 
effects  which  arise  from  the  different  thermal  expansion  coefficients  of  the  substrate 
and  the  encapsulant.  In  the  paper  by  Guido  et  al.  a  stress  field  at  the  edge  of  the 
dielectric  mask  was  determined  to  be  the  cause  of  enhanced  interdiffusion.  Kasahara  et 

also  showed  that  the  rate  of  silicon  diffusion  into  GaAs  was  dependent  on  strain 
due  to  the  thickness  of  the  encapsulant. 

Controlling  the  interdiffusion  rate  via  a  technique  which  generates  defects  at  the 
surface  raises  the  question  of  how  fast  those  defects  are  transported  to  the  layered 
structure.  Guido  et  al.  showed  that  the  interdiffusion  rate  of  a  heterostructure  is 
dependent  on  the  depth  of  a  structure  from  the  surface.  The  structures  were  annealed 
in  2  atm  As^  overpressure  at  825°C,  and  the  diffusion  coefficient  for  theKi;,  was 
D  =  5  X  1 0 ■  cm^s-i.  A  similar  result  was  reported  by  Koteles  et  al.  with  samples 
that  were  capped  with  SiOj  and  rapid-thermal  annealed  in  a  nitrogen  atmosphere  at  800 
to  900°C  for  15  seconds. P*i  They  showed  that  a  coupled  quantum  well  exhibited 
asymmetric  interdiffusion  after  the  annealing  step,  and  proposed  that  it  was  due  to  the 
unidirectional  diffusion  of  vacancies  generated  near  the  surface. 
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One  way  to  control  the  location  from  which  the  defects  are  generated  is  to  implant 
constituent  atoms  into  the  semiconductor  to  generate  internal  "overpressures,"  By 
using  the  same  atoms  as  the  host  for  implantation,  impurity  effects  in  the  final  device 
are  avoided.  The  intermixing  mechanisms  depend  on  defects  from  both  excess  atoms 
(interstitials)  and  damage  from  collisions.  Hirayama  et  performed  a  study  of 
gallium  implantation  into  GaAs/AlGaAs  in  which  they  found  the  intermixing  to  be  a 
function  of  ion  dose.  The  structures  studied  were  within  60  nm  of  the  surface,  and 
implant  energy  was  100  keV.  This  energy  corresponded  to  about  35  nm  penetration 
depth,  so  the  ions  were  implanted  into  the  structures.  For  low  doses  (less  than 
10*^  cm-2),  both  implant  damage  and  stoichiometric  disturbances  were  important,  while 
implantation  damage  dominated  for  higher  doses.  The  ion  dose  used  by  Hirayama  was 
from  10*2  to  10**  cm-^,  most  of  which  was  in  the  quantum  well  structure,  making  the 
excess  gallium  concentration  on  the  order  of  KP®  cm At  that  level,  the  Group-Ill 
interstitials  (primarily  and  Ga^,  antisite  defects  are  the  most  probable  sources  for 
enhanced  intermixing. 

Cibert  et  al.  did  a  similar  study  in  which  Ga+  ions  were  implanted  at  210  keV 
with  doses  of  5  x  1 0  cm-2  to  5  x  1 0  cm-*  which  resulted  in  an  excess  gallium 
concentration  of  6  x  i  o  cm*^  to  6  x  i  o  cm-*,  respectively,  centered  on  the  heteros¬ 
tructure  under  study, Subsequent  rapid  thermal  annealing  (RTA)  steps  at  950°C 
showed  an  increase  in  the  Al-Ga  interdiffusion  coefficient  of  up  to  two  orders  of 
magnitude  compared  to  the  un-implanted  regions  at  the  beginning  of  the  RTA,  After 
the  first  15-second  RTA  step,  the  interdiffusion  coefficient  in  the  implanted  structures 
began  to  decrease  again.  The  implantation  generated  vacancies,  interstitials,  and  defect 
complexes  which  were  annealed  out,  thereby  reducing  the  defects  responsible  for  the 
enhanced  interdiffusion. 


-14- 


While  the  previous  examples  of  ion  implantation  deal  with  ions  (and  hence 
damage)  placed  at  the  layer  interfaces,  Elman  et  al.  studied  the  effect  of  ions  implanted 
away  from  the  structure,  In  this  experiment,  the  layers  were  buried  approximately 
3(X)0  A  from  the  surface  and  low  energy  As+  ions  were  implanted  approximately  160  A 
from  the  surface  using  a  low  implant  energy  of  35  keV  in  doses  of 5  x  i  o  cm-^  to 
5x10*®  cm-2.  Ion  doses  of  2  x  10*“*  cm-^  produced  the  largest  Al-Ga  interdiffusion 
coefficient,  which  was  attributed  to  the  generation  of  vacancies  which  diffused  down  to 
the  layers.  Above  that  dose,  defect  complexes  were  generated  in  the  implanted  region. 
These  complexes  had  low  mobilities  compared  to  point  defects,  so  they  didn't  reach  the 
layers  fast  enough,  and  the  interdiffusion  coefficient  was  gradually  reduced  as  the  dose 
was  increased. 

In  order  to  create  useful  devices  for  integrated  optics,  a  method  of  selectively 
disordering  regions  of  quantum  semiconductor  structures  without  creating  light¬ 
scattering  or  absorption  centers  must  be  developed.  The  underlying  task  in  selective 
disordering  is  to  locally  control  the  concentration  of  point  defects  around  the  quantum 
well.  There  are  several  ways  to  do  this.  The  first  method  is  impurity-enhanced 
diffusion,  either  by  thermal  diffusion^'’*!  or  ion  implantationP’>^i  of  impurities  with  the 
regions  to  be  preserved  masked.  In  impurity-induced  disordered  structures,  however, 
the  optical  absorption  due  to  impurities  is  high  compared  to  undoped  heterostructures. 
Although  this  is  not  a  problem  for  lasers,  which  are  lossy  by  their  very  nature,  losses 
must  be  reduced  in  passive  waveguides  in  order  to  make  integrated  optics  a  viable 
technology.  For  this  reason,  the  use  of  host-type  atoms  to  control  intermixing  such  as 
applying  local  gallium  or  arsenic  overpressure  needs  further  study.  By  using  these 
constituent  atoms  as  an  overpressure  or  implant  source,  impurity  effects  in  the  final 
device  are  avoided. 
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Selective  layer  intermixing  (SLI)  using  IILD  is  well  documented  (see  the  review 
by  Burnham  et  al.  for  instance),  but  work  on  impurity-free  SLI  is  only  beginning. 

In  early  work,  Dq)pe  et  a/.P^i  demonstrated  a  stripe-geometry  laser  with  active  regions 
laterally  defined  by  defect  diffusion.  In  that  case,  ll-?4m  SiN,^  stripes  were  used  to 
protect  the  surface  from  SiOj  which  was  deposited  on  the  surface.  During  the 
annealing  process,  gallium  vacancies  were  generated  at  the  GaAs-SiOj  interface,  which 
diffused  down  to  a  superlattice.  These  regions  were  intermixed  more  than  under  the 
stripes,  resulting  in  a  higher  bandgap  than  under  the  SiN,,  thereby  confining  the 
carriers  and,  to  some  extent,  the  optical  field  in  the  active  region. 

Gallium  or  arsenic  ions  can  also  be  used  to  delineate  structures  either  by  masking 
and  implanting  ions  into  the  surface,  or  by  direct-writing  of  the  ions  using  a  focused  ion 
beam.  Hirayama  et  al.  claimed  to  have  fabricated  selectively  disordered  quantum  wells 
with  lateral  dimensions  of  the  preserved  strips  small  enough  to  quantize  the  carriers  in  a 
second  dimension,  thus  creating  a  quantum  wire.i'*^]  They  used  a  1(X)  keV  gallium 
focused  ion  beam  to  write  1(XX)  A  lines  separated  by  20(X)  A  using  doses  from 
1  X  1 0  cm-2  to  2  X  1 0  cm-2.  A  thermal  annealing  stq)  at  750  to  800'’C  for 
60  minutes  resulted  in  a  lateral  spreading  of  the  disordered  region  of  several  hundred 
angstroms  in  each  direction.  Photoluminescence  excitation  was  used  to  identify  the 
quantized  states  resulting  from  the  lateral  confinement. 

2.1.2  Intermixed  Heterostructure  Potential  Wells 

The  shape  of  the  potential  of  a  heterostructure  is  defined  by  the  bandgaps  of  the 
materials  comprising  its  layers.  For  instance  a  quantum  well  defined  by  a  thin  GaAs 
layer  surrounded  by  AlAs  barriers  has  a  potential  well  in  the  conduction  band  which 
has  a  depth  defined  by  the  difference  between  thef-bandgaps  of  the  GaAs  and  AlAs, 
scaled  by  the  band  offset  between  the  two  materials.  The  valence  band  also  has  a 
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potential  well  for  the  holes  with  a  depth  given  by  the  bandgap  difference  scaled  by  the 
valence  band  offset.  When  the  interfaces  of  these  materials  are  intermixed  by  the 
processes  described  above,  they  become  graded  in  composition,  which  changes  the 
abrupt,  square  quantum  well  to  a  well  with  a  graded  potential  as  depicted  in  Figure  2.4. 
The  shape  of  this  modified  well  depends  on  the  diffusion  coefficient,  which,  in  general, 
is  a  function  of  both  the  temperature  and  the  composition  of  the  diffusing  atoms. 

The  concentration  profile  of  the  intermixed  quantum  well  can  be  found  by  solving 
Pick's  Second  Law,  which  is  given  byt**^! 


SCjzJ) 

dt 


a 

dz 


D{C,T) 


az  ) 


(2.3) 


where  C  ( 2 ,  t )  is  the  composition  of  aluminum  in  the  2:  direction  (perpendicular  to  the 
layers),  t  is  time,  Z)  (  C ,  T )  is  the  diffusion  coefficient  and  T  is  the  diffusion 
temperature.  In  general,  the  derivative  of  D  must  be  taken  since  £>  is  related  to  2: 
through  the  concentration.  If  this  dependence  is  absent,  then  Equation  (2.3)  becomes 


3C(2) 

dt  az"  ■  (2.4) 

The  solution  to  Equation  (2.4)  for  a  single  interface  between  two  semi-infinite  layers  is 
the  complementary  error  function  (er /c  ),  which  is  given  by 


C(z.n-^C,erfc^^^ 


(2.5) 


where  C  o  is  the  initial  aluminum  composition  (or  the  difference  in  composition  if  the 


two  layers  have  different  aluminum  compositions). 

For  a  quantum  well  with  two  interfaces  close  enough  for  their  composition 
profiles  to  cross,  the  results  from  the  two  semi-infinite  layers  are  added  with  the  origins 
of  each  case  displaced  by  half  the  width  of  the  well,  L  The  composition  across  a  well 
centered  at  z  =  0  is  then 
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Figure  2.4  Energy  diagram  of  an  intermixed  GaAs/AlGaAs  quantum  well.  The  dotted 
lines  represent  the  as-grown  potential  well  with  e,-hh,  transition  labeled  as  The 
solid  line  is  the  intermixed  well  with  the  new  e,-hh,  transition,  Eq^'. 
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C{z,t)  =  C, 


(2.6) 


1  --er/c 


(— 


For  a  superlattice,  C(  z ,  t)  is  a  superposition  of  all  the  interfaces  in  the  structure. 

While  the  diffusion  coefficient  is  taken  by  many  researchers  as  a  constant  with 
respect  to  aluminum  composition,  Chang  and  Koma  showed  that  there  is  such  a 
dependence,  and  their  data  is  reproduced  in  Figure  2,5. p*i  If  the  appropriate  conditions 
are  met,  the  Boltzmann-Matano  method  can  be  used  to  find D( Once  this  known, 
the  composition  profile  can  be  found  using  a  first-order  Taylor  series  expansioni**i 
where 


+  =  +  — [  D(2) 


aC(2.to) 

dz 


At 


(2.7) 


This  is  solved  numerically  for  each  position  of  2,  then  incremented  by  A  t  until  the  total 
diffusion  time  is  reached. 

With  the  composition  profile  determined,  the  potential  well  profile  is  obtained  by 
dividing  the  well  into  many  layers,  or  regions,  each  of  which  has  a  bandgap 
characteristic  of  its  average  composition.  The  bandgaps  are  scaled  by  the  appropriate 
band  offset  to  yield  a  good  approximation  of  the  intermixed  well  potential.  The 
remaining  task  is  to  determine  the  electronic  properties  of  the  new  quantum  well. 


Quantum  Wells 


Solutions  of  bound  states  in  finite,  square  quantum  wells  require  either  graphical 
or  numerical  techniques.  One  such  method  was  developed  by  Kolbas  and  Holonyak  in 
which  energies  are  found  by  numerically  solving  the  Schrodinger  Wave  Equation  using 
the  boundary  conditions  of  the  wavefunction  and  its  derivative  at  the  interfaces,  This 
method  also  provided  information  about  the  carrier  envelope  wavefunctions,  and  could 
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Figure  2.5  Al-Ga  interdiffusion  coefficient  as  a  function  of  temperature  and  aluminum 
composition.  The  solid  lines  indicate  the  range  of  the  interdiffusion  coefficient  between 
GaAs  (top  line)  and  AlAs  (^ttom  line).  Each  range  at  a  specific  temperature  is 
indica^  by  a  bar,  and  the  interdiffusion  coefficient  as  a  function  of  aJuminum 
composition  (upper  abscissa)  is  shown  by  the  dotted  lines. 
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be  extended  to  multiple  wells  with  varying  well  depths  and  barrier  heights.  For 
solutions  to  disordered  quantum  wells,  the  potential  well  shape  must  be  approximated 
by  a  finite  number  of  layers  with  constant  potentials. 

To  calculate  the  bound-state  energy  analytically,  Camras  et  used  a  modified 
Poshel-Teller  potentialf^*i  of  the  form 


1/(2)- - :X(\-1) - 5 -  . 

2m  cosh^(a2)  ^2  g 

Where  m*  is  the  effective  mass,  his  Planck’s  constant,  and  the  other  variables  are 

defined  below.  Using  the  well  depth  from  the  error  function  profile  as 
1/(2  =  0)  =  -  AF,  Xis  determined  by 

X  =i+  /  1  ,  2m*  A  £ 

2  \  4  (2.9 

The  other  parameter,  a  ,  is  chosen  such  that  the  Poschel-Teller  potential  best  fits 
the  error  function  profile.  The  energy  of  the  n*  bound  state  is  then  given  as 


h^CL^  ■} 

7(^-n)  ;  n<\  . 

2^  (2.10) 
The  advantage  of  this  approximation  is  that  it  yields  an  analytical  solution  for  the 
energies.  The  problem  with  this  approach  is  that  the  Poschel-Teller  potential  is  a  good 
fit  only  for  a  well  that  has  been  intermixed  to  a  great  extent.  For  a  well  that  has  only 
been  slightly  intermixed,  a  more  precise  method  is  required. 

In  previous  work,  a  "shooting"  methodi^’i  was  employed  to  calculate  the  bound 
state  of  an  error  function  quantum  well.P’i  Hsieh  used  this  model  where  the 
wavefunction  corresponding  to  a  particular  energy  was  found  by  numerically 
integrating  the  Schrodinger  Wave  Equation.  The  energy  was  adjusted  until  the 
wavefunction  generated  using  the  boundary  conditions  of  the  left-side  of  the  well 
matched  the  wavefunction  using  the  right-side  boundary  conditions,  resulting  in  an 
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energy  eigen  state.  In  this  model,  the  wave  equation  includes  a  factor  which  is  a 
measure  of  the  classical  nature  of  the  square-well  potential,  and  Hsieh  assumed  that 
the  square-well  factor  was  a  good  enough  approximation  for  the  error  function  profile. 
This  factor  was  a  function  of  well  width,  so  as  a  well  becomes  more  disordered,  the 
well  width  is  modified  and  this  factor  should  be  adjusted  accordingly. 

Another  problem  with  the  shooting  method  was  discovered  with  ultrathin  wells. 
Since  the  wavefunction  of  such  a  well  has  a  large  spatial  extent,  and  this  entire  space  is 
numerically  evaluated,  the  actual  well  region  accounts  for  only  a  small  part  of  the  total 
space.  In  order  to  keep  the  total  number  of  divisions,  and  hence  the  computation  time, 
reasonable,  few  divisions  are  left  in  the  actual  well  which  reduces  the  accuracy  of  the 
model. 

SupeHattices 

For  superlattices  the  Kolbas-Holonyak  approach  is  well-suited  for  a  finite  number 
of  coupled  square  wells  which  have  few  divisions  per  period.  To  handle  an  intermixed 
superlattice,  each  period  would  have  to  be  divided  into  many  regions.  This  further 
limits  the  number  of  coupled  wells  the  Kolbas-Holonyak  model  can  accommodate  due 
to  practical  computation  considerations.  Therefore,  an  alternate  approach  is  required 
using  a  periodic  potential.  Two  methods  which  have  been  used  for  square- well 
superlattices  are  the  tight-binding  approximation  and  the  Kronig-Penney  model. 

The  tight-binding  calculations  for  superlattices  in  the  envelope-function 
approximation  have  been  limited  to  a  simple  square  well  basis. A  periodic 
potential,  t/  ( z ),  is  used  where  U  =  U  ra)  with  a  being  the  length  of  the 
period  and  r  an  integer. 

With  the  single  well  potential  and  bound  state  energy  given  by  ( z  )  and  E  „, 
respectively,  the  dispersion  relation  is  given  by 
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single  well  calculation  described  earlier,  and  q  is  the  superlattice  wavevector.  With 
nearest  neighbor  interactions  only  (r  =  0 .  ±  1),  this  reduces  to 

£■  =  £„-a-2Pcos(ga)  ,  (2.12) 

where  a  describes  the  integral  in  Equation  (2.11)  for  r  =  0,  and  is  the  shift  to  £„  due  to 

the  coulombic  interaction  of  the  coupled  wells.  The  integral  forr  =  ±  1  is  described  by 
|3,  which  is  the  exchange  energy  that  accounts  for  the  broadening  of  the  bound  state  into 
an  energy  band. 

While  the  envelope  function  approximation  is  valid  as  long  as  the  spatial  extent  of 
the  wavefunction  is  large  compared  to  the  atomic  pacing,  it  becomes  invalid  when  a 
superlattice  period  is  only  a  few  atomic  layers.  For  thin  single  wells,  the  wavefunction 
extends  far  into  the  barriers,  making  it  large  overall.  But  in  a  short-period  superlattice, 
the  wavefunction  is  periodic  with  the  same  period  as  the  structure.  In  this  case,  an 
atomic  tight-binding  calculation  is  more  appropriate,  such  as  that  done  by  Rucker  et 

Another  method  to  find  the  energy  dispersion  relation  for  a  square-well 
superlattice  is  the  Kronig-Penney  model,  In  the  effective  mass  approximation  this 
model  gives  an  exact  solution  to  the  periodic  square-well  potential,f”i  but  for  more  than 
two  layers,  the  usual  method  of  solution  becomes  impractical.  An  approach  described 
by  Yuh  and  Wangi*«i  defines  a  2  x  2  transfer  matrix  for  each  layer  in  the  basis,  given  by 

m* 

cosC/Cjli)  sin(kjl,) 

ki 

- :sin(/c,f,)  cos(/CiZ.<) 

rrii 
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where  £  denotes  the  layer,  If  is  the  layer  thickness,  and 


(2.14) 


The  transfer  matrices  of  each  region  are  multiplied  together  to  give  the  total  matrix,  S. 
Using  r(2),  a  vector  with  elements  Ti  (2)  =  Tp(2)  andr2(2)  -  the 

eigenvalues  of  S  are  found  using 

r(r  +  a)  =  sr(2)  =  \ir(2)  .  (2.15) 

Along  with  the  Bloch  theory,  which  givesX.  =  e**’®,  we  find 


Sii  +  522  =  e“'“  +  e‘'‘'“  =  2cos(ga)  .  (2-16) 

where  q  is  the  superlattice  wavevector.  If  a  simple  two-layer  basis  is  used,  the  result  is 

the  familiar  Kronig-Penney  solution. 


2.2  Optical  Properties  of  Semiconductor  Quantum  Structures 

The  optical  properties  of  a  semiconductor  are  contained  in  the  complex  dielectric 
constant,  €(00),  which  is  given  by 

e(tw) -€  1  (00)  162(00)  .  (2.17) 

Once  either  the  real  or  imaginary  part  is  known,  the  other  part  can  be  found  using  the 
Kramers-Kronig  transformations.  J*’-**! 

The  usual  method  for  determining  the  dielectric  properties  below  the  bandgap  is  to 
measure  the  reflectance  at  the  surface  to  first  obtain  the  index  of  refraction.  This  was 
performed  for  AlGaAs  by  Casey  et  where  the  index  was  found  for  energies  below 
as  well  as  slightly  above  the  bandgap.  The  extinction  coefficient  was  neglected  since  it 
was  still  small  and  only  gave  an  error  of  about  1  %  in  the  worst  case.  Aspnes  et  al. 
measured  the  complex  reflectance  ratio  for  various  angles  using  a  rotation  analyzer 
spectroellipsometer,  from  which  they  calculated  the  real  and  imaginary  parts  of  the 
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dielectric  constant. 

Several  empirically  determined  formulae  exist  for  the  index  of  refraction  of 
AlGaAs  below  the  bandgap,i**-“i  one  of  which  (given  by  Hunsperger)  is 

-D(.x)Kl  . 

Ko  ^2.18) 

where  AT  is  the  aluminum  fraction  in  the  Al,Ga,.,As  alloy,  is  the  free-space 
wavelength  in  microns,  and 

A=  10.9060- 02. 92>: 

5  =  0.97501 

C  =  (0.52886-0. 735x)^  :x<0.36 
C  =  (0.30386 -0.105jc)^  ;jc>0.36 
Z)  =  0.002467(l  +  1.41X)  . 

Recently,  Jensen  and  Jensen  presented  a  model  in  which  the  real  index  near  the  band 
edge  can  be  calculated  directly  from  experimentally  available  band  structure 
parameters.!^! 

To  calculate  the  dielectric  function,  one  usually  starts  with  the  imaginary  part, 

€ ' '  ( uo ),  which  is  proportional  to  the  square  of  the  optical  matrix  elemem  and  the  joint 
density  of  the  states  for  the  electronic  transitions. Once  6  "  is  calculated,  g  '  is 
found  using  the  Kramers-Kronig  transformations,  then  the  index  of  refraction, 
extinction  coefficient,  or  absorption  coefficient  may  be  determined. 

An  alternate  method  for  calculating  the  dielectric  constant  is  to  treat  the  electronic 
system  as  a  collection  of  oscillators  and  add  the  contributions  of  each.  In  this  manner, 
both  G  'and  e  "  are  calculated  separately  since  the  energy  spectrum  is  too  narrow  to 
accurately  apply  the  Kramers-Kronig  relations.  Kahen  and  Leburton  used  such  an 
approach!2°!  where  they  calculated  the  band  structure  using  the  Kane  model,!^^  then  used 
the  calculated  properties  of  the  band  edges  for  the  characteristics  of  the  oscillators. 
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This  method  was  later  applied  to  superlattices  and  multiple  quantum  wells  in  which  the 
energy  dispersion  in  the  quantized  direction  was  given  by  a  tight-binding-like 
expression.  The  oscillator  approach  is  well  suited  for  calculating  the  dielectric 
constants  of  quantum  wells  as  performed  by  Tsu  and  Ioriattii“J  since  their  energy  levels 
are  discrete. 

In  the  oscillator  approach,  the  complex  dielectric  constant  for  a  system  of  N  is 


:(Go)  =  eo  +  7V  V  - - - p- - 

Y^o(oof-<Ai2)-i2Yy 


CO 


(2.19) 


The  dielectric  function  is  then  the  sum  of  the  dielectric  responses  of  all  the  oscillators  in 
the  atomic  system,  each  with  a  resonant  frequency,  co  j ,  oscillator  strength  /  j ,  and 
damping  coefficient  y  j.  The  oscillator  strength  is 

- |M|" 

^2.20) 

where  |  Af  |  ^  is  the  momentum  matrix  element  of  the  oscillator  squared. 


The  properties  of  the  quantum  well  are  then  dependent  on  the  oscillator 
frequency,  its  damping  coefficient  and  the  matrix  element.  The  frequency  can  be 
calculated  by  solving  for  the  energy  levels  in  the  quantum  well.  This  should  also  yield 
the  envelope  wavefunctions  of  the  electron  and  hole  states,  which  is  useful  for 
calculating  the  matrix  element. 

The  matrix  element  for  a  single  oscillator  between  two  states,  m  and  n,  with  a 
single  polarization  is 

where  <  n  |  m  >  is  the  overlap  of  the  electron  wavefunction  in  state  n  and  hole 
wavefunction  in  state  m.  /’mn  is  the  interband  matrix  element  including  the  selection 
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rules,  and  the  delta  function  conserves  momentum  in  the  directions  parallel  to  the  layers 
in  the  quantum  well.  The  inteiband  matrix  element  in  m-V  semiconductors  for 
electron  to  heavy-hole  (TE,  transverse  electric)  transitions  at  2  *  Oist®®i 


where  the  factor  of  \  comes  from  5  for  the  average  over  all  directions  ^  for  bulk)  and  \ 

for  the  conservation  of  spint^^i  (since  the  transition  can  only  occur  to  the  band  with  the 
same  spin).  P  0  is  found  usingi*^^ 

2  _  ^0^  of  ^  ^  E  g-*-  A 

2  U:  Je.^Ia  (2.23) 

where  mg  is  the  ffee-electron  mass,  is  the  bandgap  energy,  m*  is  the  conduction- 

band  effective  mass,  and  A  is  the  split-off  energy  for  the  split-off  valence  band.  For 
GaAs  at  room  temperature,  the  squared  matrix  element  is  estimated  by  Casey  and 
PanishP‘1  as 

1.33m,f,  .  2.24) 

the  quantum  well  value  being  3/2  |  Mbu,*  |^.  The  value  is  generally  considered  to  be 

constant  between  the  materials  of  a  heterostructure. 

The  value  for  the  damping  coefficient  can  be  measured  from  luminescence  data,^i 
where  2 y  =  Aou  is  the  full  width  at  half  maximum  (FWHM)  of  the  luminescence 
spectrum.  The  line- width  broadening  is  due  to  the  finite  lifetimes  of  the  electronic 
states  and  is  related  to  the  lifetime  through  the  uncertainty  relation  ft  A  cox  •  ft,P^i 
where  x  is  the  lifetime  of  the  excited  state.  The  lifetime  can  be  either  measured  or 
calculated,  but  if  calculated,  both  radiative  and  non-radiative  processes  should  be 
considered. 
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3  Experimental  Methods 


3.1  Crystal  Growth 

The  structures  described  in  this  work  were  grown  by  molecular  beam  epitaxy 
(MBE).  The  MBE  system  used  is  a  commercial  Varian  360.  Since  the  emphasis  of  this 
study  is  to  ascertain  the  effects  of  thermal  processing  on  the  optical  properties  of 
quantum  wells,  the  details  of  crystal  growth  are  beyond  the  scope  of  this  thesis.  The 
interested  reader  is  referred  to  an  excellent  description  of  the  system  by  Andersen, 
but  the  sample  preparation  and  growth  conditions  of  the  samples  grown  for  this  work 
are  specified  below. 

The  samples  are  prepared  by  degreasing,  etching,  and  mounting  them  on 
molybdenum  blocks.  Unless  otherwise  stated,  all  substrates  are  semi-insulating  (SI), 
(100)-oriented  GaAs  wafers.  The  substrates  are  degreased  by  rinsing  for  five  minutes 
each  in  trichloroethane  (TCE),  acetone,  and  methanol  in  order,  with  a  final  rinse  in 
deionized  water  (DI-HjO).  The  surface  is  then  etched  with  a  solution  of  7: 1 : 1 
H2S04:H202:H20.  The  substrates  are  then  mounted  on  the  molybdenum  blocks  using 
molten  indium,  and  immediately  loaded  into  the  preparation  chamber  of  the  MBE  via  a 
load-lock.  The  samples  are  then  heated  to  outgas  the  surface,  and  just  before 
beginning  the  growth  of  the  epitaxial  layers,  the  wafers  are  heated  to  575  in  the 
growth  chamber  for  five  minutes  under  an  arsenic  overpressure  to  remove  any  oxide 
that  may  have  formed  on  the  surface. 

The  growth  conditions  consist  of  a  combination  of  molecular  beam  fluxes,  which 
are  measured  before  and  after  layer  growth,  and  substrate  temperature,  which  is 
monitored  during  growth.  The  growth  rates  and  alloy  compositions  are  calibrated  using 
calibration  samples  just  prior  to  the  growth  of  the  desired  structures.  The  crystal 
quality  is  determined  using  reflection-electron  diffraction  and  is  verified,  along  with 
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material  composition,  using  photoluminescence  spectroscopy.  The  substrate 
temperature  is  200®C  for  low-temperature  growth  of  GaAs.  For  the  Al,Gai.^s  the 
temperature  ranged  from  640°C  for  GaAs  to  700®C  for  ALAs,  and  the  Ino  jGa^  gAs  is 
grown  at  520°C. 

3.2  Optical  Characterization 

The  optical  characterization  performed  in  this  study  was  photoluminescence 
spectroscopy.  Two  techniques  used  were  standard  (time-integrated)  photoluminescence 
(PL)  in  which  no  temporal  information  was  collected,  and  time-resolved  photolumines¬ 
cence  (TPL)  in  which  temporal  data  was  taken  with  a  resolution  of  about  IS  ps.  Both 
PL  and  TPL  were  used  to  examine  the  characteristics  of  the  spontaneous  emission  of 
optically  excited  samples  in  order  to  probe  their  electronic  structures. 

The  PL  process  shown  in  Figure  3.1  begins  with  the  creation  of  electron-hole 
pairs  in  a  semiconductor  via  the  absorption  of  photons  of  energy  higher  than  that  of  the 
electronic  bandgap  of  the  semiconductor.  This  process  is  refened  to  as  photopumping 
the  sample.  The  scattering  times  of  the  electrons  (holes)  within  the  conduction 
(valence)  band  are  much  shorter  than  the  lifetime  for  the  electrons  and  holes  to 
recombine,  so  they  quickly  scatter  to  the  conduction  (valence)  band  edge  before 
recombining.  Momentum  and  energy  is  conserved  as  the  carriers  thermalize  to  the 
band  edges  through  the  emission  of  phonons.  The  carriers  then  recombine  across  the 
bandgap,  emitting  photons  equal  in  energy  to  that  of  the  gap,  thereby  conserving 
energy.  If  these  photons  are  not  re-absorbed,  then  they  will  escape  from  the  sample, 
where  they  can  be  collected  and  detected.  In  the  case  of  TPL,  the  excitation  pulsewidth 
must  also  be  much  shorter  than  the  recombination  time  in  order  to  observe  the  temporal 
characteristics  of  the  electronic  system. 
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Figure  3.1  The  photoluminescence  process.  The  sample  is  photopumped  with 
photons  of  energy  h  v  p  (1)  creating  an  electron-hole  pair.  The  carriers  then 
thermalize  to  the  band  edges  (2)  where  some  recombine,  giving  off  photons  of 
energy  /i  v  (3)  while  others  diffuse  toward  the  well  (4),  The  carriers  are  then 
captured  by  the  well  (5)  where  they  thermalize  to  the  bound  states  hrom  which 
they  recombine  to  give  off  photons  of  energy  v  (6). 
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The  detection  apparatus  measures  the  intensity  of  the  emitted  photons  as  a 
function  of  wavelength,  which  is  easily  converted  to  energy.  In  the  case  of  TPL,  the 
temporal  characteristics  of  the  photon  emission  process  is  measured  as  well.  The  PL 
spectrum  contains  information  about  the  electronic  band  structure  of  the  structure,  and 
the  time-resolved  data  from  TPL  adds  to  it  information  about  the  carrier  recombination 
lifetimes. 

3.2.1  Standard  Photoluminescence  (PL)  Spectroscopy 

The  standard  (or  time-integrated)  PL  system  used  a  Coherent  Inova  90  argon-ion 
laser  with  a  Coherent  7208  cavity  dumper  to  pump  the  samples.  The  laser  is  capable  of 
5  W  of  continuous  wave  (CW)  laser  power  at  the  wavelength  X.  =  5145  A,  and  is 
capable  of  much  higher  peak  power  when  used  in  pulsed  mode.  The  pulse  rate  is 
selected  using  a  Model  7200  Cavity  Dumper  Driver,  and  the  pulse  width  is 
approximately  4  ns  full-width-at-half-maximum  (FWHM).  The  lower  power  CW 
operation  allows  the  photopumped  electron-hole  pairs  in  the  sample  to  thermalize  and 
recombine  before  the  energy  band  fills  up.  Conversely,  the  pulsed  mode  concentrates 
more  photons  on  the  sample  during  the  pulse  so  that  radiative  recombination  from 
higher  energy  features  can  be  detected,  such  as  the  barrier  energy  in  a  quantum  well 
structure. 

The  samples  are  cleaved  pieces  of  the  original  wafer,  and  measure  from  a  few 
millimeters  to  as  much  as  a  centimeter  on  a  side,  and  are  mounted  on  a  heat-sinking 
"plug.”  The  plug  is  a  solid  copper  cylinder  with  an  indium  cap  soldered  on  top  of  it 
which  is  pressed  flat  to  accept  the  sample.  Tbe  sample  is  either  pressed  into  the  indium 
or  is  held  in  contact  to  the  indium  with  radio-TV  service  cement.  The  plug  is  then 
bolted  to  a  copper  block  which  is  attached  to  the  bottom  of  a  stainless  steel  dewar.  The 
dewar  then  fits  on  top  of  a  lighthouse,  or  a  vacuum  chamber  with  three  quartz  windows 
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which  allow  optical  access  to  the  isolated  sample.  The  lighthouse  is  evacuated  with  a 
mechanical  pump  to  less  than  20  mTorr,  and  liquid  nitrogen  is  poured  into  the  dewar  to 
cool  the  sample  to  77K. 

Light  from  the  laser  is  focused  to  a  spot  size  of  about  20  %m  to  100  %m  on  the 
sample's  surface,  which  is  held  approximately  20°  from  normal  incidence  to  the  pump. 
The  emitted  light  is  collected  90°  for  the  direction  of  the  input  beam  (about  70°  from 
the  sample),  the  pump  light  is  filtered  out,  and  the  remaining  light  is  focused  through 
the  input  slit  of  a  0.5-m  Jerrel  Ash  Model  82-020  spectrometer  with  a  dry-iced-cooled 
S-1  photomultiplier  tube  at  1000  V  to  detect  the  spectrally  dispersed  signal.  The  output 
of  the  photomultiplier  is  amplified  by  a  Keithly  602  electrometer  which  drives  a  chart 
recorder  to  record  the  data. 

3.2.2  Time-Resolved  Photoluminescence  (TPL)  Spectroscopy 

The  basic  apparatus  for  TPL  is  similar  to  that  for  standard  PL  except  for  the 
addition  of  a  streak  camera  to  measure  the  time  response  of  the  optical  emission. 
FoUowing  is  a  brief  description  of  the  TPL  system.  A  more  detailed  description  with 
the  theory  of  operation  is  given  by  Benjamin.P«i 

The  excitation  source  in  this  case  is  a  Spectra  Physics  Model  37SB  dye  laser 
which  emits  tunable  laser  radiation  in  5-  to  6-ps  pulses  at  an  82  MHz  repetition  rate. 
The  wavelength  range  of  the  dye  laser  is  62Q0  to  6900  A  using  DCM  dye,  which  is 
synchronously  pumped  by  a  frequency-doubled  Model  3460  Nd:YAG  laser.  The 
Nd:YAG  is  mode-locked  to  produce  120-ps  pulses  at  1064  nm,  which  is  then 
frequency-doubled  to  85-ps  pulses  at  532  nm  by  passing  through  a  KTP  crystal.  The 
final  output  of  the  laser  system  is  monitored  using  a  Spectra  Physics  Model  409 
autocorrelator,  which  displays  the  pulsewidth  on  an  oscilloscope. 
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Sample  plugs  are  mounted  on  the  cryotip  of  an  Air  Products  closed-cycle  helium 
refrigeration  system  (Model  DE202  expander  and  Displex-type  compressor).  The 
refrigeration  system  is  capable  of  holding  the  sample  plug  at  between  15K  and  300K, 
but  all  data  in  this  study  was  taken  at  77K.  The  excitation  source  is  focused  on  the 
sample  through  a  quartz  window  in  a  vacuum  shroud  that  isolated  the  sample,  which  is 
held  at  approximately  1(H  Torr  via  a  mechanical  pump.  The  sample  luminescence 
emitted  through  a  second  quartz  window  is  collected,  collimated  and  focused  on  a 
spectrometer  slit  using  a  dual-lens  system. 

For  spectral  analysis,  an  Instruments  SA  DHR320S  spectrometer  with  a 
1200-grooves-per-millimeter  diffraction  grating  is  used.  It  employs  a  double-grating 
configuration  to  eliminate  temporal  broadening  due  to  the  optical  path  length  difference 
existing  in  single-grating  spectrometers,  which  can  be  on  the  order  of  the  carrier 
lifetimes  of  the  samples  being  studied. 

The  output  of  the  spectrometer  is  collimated  and  focused  onto  the  input  slit  of  a 
Hamamatsu  C1S87  streak  camera.  The  light  is  then  focused  on  an  S-1  photocathode 
where  it  is  converted  to  an  electron  image.  The  electrons  are  accelerated  by  a  mesh 
electrode  into  a  deflecting  electric  field,  A  high-speed-ramp  voltage,  which  is  triggered 
by  the  excitation  source  via  a  high-speed  detector,  is  applied  to  the  deflecting  plates, 
which  sweeps  the  electron  beam  from  top  to  bottom.  The  electron  image  is  amplified 
by  a  microchannel  plate,  then  strikes  a  phosphor  screen,  which  displays  temporal 
dispersion  vertically,  and  spectral  dispersion  horizontally. 

The  phosphor  image  is  sent  to  a  Hamamatsu  CIOOO  High  Sensitivity  Video 
Camera,  which  passes  its  video  signal  to  a  Hamamatsu  C1098  Temporal  Analyzer. 

The  temporal  analyzer  displays  the  screen  image  on  a  video  monitor,  where  a  vertical 
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window  representing  a  narrow  spectral  band  can  be  selected  for  analysis.  The  intensity 
versus  time  data  centered  at  various  wavelengths  is  then  stored  in  a  Hewlett  Packard 
310  computer  for  analysis  and  presentation. 

The  time  resolution  of  the  streak  camera  is  limited  by  the  time  spread  of 
photoelectrons,  spatial  spread  of  the  input  slit  and  resolution  of  the  readout  system,  and 
the  deflection  electric  field.  The  total  resolution  is  limited  further  by  the  time  spread  of 
the  collection  optics,  and  that  of  the  trigger  jitter  for  low-level  signals  where  multiple 
integrations  are  performed.  The  total  resolution  of  the  detection  system  is 
approximately  15  ps,  which  is  sufficient  for  measuring  the  temporal  responses  of  the 
samples  in  this  study,  which  range  from  about  150  ps  to  over  1000  ps. 

3.3  Thermal  Annealing 

Thermal  annealing  of  the  quantum  well  structures  results  in  intermixing  the 
materials  of  the  heterostructures  at  their  interfaces.  This  process  modifies  the  shape  of 
the  potential  well,  which  in  turn  changes  the  optical  and  electronic  properties  of  the 
structure.  This  phenomenon  is  the  general  focus  of  this  thesis,  and  will  be  addressed  in 
detail  in  later  sections.  The  overall  process  to  anneal  a  sample  is  to  clean  it  and  any 
source  of  overpressure,  seal  it  with  the  overpressure  source  in  an  evacuated  ampoule, 
heat  it  to  the  annealing  temperature  for  a  predetermined  time,  cool  it,  and  remove  it 
from  the  ampoule. 

The  first  step  in  preparing  the  samples  is  to  remove  any  indium  that  is  on  the  back 
side  of  the  wafer.  This  is  done  by  mounting  the  sample  top-down  on  a  pre-cleaned 
glass  slide  using  melted  black  (apiezon)  wax.  Large  pieces  of  indium  are  then  scraped 
off,  and  the  remaining  indium  is  etched  off  in  HCl  for  2-3  hours.  The  slides  are  rinsed, 
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then  the  wax  is  re-melted,  and  the  samples  are  immersed  in  TCA  to  dissolve  the 
remaining  wax.  The  cleaning  step  for  the  samples  consists  of  standard  degreasing  steps 
described  in  the  section  on  Crystal  Growth. 

For  anneals  using  arsenic,  a  solid  piece  of  arsenic  is  loaded  into  a  quartz  tube, 
which  is  evacuated  to  10-*  Torr.  The  arsenic,  which  is  dark  from  arsenic  oxide,  is  then 
heated  in  excess  of  the  sublimation  temperature  using  a  hydrogen-oxygen  torch.  The 
oxide  is  driven  from  the  surface,  along  with  the  top  layers  of  pure  arsenic,  leaving 
behind  pure  solid  arsenic  with  a  shiny,  silver-colored  surface.  The  sublimated  material 
condenses  at  the  opposite  end  of  the  ampoule,  which  is  kept  cool  with  a  towel  which  is 
soaked  with  cool  tap  water.  When  the  sample  and  the  ampoule  to  be  used  in  the  anneal 
is  ready,  the  arsenic  is  removed  from  the  vacuum,  cut  and  weighed,  and  loaded  along 
with  the  prepared  sample  into  the  ampoule  in  which  the  annealing  will  take  place.  This 
way,  the  clean  arsenic  is  exposed  to  the  air  for  a  minimal  amount  of  time. 

For  anneals  using  a  gallium  overpressure,  the  gallium  source  must  be  crystallized 
first,  since  it  is  a  liquid  near  room  temperature  with  a  melting  point  of  29“C.  First,  a 
beaker  of  gallium  (99.9999%)  is  warmed  to  ensure  it  is  a  liquid,  and  the  top  layer  is 
skimmed  to  remove  any  debris  that  may  have  fallen  or  adsorbed  onto  the  surface  during 
storage.  A  pre-cleaned  microscope  slide  is  then  dipped  into  liquid  nitrogen  to  cool  it  to 
well  below  the  melting  point  of  the  gallium.  The  slide  is  then  immersed  in  the  liquid 
gallium,  where  the  gallium  crystallizes  on  the  slide.  If  the  slide  is  removed  from  the 
liquid  just  before  the  crystals  warm  and  start  to  re-melt,  large  gallium  crystals  can  be 
formed,  measuring  up  to  about  5  mm  long  and  weighing  as  much  as  a  gram.  These 
crystals  are  hard  and  brittle,  which  makes  them  easy  to  handle  and  cleave  into  pieces  of 
the  desired  size.  As  with  the  arsenic,  the  gallium  sources  are  loaded  into  an  ampoule 
with  the  sample  and  evacuated  as  soon  as  possible. 
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Both  sources  are  cut  with  clean  stainless  steel  razor  blades,  then  weighed  on  a 
Seederer-Kohlbusch  balance.  If  needed,  the  sources  are  re-cut  and  weighed  several 
times  until  the  required  weight  is  obtained.  For  arsenic,  that  weight  is  determined  by 
the  desired  overpressure  (P),  the  annealing  temperature  (T),  and  volume  of  the  sealed 
ampoule  (V).  Assuming  that  all  the  solid  arsenic  becomes  a  vapor  at  the  annealing 
temperature,  and  at  high  temperatures  (800  to  the  AS4  pressure  is  much  higher 

than  the  ASj  pressure,  then  the  overpressure  is  composed  of  AS4  with  a  molar  weight 
(M)  of  4(74.922)  g-mole  *.  The  weight  (W)  of  the  solid  source  is  obtained  using  the 
Ideal  Gas  Law, 

PV  =  nRT  .  (3.1) 

where  R  is  the  gas  constant  and  n  =  W/M.  Because  the  smallest  tare  weight  available 
in  our  lab  is  5  mg,  the  length  of  the  ampoule  is  easier  to  adjust  than  the  weight  in  order 
to  get  the  desired  overpressure.  In  this  study,  all  anneals  involving  arsenic 
overpressure  were  done  at  850®C  (1123K)  with  1  atm  AS4.  Using  a  15  mg  source  and 
an  ampoule  with  an  inside  diameter  of  7.5  mm,  the  length  of  the  sealed  ampoule  is 
10.4  cm.  This  size  is  small  enough  to  handle  easily,  and  large  enough  to  minimize  the 
error  in  the  assumption  that  the  ampoules  are  cylinders  with  flat  ends. 

For  gallium  overpressure  experiments,  the  vapor  pressure  in  the  ampoule  is 
difficult  to  ascertain,  since  part  of  the  source  remains  a  liquid  at  the  annealing 
temperature.  Since  the  intent  of  such  an  experiment  in  the  context  of  this  thesis  was  to 
minimize  the  interdiffusion  rate  at  the  interfaces,  data  reported  by  HsiehP’i  was  used  to 
determine  the  source  weight  for  the  minimum  interdiffusion  coefficient.  In  his  work, 
Hsieh  found  that  a  source  weight  of  10  mg  in  an  ampoule  with  a  volume  between 
3-4  cm^  yielded  the  minimum  interdiffusion  coefficient  of  about2  x  1  O' cm^s  *.  In 
this  work,  the  volume  used  was  about  1.5  times  that  used  by  Hsieh,  so  a  source  weight 
of  15  mg  was  used  to  yield  the  smallest  interdiffusion  coefficient.  Although  this 
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estimation  is  rather  crude,  halving  or  doubling  the  source  weight  resulted  in  less  than  a 
50%  increase  in  the  interdiffusion  coefficient.  Also,  gallium  overpressure  was  used 
only  in  a  proof-of-concept  experiment  where  order-of-magnitude  differences  between 
gallium-,  arsenic-,  and  no-overpressure  anneals  were  being  sought. 

The  ampoules  into  which  the  samples  and  sources  were  loaded  are  fashioned  from 
9-mm  (outside  diameter)  quartz  tubes.  The  tubes  are  approximately  15  inches  long 
with  one  end  sealed  using  a  hydrogen-oxygen  torch.  The  insides  of  the  tubes  are  etched 
with  1:10  HFiDI-HjO  for  10  minutes,  and  rinsed  with  DI-HjO.  The  tubes  are  dried 
and  stored  in  a  bakeout  oven  at  200®C.  When  needed,  a  tube  is  removed  from  the 
oven,  the  samples  and  source  are  loaded  into  it,  and  it  is  loaded  onto  the  vacuum  station 
to  be  evacuated. 

Rough-pumping  at  the  vacuum  station  is  accomplished  by  a  mechanical  pump 
which  achieves  vacuums  of  about  40  mTorr.  At  that  point,  the  mechanical  pump  is 
isolated  from  the  ampoule,  and  a  diffusion  pump  evacuates  the  ampoule  to 
approximately  10*  Torr,  where  the  ampoule  can  be  sealed. 

The  ampoules  are  sealed  in  several  steps.  First,  any  water  that  has  adsorbed  onto 
the  inside  surface  of  the  ampoule  is  driven  off  using  a  hydrogen  flame.  Next,  if  a 
source  is  being  used,  a  dimple  is  made  with  the  torch  (now  mixed  with  oxygen)  near 
the  end  of  the  tube.  The  dimple  keeps  the  source  and  samples  separated  while  loading 
it  into  the  furnace.  Finally,  the  tube  is  collapsed  at  the  appropriate  length  from  the  end 
until  the  ampoule  is  closed  off  from  the  rest  of  the  tube.  The  torch  is  then  used  to 
separate  the  ampoule  from  the  tube. 

After  carefully  shaking  the  ampoule  so  that  the  source  and  samples  are  separated 
by  the  dimple  and  samples  are  top-side-up,  it  is  loaded  into  a  Lindberg  Hevi-Duty  36" 
3-zone  furnace  which  has  been  set  to  the  proper  temperature.  The  ampoules  are  slid 
into  the  center  zone  where  the  temperature  is  monitored  by  a  chromel-alumel 
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thermocouple  connected  to  an  Omega  Model  660  Thermocouple  Thermometer.  After 
annealing  for  30  minutes  to  S  hours,  the  ampoule  is  removed  from  the  furnace  by 
pulling  it  to  the  edge  and  grasping  it  with  a  wet  towel.  The  wet  towel  quenches  one 
end  of  the  ampoule  (without  the  sample)  onto  which  the  vapors  condense.  The  rest  of 
the  ampoule  is  then  quenched  under  tap  water.  The  ampoules  are  then  cut  and  the 
samples  are  removed  for  characterization. 

3.4  Device  Processing 

Device  processing  for  this  work  consisted  of  fabricating  stripes  layers  of  either 
SiN,  or  GaAs  epitaxial  layers  grown  at  low  substrate  temperature  (LT-GaAs). 

The  SiN,  film  is  deposited  on  GaAs  wafers  by  plasma-enhanced  chemical  vapor 
deposition  using  an  LFE  Corporation  Model  PND-310  plasma  nitride  deposition 
system.  Surface  oxide  is  removed  hrom  the  samples  prior  to  nitride  deposition  by 
etching  in  HCl  for  IS  seconds,  rinsing  wirh  DI-HjO  for  2  minutes,  dipping  in 
NH^OHiDI-HjO  (1:15)  for  20  seconds,  and  then  rinsing  in  DI-HjO  for  another 
2  minutes.  The  samples  are  then  loaded  into  the  PND-301  vacuum  chamber  where  they 
are  heated  to  300®C  in  a  vacuum  of  less  than  20  mTorr.  SiH4  (1.7%,  diluted  with 
argon)  is  introduced  with  ultra  high  purity  Nj  into  the  chamber  for  deposition.  A 
1(X)  W  radio  frequency  (RF,  13.56  MHz)  plasma  is  generated  in  the  chamber  for  5 
minutes.  Under  ideal  conditions,  a  1000  A-thick  Si3N4  film  is  deposited  on  the  surface 
of  the  wafers. 

The  quality  of  the  SiN,  is  checked  using  a  computer-controlled  ellipsometer.  The 
ellipsometer  measures  the  thickness  of  the  film  and  the  index  of  refraction,  which 
should  be  as  close  as  possible  to  2.02.  The  film  is  then  checked  under  an  optical 
microscope  to  ensure  that  the  surface  is  clear  of  defects. 
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Standard  photolithography  is  used  U)  define  stripes  on  both  the  SiN,  Hlms  and 
LT-GaAs.  The  samples  are  pre-baked  at  120®C  for  15  minutes  to  dehydrate  the 
surface.  Each  sample  is  loaded  on  a  Headway  spinner  fitted  with  a  0.5-inch  vacuum 
chuck.  One  drop  of  AZ  1400-27  positive  photoresist  is  deposited  at  the  center  of  the 
sample,  which  is  then  spun  at  4000  rpm  for  40  seconds,  resulting  in  a  photoresist  layer 
1.35  ^m  thick.  The  samples  are  then  soft-baked  at  90"C  for  30  minutes  to  remove 
solvent  from  the  photoresist. 

After  it  is  cool,  each  sample  is  loaded  into  a  Karl  Suss  M4  mask  aligner  with  the 
appropriate  photoresist  mask.  There,  the  photoresist  is  exposed  to  ultraviolet  light 
(310  nm)  for  90  seconds.  After  exposure,  the  sample  is  removed  and  the  photoresist  is 
developed  in  one-part  Shipley  351  developer  diluted  in  five-parts  DI-HjO  for 
24  seconds,  followed  by  a  rinse  in  DI-HjO.  The  patterns  are  inspected  under  an  optical 
microscope  fitted  with  a  yellow  filter  to  check  the  quality  of  the  developing,  and,  if 
satisfactory,  the  photoresist  is  hard-baked  at  120‘*C  for  15  minutes. 

The  exposed  SiN,  stripes  are  etched  in  a  Technics  Model  PEII  reactive  ion  plasma 
etcher.  The  samples  are  loaded  into  the  chamber  which  is  evacuated  to  less  than 
0.4  Torr,  then  CF4  (1.25  Torr)  and  Oj  (0.05  Torr)  are  introduced.  A  1(X)-W  RF 
plasma  is  generated  in  which  the  exposed  SiN^  etches  for  3  minutes  (500-l(XX)  A  per 
minute).  After  etching,  the  photoresist  is  removed  in  acetone  with  ultrasonic  agitation 
for  5  minutes. 

The  stripes  in  LT-GaAs  are  delineated  with  a  selective  wet-chemical  etch  which 
stops  at  AlGaAs.  The  etchant  is  one  drop  of  NH40H:DI-H20  (1:20)  in  HjOjiDI-HjO 
(22: 100).  Small  amounts  of  H202:DI-H20  are  then  added  to  adjust  the  pH  to  7.0.  The 
etch  rate  is  approximately  40  A  per  minute  with  constant  stirring,  and  the  LT-GaAs  is 
completely  removed  when  the  stripes  take  on  a  blue  color  characteristic  of  AlGaAs. 

The  photoresist  is  then  removed  in  acetone  as  above. 
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3.5  Hall  Mobility  Measurements 

The  electron  mobility  in  modulation-doped  field  effect  transistor  (MODFET) 
structures  is  measured  using  the  Vander  Paw-Hall  technique.  The  samples  are  cleaved 
into  1-cm  squares,  and  small  indium  dots  are  placed  on  the  comers  of  the  surface. 

Each  sample  is  then  loaded  on  a  heater  element  in  a  vacuum  chamber  which  is  purged 
with  nitrogen  three  times.  With  a  small  flow  of  nitrogen  passing  through  the  chamber 
to  keep  it  at  1  atm  pressure,  the  sample  is  heated  to  400®C  for  4  minutes  to  drive  the 
indium  into  the  surface  and  make  on  ohmic  contact  with  the  two-dimensional  electron 
gas  in  the  MODFET  structure.  After  it  is  cool,  the  sample  is  removed  from  the 
chamber  and  the  ohmic  contacts  are  checked  with  a  curve  tracer. 

Gold  wire  is  soldered  to  the  indium  dots  and  the  sample  is  loaded  into  the  Hall 
mobility  measurement  system.  The  system  consists  of  a  Varian  Model  V  4005 
electromagnet  driven  by  a  Varian  Fieldial  Mark  I  field  regulated  magnet  power  supply. 
Electric  current  is  supplied  by  a  Keithley  Model  225  current  source,  and  voltages  are 
measured  by  an  HP  3468A  multimeter.  The  position  and  polarity  of  the  current  source 
is  manually  changed  with  a  switch  box,  and  the  polarity  of  the  magnetic  field  is 
manually  switched  from  the  power  supply. 

The  measurements  are  made  at  77K  by  pouring  liquid  nitrogen  into  the  sample 
chamber.  The  applied  magnetic  field  is  3  kG,  and  the  current  source  ranges  for  1  ^A  to 
1  mA.  The  voltages  are  read  from  the  multimeter  and  entered  on  a  data  sheet.  The 
data  sheet  contains  formulae  which  reduce  the  data  to  sheet  resistance,  sheet  charge 
density,  and  mobility. 

In  this  study,  the  effects  of  annealing  on  the  mobility  was  studied,  so  the  indium 
dots  were  cleaved  off  to  avoid  the  effects  the  indium  would  have  in  the  annealing 
process.  After  annealing,  the  samples  were  6.5  mm  square,  and  the  mobility 
measurement  was  performed  again.  To  compensate  for  differences  in  the  mobility  due 
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to  the  different  sized  samples,  several  samples  were  measured  at  1  cm  and  6.S  mm 
without  annealing,  and  the  mobilities  for  the  annealed  samples  were  normalized  to  the 
pre-annealed  condition. 

All  the  above  experiments  were  performed  at  the  various  facilities  within  Daniels 
Hall  on  the  North  Carolina  State  University  campus.  The  results  of  these  experiments 
are  presented  in  the  following  chapters. 
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4  Heterostructure  Intermixing 


Before  the  effects  of  thermal  disordering  of  quantum  wells  and  superlattices  can 
be  characterized,  an  accurate  method  for  computing  the  electronic  structure  of 
appropriately  shaped  potential  wells  must  be  developed.  This  chapter  describes  the 
computational  models  developed  in  this  work  which  were  used  in  the  studies  discussed 
in  subsequent  chapters. 


4.1  Error  Function  Model  for  Intermixed  Quantum  Wells 


The  most  convenient  model  of  an  intermixed  layer  interface  uses  the  error- 
function-shaped  composition  profile.  This  solution  assumes  no  aluminum  composition 
dependence  in  the  Al-Ga  interdiffusion  coefficient.  The  computer  algorithm  for  this 
model  which  was  originally  used  by  this  group  was  replaced  with  a  faster,  more 
accurate  one  in  this  work.  The  shape  of  the  potential  well  is  a  function  of  the 
aluminum  composition,  x,  which  is  a  superposition  of  the  error  function  profiles  of 
each  interface.  Each  interface  is  treated  as  the  junction  of  two  semi-infinite  layers,  then 
the  profiles  of  each  interface  are  added  to  give 


C{z,t)  =  C„ 


1  +-er/cl 


(4.1) 


where  xis  the  displacement  measured  in  the  direction  perpendicular  to  the  layers,  t  is 
the  annealing  time,  Co  is  the  initial  aluminum  concentration,  I ,  is  the  initial  well  width 
(centered  at  z  =  Q)  and  £>is  the  interdiffusion  coefficient.  The  result  is  depicted  in 
Figure  4.1. 

In  order  to  calculate  the  bound  states,  the  intermixed  potential  is  divided  into 
several  regions,  each  described  by  a  constant  potential  resulting  from  the  average 
aluminum  composition  of  that  region.  At  this  point,  the  well  is  approximated  by  a 
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figure  4.1  Composition  profile  of  an  intermixed  quantum  well.  The  as-grown  well  is 
represented  by  the  heavy  dasiied  line.  After  anne.:Jing,  the  composition  profile  of  the 
bsurriers  sprees  out,  and  is  represented  by  the  light  dashed  curves.  These  profiles  are 
added,  resulting  in  the  composition  profile  given  by  the  heavy  solid  line. 
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stair-step  potential  as  seen  in  Figure  4.2.  This  information  is  then  passed  to  a  main 
program  in  the  form  of  region  widths  and  compositions.  The  bound  state  energy  levels 
are  then  found  by  solving  the  Schrddinger  Wave  Equation  using  the  boundary 
conditions  of  the  wavefimction  and  its  derivative  at  the  region  interfaces. 

For  a  small  number  of  regions  in  the  well,  the  stair-step  approach  results  in  an 
error,  which  is  reduced  as  the  number  of  steps  is  increased.  To  keep  that  number 
reasonable,  an  iterative  approach  is  used  in  which  the  bound  states  are  calculated 
several  times,  doubling  the  number  of  regions  in  the  well  for  each  iteration.  After  a 
few  iterations,  the  values  start  to  converge  to  the  answer  which  represents  an  infinite 
number  of  regions.  An  algorithm  was  employed  which  extrapolates  the  first  few 
iterations  to  that  final  value, 

To  speed  up  the  convergence,  a  method  was  developed  which  concentrates  more 
regions  where  the  composition  varies  quickly  (i.e.,  near  the  original  heterojunction 
interfaces)  rather  than  dividing  the  well  into  evenly  spaced  regions.  This  allows  the 
calculation  to  converge  to  the  bound  state  more  quickly,  as  shown  in  Figure  4.3.  The 
resulting  speed  advantage  is  exploited  when  bound  states  of  many  wells  must  be 
computed.  Because  this  method  is  also  more  accurate  for  a  finite  number  of  regions,  it 
yields  a  closer  approximation  to  the  actual  bound  state  than  a  well  evenly  divided  into 
the  same  number  of  regions.  This  new  method  is  especially  useful  for  occasions  when 
a  single  iteration  is  used. 

The  transition  energies  between  the  conduction  band  electron  state  (ej  and  the 
valence  band  states  for  heavy  holes  (hhj  or  light  holes  OhJ  can  be  computed  using  the 
potential  wells  defined  by  the  bandgaps  of  the  material  in  each  region.  The  bandgap  of 
each  region  of  Al,Ga,.,As  is  given  by 


E 


gap 


1.519+  1.247X- 


5.405X 

T-^204 


(4.2) 
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Figure  4.2  Stair-step  approximation  to  the  intermixed  quantum  well.  The  initial 
approximation  represents  the  intermixed  well  with  8  regions.  Improvements  to  the 
csdculation  of  the  bound-state  energies  are  achieved  as  Uie  number  of  regions  is 
doubled.  The  last  frame  shows  the  continuous  profile  with  the  projected  solution  to  the 
bound-state  energy  for  an  infinite  number  of  regions. 
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Figure  4.3  Convergence  of  the  calculated  bound  state  for  an  intermixed  quantum  well. 
The  dotted  line  represents  the  convergence  for  a  well  which  is  divided  into  uniform 
regions,  while  the  solid  line  shows  the  improvement  of  the  calculation  when  the  regions 
are  concentrated  in  areas  where  the  composition  changes  rapidly  in  z. 
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and  60%  of  the  bandgap  difference  between  adjacent  regions  is  assumed  to  be  in  the 
conduction  band.  This,  along  with  the  appropriate  product  of  the  interdiffiision 
coefficient  and  time,  defines  the  potential  well  shape.  The  bound  states  for  the 
electrons,  heavy  holes,  and  light  holes  are  found  separately  using  the  following  carrier 
effective  masses, 

m*  =  (0.0665  + 0.0835x‘)mo  , 
m;ft  =  (0.45+0.302A:)m,  . 

=  (0.08  + 0.057x)mo  . 

The  gap  between  the  wells  is  added  to  the  electron  and  hole  bound  state  energies  in 
question  to  give  the  total  transition  energy  between  the  conduction  and  valence  band 
states,  which  is  easily  converted  to  a  wavelength. 

Using  this  method,  the  wavelength  shift  due  to  intermixing  versus  the  well  width 
for  various  times  and  temperatures  was  calculated  and  is  presented  in  Figure  4.4.  For 
thick  wells,  the  shift  was  small,  and  as  the  width  decreased,  the  shift  increased.  This 
can  be  explained  as  follows:  as  the  well  intermixes,  the  bottom  of  the  well  begins  to 
fill  from  the  edges.  This  causes  the  lowest  state  to  be  squeezed  upward,  which  occurs 
faster  for  thinner  wells.  But  for  wells  which  were  less  than  a  few  monolayers  thick,  the 
shift  became  small  again,  and  in  fact,  the  shift  from  single  monolayer  wells  was  almost 
negligible! 

The  reason  for  the  reduced  amount  of  shift  is  primarily  because  the  lowest  bound 
state  is  near  the  top  of  the  well,  which  is  actually  getting  wider  as  the  well  becomes 
disordered.  This  allows  the  state  to  relax  to  a  lower  energy  with  respect  to  the  bottom 
of  the  well.  On  the  other  hand,  the  bottom  is  moving  up  in  energy  as  the  well  fills  in. 
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Well  Width  (A) 


Figure  4.4  Transition  wavelength  shifts  for  intermixed  quantum  wells  as  a  function  of 
well  width.  For  large  well  widths,  the  shift  is  small  since  the  bound  states  are  near  the 
bottom  of  the  well,  which  takes  a  long  time  to  fill  as  aluminum  diffuses  in  from  the 
barriers.  For  vei^  thin  wells,  the  bound  state  is  pinned  near  the  tc^  of  the  well,  so 
intermixing  has  little  effect.  Wells  10-50  A  thick,  however,  exhibit  large  shifts  since 
the  bound  state  is  sensitive  to  the  change  in  the  potential  well  as  aluminum  diffuses  in 
from  the  barriers. 
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These  competing  mechanisms  yield  a  small  net  shift  in  wavelength.  This  result  is 
important  is  designing  optoelectronic  devices  which  will  be  subjected  to  high 
temperature  processing. 


4.2  Compositional  Dependence  on  the  Interdiffusion  Coefficient 


While  most  of  the  recent  work  on  intermixed  GaAs/AlGaAs  heterostructures  uses 
an  error-function-shaped  composition  profile,  the  early  data  by  Chang  and  Komai^^i 
clearly  shows  a  dependence  of  the  interdiffusion  coefficient  on  the  aluminum 
composition  in  the  AlGaAs  alloy.  They  used  Auger-electron  profiling  and  the 
Boltzmann-Matano  method  of  determining  the  composition-dependent  interdiffusion 
coefficient,  D  ( x .  7 ),  to  get 


7>(x.7)-Z)„(x)e 


92e'®-^*e 


cm 


(4.4) 


where  x  is  the  mole  fraction  of  aluminum,  k  b  is  Boltzmann's  constant,  7  is  the 
temperature  in  Kelvins,  and  Q  is  the  activation  energy  in  eV. 

To  model  the  composition  profile  resulting  from  a  composition-dependent 
interdiffusion  coefficient,  the  error  function  solution  to  Pick's  Second  Law  is  no  longer 
valid.  This  is  because  the  interdiffusion  coefficient  is  no  longer  independent  of 
displacement,  and  Pick's  Law  must  be  written  more  generally  as 


dC{z,t)  d 


dt 


--I  D(C) 


dCjzA) 

dz 


(4.5) 


Instead,  the  solution  is  found  using  the  first-order  Taylor  series  expansion. 


C(2,t„  + Af)  =  C(2,t  J  +  — 1  D{z) 


dC^zAo) 

dz 


Af 


(4.6) 
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where  t » is  the  initial  time,  A  t  is  the  time  increment,  and  £>is  the  interdiffusion 


coefficient  that  depends  on  composition,  and  hence,  depends  on  displacement  as  well. 

In  this  case,  the  initial  well  profile  is  represented  by  many  evenly-spaced  points, 
each  with  a  composition  C.  Starting  at  °  Oand  using  an  appropriate  interdiffusion 
coefficient,  the  value  of  C  for  each  point  at  to  A  f  is  found.  Resetting  t  o  to  the  current 
time  (to  A d,  the  algorithm  is  repeated.  This  continues  until  the  time  is  incremented  to 
the  total  diffusion  time.  As  an  example,  the  interdiffusion  of  a  40-A  Al^ ^sGa^  55AS 
quantum  well  for  three  hours  at  850°C  was  modeled  using  Equation  (4.4)  for 
interdiffusion  coefficient.  The  results  appear  in  Figure  4.5.  Also  shown  for 
comparison  in  Figure  4,5  is  the  case  for  a  constant  interdiffusion  coefficient-error 
function  profile  using  Equation  (4.4)  withx  =  0.  The  difference  in  the  wavelengths  for 
the  e,-hh,  transitions  was  15  A. 

Since  the  overpressure  conditions  and  material  quality  were  unknown  in  the 
experiments  leading  to  Equation  (4.5),  an  experiment  was  designed  to  determine;  (1)  if 
the  concentration  dependence  is  critical  in  this  work,  and  (2)  the  magnitude  of  that 
dependence.  Since  the  sizes  of  layers  used  here  were  too  small  to  employ  a  method 
such  as  Auger-electron  spectroscopy,  the  Boltzmann-Matano  method  could  not  be 
applied  ,  and  so  another  method  was  needed  to  determine  the  concentration 
dependence. 

If  the  concentration  dependence  is  exponential,  then  the  diffusion  coefficient  is  a 
function  of  the  gallium  and  aluminum  self-diffusion  coefficients,  Z?ca  and  Z)ai, 
respectively,  and  is  given  byf’*! 


Z?(x)  =  Z?„e  ^ 

where 


(4.7) 


P  =  ln 


DoJ 


(4.8) 
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Figure  4.5  Comparison  of  error  function  and  composition-dependent  diffusion  of  a 
GaAs/AlGaAs  quantum  well.  The  dashed  line  represents  the  conduction  band  potential 
well  using  a  constant  interdiffusion  coefficient,  while  the  solid  line  results  from  a 
composition-dependent  interdiffusion  coefficient.  The  initial  well  width  was  40  A  and 
the  barriers  had  35%  aluminum.  Diffusion  at  850°C  for  3  hours  was  simulated,  and 
the  results  showed  a  15-A  difference  in  transition  wavelength  shifts. 
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and  I>  a  is  a  constant  for  a  given  temperature.  While  D  c,  is  well  documented.P^i  very 
little  has  been  published  on  D  ai. 

To  determine  if  a  significant  aluminum  composition  dependence  exists,  an 
annealing  experiment  was  designed  in  which  the  relative  composition  between  the  wells 
and  barriers  was  held  constant,  while  the  absolute  values  were  varied.  If  the 
interdiffusion  coefficients  were  different  for  varying  aluminum  compositions  using  an 
error  function  analysis,  then  further  investigation  would  be  warranted  to  determine  the 
actual  composition  dependence  on  the  interdiffusion  coefficient. 

Three  structures  were  grown  with  a  nominally  10-monolayer  (ML)  (approximately 
28  A)  quantum  well  of  Al,jGa|.,As  with  AlyGa,.yAs  barriers.  The  nominal  composition 
values  of  (x,y)  over  (0,0.2),  (0. 1,0.3),  and  (0.2,0. 4)  for  the  three  structures.  The  basic 
structure  is  shown  in  Figure  4.6.  A  1000-A  well  separated  500  A  from  the  10-ML  well 
was  included  in  each  to  measure  the  actual  well  composition,  since  the  n=l  bound  state 
energy  of  such  a  well  is  negligible  and  does  not  change  appreciably  during  the 
annealing  steps  in  this  study.  After  the  wafers  were  grown,  the  compositions  were 
checked  using  photoluminescence,  and  the  (0,0.2)  wafer  was  close  to  the  nominal 
values,  but  the  (0. 1,0.3)  wafer  was  closer  to  (0.1,0.35).  No  bound  state  or  barrier 
luminescence  was  observed  from  the  (0.2, 0.4)  wafer,  and  it  was  regrown  only  to 
produce  the  same  results.  This  structure  was  therefore  left  out  of  the  experiment, 
which  proceeded  with  only  the  (0,0.2)  and  (0.1,0.35)  structures. 

Several  samples  were  taken  from  each  wafer.  Their  compositions  and  layer 
thicknesses  were  determined  before  annealing.  The  former  was  accomplished  by 
measuring  the  barrier  and  1000-A  well  photoluminescence  (PL)  sprctra.  The  layer 
thicknesses  were  determined  by  calculating  the  bound  state  of  a  quantum  well  with  the 
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Figure  4.6  Structure  to  test  compositional  dq)endence  of  the  interdiffiision  coefficient. 
The  basic  structure  shown  above  contains  a  10-ML  quantum  well  of  which  the 
transition  energies  are  highly  sensitive  to  intermixing.  A  1000-A-thick  well  is  included 
to  monitor  the  initial  composition  of  the  well  material,  and  is  isolated  from  the  thin  well 
by  a  500-A  barrier.  The  well  material  is  Al,jGa,.jAs,  while  the  barriers  are  AlyGa,.yAs 
with  (x,y)  equal  to  (0,0.2)  and  (0.1,0.35)  for  the  interdiffusion  experiment. 
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measured  compositions,  and  then  varying  the  width  of  the  layers  until  the  bound  state 
energy  matched  that  of  the  measured  bound  state  luminescence.  The  samples  were  then 
annealed  at  850®C  in  1  atm  AS4  overpressure  for  0.5,  1.75,  4,  6,  and  7  hours. 

The  new  bound  states  of  the  intermixed  quantum  wells  were  measured  by  photolu¬ 
minescence.  Using  the  parameters  measured  prior  to  the  annealing  step  and  the  error 
function  model  described  above,  the  interdiffusion  coefficient  was  adjusted  until  the 
calculated  bound  state  matched  the  measured  bound  state.  The  resulting  interdiffusion 
coefficients  versus  annealing  time  are  pre^nted  in  Figure  4.7.  The  error  analysis  of 
the  data  is  based  on  a  10- A  error  in  the  photoluminescence  wavelength,  which 
corresponds  to  ±0.003  deviation  in  the  aluminum  composition  for  the  barriers  and  the 
AIq  ,Ga(,9As  well.  For  comparison,  the  same  structures  were  modelled  numerically  and 
a  simulation  of  the  annealing  process  was  performed  using  the  composition-dependent 
interdiffusion  coefficient  given  in  Equation  (4.4).  To  maintain  consistency  of  method, 
the  resulting  shifts  in  bound  state  energies  were  used  in  an  error-function  analysis  and 
plotted  on  the  same  graph. 

The  first  thing  to  note  in  Figure  4.7  is  that  the  error  bars  for  both  sets  of 
experimental  data  overlap,  which  indicates  that  the  differences  due  to  the  aluminum 
composition  are  nearly  immeasurable.  On  the  other  hand,  the  trend  of  the  (0.1,0.35) 
structure  being  consistently  higher  than  the  (0,0.2)  samples  matched  the  predicted 
trends  in  the  calculated  values.  The  upward  trend  in  the  experimental  data  may  also 
indicate  a  composition  dependence,  since  the  average  aluminum  composition  of  the  well 
increases  as  aluminum  diffuses  in.  The  calculated  data  shows  this  trend  also,  but  it  is 
very  slight,  which  means  that  the  magnitude  of  the  upward  slope  is  probably  due  more 
to  temperature  gradients  and  overpressure  fluctuations  in  the  ampoule  than  it  is  due  to 
composition.  Note  that  the  data  for  each  time  were  produced  by  annealing  both 
samples  in  the  same  ampoule,  so  they  were  subjected  to  the  same  conditions. 
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Figure  4.7  Composition-dependent  Al-Ga  interdiffusion  coefficient  as  a  function  of 
annealing  time,  llie  calculated  interdiffusion  coefficient  for  measure  PL  shifts  is 
presented  for  a  10-ML  GaAs/AL  jCag  gAs  well  (triangles)  and  a  10-ML 
Alg  iG^  9As/Alo  3jGao  ^As  well  (squai^).  The  lines  at  the  top  are  for  simulated 
annealing  using  a  comfrasitionally-^ependent  interdiffusion  coefficient  and  using  the 
resulting  energy  shifts  in  an  error-function  analysis  to  find  an  average  interdiffusion 
coefficient.  The  star  corresponds  to  the  Alo  ,G^  9As/AlojjGa(,^As  well  and  the  asterisk 
corre^nds  to  the  GaAs/Alg  ^Gag  gAs  well.  The  cross-hatched  areas  in  the 
experimental  curves  result  from  ah  error  analysis  assuming  a  10-A  error  in  the  PL  data. 
The  overlap  indicates  that  the  compositional  dependence  cannot  be  determined  reliably. 
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Since  the  data  does  not  show  a  strong  compositional  dependence,  further  attempts 
to  accurately  determine  the  composition-dependent  interdiffusion  coefficient  were 
deemed  unnecessary.  The  reasonable  conclusion  of  this  study  then,  is  that  the 
compositional  dependence  of  the  interdiffusion  coefficient  is  negligible  for  annealing 
processes  at  850 °C  in  1  atm  AS4  overpressure.  Based  on  the  data  by  Chang  and  Koma, 
however,  the  dependence  is  greater  at  higher  annealing  temperatures,  and  a  similar 
study  at  higher  temperatures  may  show  a  more  distinct  dependence,  but  such  an 
investigation  is  beyond  the  scope  of  this  thesis. 


4.3  Model  for  Intermixed  Superlattices 


For  calculating  the  dispersion  relationships  of  superlattices,  the  Kolbas-Holonyak 
method  described  in  Section  2. 1  is  adequate  for  a  superlattice  with  a  small  number  of 
periods,  each  divided  into  only  a  few  regions.  An  intermixed  well  requires  many 
regions  to  be  accurately  modeled,  so  many  periods  of  coupled,  intermixed  wells  would 
quickly  become  impractical  with  their  approach.  Therefore,  an  alternate  approach  was 
needed  which  used  a  periodic  potential  based  on  an  error  function  composition  profile. 
The  tight-binding  approach  was  examined  first. 

The  energy  dispersion  relation  for  the  one-dimensional  tight  binding  calculation  in 
the  envelope  function  approximation  was  given  in  Section  2.1,  and  is  repeated  here  for 
convenience.  It  is 

f  =  X  r"Y;(2-ra)[I/(2)-K(2)]Y„C2)d2  . 

r—  (4,9) 

where  is  the  calculated  bound  state  for  a  single,  interdiffused  quantum  well,  r  is  an 


integer,  q  is  the  superlattice  wave  vector,  a  is  the  period  of  the  superlattice,  Y  „  is  the 
wavefunction  of  the  bound  state,  ( 2)  is  the  periodic  potential  defined  by  the 
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superposition  of  all  contributing  error-function-shaped  layer  interfaces,  and  k  ( 2 )  is  the 
potential  of  the  single  interdiffused  quantum  well.  With  nearest-neighbor  interactions 
only  (r  =  0 ,  ±  1),  this  reduces  to 

£  =  £•„- a- 2|3cos((7a)  .  (4-10) 

As  long  as  the  overlap  of  the  wavefiinctions  from  adjacent  wells  is  small,  this 
method  provides  a  good  approximation  to  the  dispersion  relation  calculated  using  the 
Kolbas-Holynak  approach  with  several  coupled  wells.  For  larger  wavefunctions 
resulting  from  smaller  dimensions,  the  solutions  depart.  Adding  next  nearest  neighbor 
terms  actually  causes  the  solutions  to  differ  even  more.  This  suggests  that  the  values  of 
the  calculated  wavefunctions  far  from  the  immediate  vicinity  of  the  single  well  are  not 
accurate  enough  to  give  consistent  solutions  for  short  periods. 

A  Kronig-Penney-type  calculation  was  examined  next.  For  this  study,  a  method 
used  by  Yuh  and  Wang,^**!  and  introduced  in  Section  2.1,  was  adapted  to  handle 
interdiffused  superlattices.  This  approach  defines  a  2  x  2  transfer  matrix  for  each  layer 
in  the  basis,  given  by 


C0S(/Cji() 

—  sinikiLi) 
rrii 


nii 

cos{kiLi) 


where  i  denotes  the  layer,  /.j  is  the  layer  thickness,  and 


(4.11) 


ki  = - r - 

ft  (4.12) 

The  transfer  matrices  of  each  region  are  multiplied  together  to  give  the  total 
matrix,  S.  Using f  ( 2 ),  a  vector  with  elements ri(2)  =  nj(2)and 
r2(^)“V'(2)/^*,  the  eigenvalues  of  S  are  found  using 

r(r-^a)  =  Sr(2)  =  Mr(2)  .  (4.13) 
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Along  with  Bloch  theory,  which  gives  X.  =  e**"",  we  find 

S„  +  S22  =  G‘''“-^e-“'“  =  2cos(ga)  .  (4-14) 

where  q  is  the  superlattice  wavevector.  If  a  simple,  two-layer  basis  is  used,  the  result 

is  the  familiar  Kronig-Penney  solution  in  the  effective  mass  approximation. 

This  model  was  used  to  calculate  the  energy  dispersion  of  a  GaAs-AlAs 
superlattice  with  35-A  wells  and  15-A  barriers.  The  results  are  shown  in  Figure  4.8 
representing  several  annealing  times  up  to  4  hours  using  an  interdiffusion  coefficient, 

£)  =  8  X  10'^’  cm^s  *  (850®C,  1  atm  AsJ.  This  structure  was  grown  and  annealed  as 
well,  and  the  PL  spectra  are  presented  in  Figure  4.9.  Notice  that  the  spectra  all  contain 
several  distinct  peaks,  and  the  second  peak  is  dominant  in  the  as-grown  sample,  as  well 
as  in  the  sample  annealed  for  one  hour. 

The  calculated  transition  wavelengths  from  the  critical  points  of  the  n= 1  electron 
subband  are  also  shown  in  Figure  4.9,  which  indicate  that  the  second  large  peak  is 
probably  due  to  the  ej-lh,  transition.  The  relative  magnitude  of  this  peak  can  be 
explained  by  examining  the  band  diagrams  of  Figure  4.8.  For  the  as-grown  and  the 
1-hour  diffused  case,  the  n=2  heavy  hole  lines  up  with  the  n= 1  light  hole  at  the 
mini-zone  edge.  As  holes  thermalized  in  the  n=2  heavy  hole  band,  some  will 
recombine  at  the  band  edge  (e'j-hhj  transition,  with  the  prime  indicating  the  band-edge 
saddle  point),  while  most  will  scatter  into  the  n=l  light  hole  band.  From  there,  they 
will  either  recombine  at  the  saddle  point  (e',-lh',),  or  thermalize  further  to  recombine  at 
q  =  0(e,-lh,). 

Time-resolved  photoluminescence  (TPL)  measurements  support  this  explanation. 
In  all  cases  above,  the  measured  lifetime  of  the  second  peak  was  longer  than  that  of  the 
first  peak,  which  is  logical  if  holes  are  being  supplied  to  the  n=l  light  hole  subband 
from  the  n=2  heavy  hole  subband  over  a  finite  time.  This  is  supported  further  by  the 


-58- 


WAVEVECTOR  (MINI-ZONES) 


WAVEVECTOR  (MINI-ZONES) 


WAVEVECTOR  (MINI-ZONES) 


WAVEVECTOR  (MINI-ZONES) 


Figure  4.8  Energy  dispersion  curves  for  an  intermixed  GaAs/ AlAs  superlattice,  ^e 
frames  above  show  the  energy  diagrams  of  a  GaAs/ AlAs  si^rlattice  grown  vdA  35-A 
wells  and  15-A  barriers  and  annealed  for  30  60,  120,  and  240  minutes  ((a),  (b),  (c), 
and  (d)  respectively)  for  the  first  two  mini-zones.  The  heavy  solid  lines  represent  toe 
electrons  and  heavy  holes,  while  toe  dashed  lines  represent  the  light  holes.  TJe  lighter 
lines  are  for  toe  as-grown  samples.  The  energy  spans  the  gap  of  AlAs,  and  the  dotted 
lines  between  toe  AlAs  and  GaAs  band  edges  are  the  positions  of  the  top  of  the 
intermixed  barriers. 
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Figure  4.9  Photoluminescence  spectra  of  an  annealed  GaAs/AlAs  superlattice.  The 
PL  spectra  of  a  GaAs/AlAs  superlattice  grown  with  35-A  wells  and  15-A  barriers  is 
shown  as-grown  (a)  and  after  annealing  for  30  minutes  (b)  and  60  minutes  (c)  at  8S0°C 
in  1  atm  As^.  The  calculated  transition  wavelengths  for  critical  points  of  the  first 
electron  subband  and  the  first  two  hole  subbands.  The  saddle  points  of  the  minibands 
are  denoted  by  a  prime. 
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measured  decay  of  the  luminescence  intensity  at  the  second-major-peak  wavelength. 
The  log  of  the  decay,  shown  in  Figure  4.10,  indicates  that  the  band  was  being  with 
supplied  with  carriers  during  the  first  few  hundred  picoseconds,  where  the  lifetime  is 
very  long  (the  slope  is  very  small  at  first).  After  about  700  ps,  the  curve  reaches  a 
constant  slope  which  yields  what  is  presumably  the  e,-lhj  lifetime  after  the  n=2  heavy 
hole  band  is  empty. 

While  both  the  PL  and  TPL  data  support  the  superlattice  model,  there  are  several 
bumps  in  the  PL  spectra  which  do  not  correspond  to  predicted  transitions  between 
critical  points  of  the  subbands.  These  features,  as  well  as  the  divergence  of  the 
calculated  and  measured  critical  points  for  longer  anneal  times  (PL  data  not  shown) 
may  be  explained  by  valence  band  mixing.r^J  If  the  zones  of  Figure  4.8  are  folded,  the 
hhj  band  and  the  Ihj  band  touch  or  cross.  At  these  points,  anticrossings  will  occur 
which  will  form  new  critical  points,  and  the  Ih,  band  will  take  on  some  characteristics 
of  the  hhj  band,  and  vice-versa.  These  effects  are  not  included  in  this  model. 
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Figure  4.10  Luminescence  decay  of  the  electron-to-light-hole  transition.  The  decay  of 
the  second  major  luminescence  p^k  of  a  GaAs/AlAs  superlattice  with  35-A  wells  and 
IS-A  barriers  is  shown  (a).  The  lifetime  is  estimated  using  a  least-squares  algorithm  to 
find  the  slope  of  the  decay.  Notice  that  the  data  does  not  become  linear  until  about 
700  ps,  in  contrast  to  a  more  typical  decay  shown  in  (b).  The  data  in  (a)  begins  with  a 
long  characteristic  lifetime,  probably  because  it  represents  the  e,-lh,  transition  with  the 
Ih,  band  being  filled  by  the  hh^  band  at  the  mini-zone  edge. 
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5  Selective  Intermixing 


In  order  to  exploit  the  various  attributes  of  intermixed  quantum  wells  and 
superlattices,  the  methods  of  producing  such  structures  must  be  controlled  so  that 
selected  areas  of  a  sample  become  intermixed,  while  other  areas  remain  preserved.  In 
this  chapter,  three  methods  of  selectively  intermixing  semiconductor  quantum  heteros¬ 
tructures  are  examined.  They  are:  (1)  using  a  focused  ion  beam  to  create  localized 
areas  of  excess  gallium  near  a  superlattice,  (2)  encapsulating  with  the  dielectric  SiN„ 
and  (3)  fabricating  features  in  GaAs  grown  at  low  substrate  temperature  to  provide 
excess  arsenic  to  local  areas  on  a  chip. 

5.1  Localized  Gallium  Overpressure 

The  use  of  a  gallium  overpressure  was  shown  by  Hsieh  et  al.  to  suppress  the 
interdiffusion  of  quantum  wells  exposed  to  high  temperatures. On  the  other  hand, 
Hirayama  et  a/. Poland  Cibert  et  showed  that  gallium  implanted  into  the  heteros¬ 
tructure  enhanced  the  interdiffusion  process.  In  this  thesis,  gallium  atoms  were 
implanted  near  the  surface,  away  from  the  heterostructure,  in  order  to  generate  local 
gallium  overpressures  which  should  suppress  the  Al-Ga  interdiffusion  at  the  layer 
interfaces. 

5.1.1  Gallium  Focused  Ion  Beam 

A  sample  was  grown  with  a  55-period  GaAs/Alo  jiGa^  gjAs  superlattice  with  35- A 
wells  and  barriers,  surrounded  by  Alg  jjGa^,  ^jAs.  The  superlattice  was  5000  A  from  the 
surface  of  the  sample.  A  gallium  focused  ion  beam  (FIB)  was  used  to  write  100  x  100  % 
m  squares  with  implant  energies  of  15,  20,  and  25  keV  in  doses  between  2  x  i  o cm-^ 
to  4  X  10'^  cm These  energies  correspond  to  implant  depths  on  the  order  of  100  A, 
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so  the  ions  were  far  enough  above  the  superlattice  not  to  cause  any  damage.  Most  of 
the  doses  were  below  10*^  crn-^,  which  is  the  threshold  at  which  implant  damage 
becomes  the  dominant  mechanism  for  generating  defects  as  reported  by  Hirayama  et 

The  photoluminescence  (PL)  spectra  of  the  patterns  were  taken  as  a  baseline.  The 
spectra  of  the  implanted  features  showed  a  very  broad  peak  for  energies  in  the  bandgap 
which  corresponded  to  luminescence  from  the  implant-damaged  surface.  This  peak  was 
absent  in  the  unimplanted  regions  between  features.  The  15-keV-implanted  features 
were  annealed  at  850®C  with  no  overpressure  for  3  hours,  and  the  20-  and  25-keV- 
implanted  features  were  annealed  at  825®C  under  the  same  overpressure  and  time 
conditions. 

The  PL  spectra  were  taken  again  after  the  annealing  step  and  compared  to  the 
pre-annealed  spectra.  The  broad  peaks  from  the  initial  implant  damage  were  gone. 

The  positions  of  the  n=l  subband  luminescence  peaks  shifted  in  wavelength  less  than 
100  A,  which  corresponds  to  an  interdiffusion  coefficient  of  2  x  1 0" ernes  ',  and  is 
commensurate  with  the  peak  shifts  expected  for  annealing  in  a  gallium  overpressure. 

Finally,  a  scan  of  the  PL  peaks  across  and  between  the  15  keV,  2.3  x  i  o  ‘°cm-2 

features  was  taken  to  compare  the  peak  shifts  of  the  implanted  and  unimplanted  areas. 
Surprisingly,  was  no  measurable  difference  was  noted.  This  means  one  of  two  things. 
First,  the  gallium  implanted  into  the  surface  had  no  effect,  and  the  results  are  due  to  the 
annealing  process  with  no  overpressure.  Second,  the  defects  generated  by  the  implant 
were  predominantly  highly  mobile  gallium  interstitials  (Iq.)  that  diffused  outward  from 
the  features  as  much  as  100  %m  or  more,  thereby  preserving  the  superlattice  around  the 
features  as  well.  Since  the  peak  shifts  agree  wiih  previous  data  resulting  from  the 
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gallium  overpressure  annealing  process, the  second  explanation  is  worthy  of  further 
investigation  when  the  opportunity  to  obtain  FIB-written  samples  again  becomes 
available. 

5.1.2  Electron  Mobility  Experiment 

The  possibility  of  using  FIB  technology  to  write  small  areas  of  gallium  to  preserve 
large  areas  of  a  quantum  structure  suggested  some  ideas  which  led  to  a  side  experiment 
involving  a  high-electron  mobility  transistor  (HEMT)  structure  (or  modulation-doped 
field  effect  transistor  (MODFET)  structure).  The  HEMT  device  relies  on  the  sharp 
interface  between  an  n-type  AlGaAs  layer  and  an  intrinsic  GaAs  layer.  The  conduction 
band  bends  as  the  Fermi  levels  of  the  two  materials  line  up,  which  causes  a  notch  to 
form  at  the  junction.  This  notch  contains  quantum  bound  states  into  which  electrons 
from  the  donors  on  the  AlGaAs  side  fall.  The  electrons  are  then  in  a  two-dimensional 
(2D)  electron  gas,  which  is  separated  from  the  ionized  impurity  scattering  sites, 
resulting  in  high  mobilities.  If  the  interface  is  not  sharp,  then  the  quantum  well  smears 
out  and  the  2D  electron  gas  is  lost. 

To  determine  the  effects  of  a  gallium  overpressure  on  an  annealed  HEMT 
structure,  an  experiment  was  designed  in  which  the  mobility  of  a  HEMT  structure  was 
measured  before  and  after  annealing  with  no-overpressure,  gallium-overpressure,  and, 
for  completeness,  arsenic-overpressure.  An  unprocessed,  three-inch  production  HEMT 
wafer  was  provided  by  the  Martin  Marietta  Corporation.  The  epitaxial  layers  consisted 
of  a  l-?4m  GaAs  buffer  grown  on  the  substrate,  followed  by  a  20- A  AIq  26030.74^5 
undoped  spacer,  a  400- A  Alg  26030.74^8  Si-doped  (/Vo  =  3x  10‘®)  layer,  and  then  a 
400-A  Si-doped  (Ao  =  4x  10*®)  GaAs  contact  layer.  This  is  a  structure  in  which  the 
channel  is  normally  on. 
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The  wafer  was  cleaved  into  1-cm  squares,  and  the  mobility  was  measured  to 
assess  its  uniformity.  The  mobilities  were  between  about  21,000  and  28,000  cmW  ^s  * 
at  77K.  Because  the  indium  used  for  the  contacts  would  contaminate  the  atmosphere 
inside  the  ampoules  during  annealing,  either  separate  samples  had  to  be  used  for  the 
annealed  structures  or  the  indium  contacts  had  to  be  cut  off.  Because  of  the  wide  range 
of  mobilities  in  the  as-grown  samples,  the  latter  method  was  chosen.  First,  however, 
several  of  the  previously  measured  samples  were  cut  down  to  6.5-mm  squares,  and 
re-measured  to  determine  if  the  size  difference  affected  the  measured  mobility.  The 
measured  mobility  of  the  smaller  samples  were  on  average  higher  by  about 
1000  cm^V-^s  *,  and  this  difference  was  taken  into  account  in  subsequent  tests. 

The  annealing  steps  were  performed  at  850°C  for  20  and  60  minutes.  The 
mobilities  of  the  arsenic-overpressure  samples  could  not  be  measured,  which  is 
probably  a  result  of  the  quantum  well  being  completely  intermixed  as  expected.  The 
measured  mobilities  of  the  samples  without  an  overpressure  for  annealing  times  of  20 
and  60  minutes  were  10.2%  and  8.0%  of  their  original  (adjusted)  values,  respectively. 
Finally,  for  the  gallium  overpressure,  the  mobility  after  20  minutes  was  10.0%  of  the 
original  value,  while  it  could  not  be  measured  after  60  minutes.  The  poorer  mobilities 
for  the  gallium  overpressure  samples  was  puzzling  at  first,  but  since  the  HEMT 
structure  is  a  surface  device,  the  mobility  was  dominated  by  the  effects  of  the  annealing 
process  on  the  surface.  In  the  gallium  atmosphere  at  high  temperature,  arsenic  atoms 
come  out  of  the  surface  faster  than  with  no  overpressure.  (The  surface  morphology  is 
commonly  preserved  by  keeping  arsenic  atoms  in  with  an  arsenic  overpressure.)  The 
defects  in  the  arsenic-poor  surface  will  then  drive  the  surface  depletion  region  deeper 
into  the  structure,  which  eventually  depletes  the  channel. 
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To  assess  the  capability  of  a  small  gallium-implanted  area  to  preserve  a  large  area 
surrounding  it,  a  more  appropriate  experiment  might  be  to  encapsulate  an  implanted 
structure  in  Si3N4  to  maintain  the  surface  morphology,  and  then  proceed  with  the 
annealing  process.  An  alternate  approach  to  encapsulation  would  be  to  use  rapid 
thermal  annealing  with  the  surface  of  the  HEMT  sample  in  contact  with  a  bare  GaAs 
wafer. 

5.2  SiN^  Encapsulation 

Another  method  which  was  investigated  for  use  in  selective  intermixing  was  to 
pattern  stripes  of  the  dielectric  SiN,  on  the  sample  surface.  This  technique  preserves 
the  stoichiometry  under  the  dielectric,  which  reduces  the  generation  of  point  defects 
during  high  temperature  processes.  The  exposed  surfaces  will  react  with  the  ambient 
overpressure  which  can  generate  enough  defects  to  enhance  the  interdiffusion  process. 

To  verify  this,  a  sample  was  prepared  which  had  a  21-period  GaAs/AlAs 
superlattice  with  79-A  wells  and  40-A  barriers.  The  superlattice  was  bounded  on  either 
side  by  Alg  jGaj^As  (l-^m  buffer  and  0.75-?4m  cap).  A  lOOO-A  SiN,  film  was  deposited 
on  the  sample,  and  lOO-^m-thick  stripes  spaced  100  %m  apart  were  fabricated  in  the  film. 
PL  spectra  were  taken  for  both  the  open  and  encapsulated  surface,  and  were  the  same 
for  both  types  of  areas. 

The  sample  was  then  annealed  at  850°C  for  3  hours  with  no  overpressure.  The 
PL  spectra  were  taken  again,  and  they  are  presented  in  Figure  5.1.  The  encapsulated 
stripes  experienced  a  wavelength  shift  in  the  luminescence  peak  of  only  40  A,  while  the 
exposed  stripes  shifted  200  A,  or  five  times  that  of  the  preserved  stripes.  The 
corresponding  interdiffusion  coefficients  were  1 .6  x  1 0*  cm^s  *  and 
6.7  X  10'*’  ernes  ',  respectively,  which  resulted  in  a  4.2:1  selectivity  ratio. 
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Figure  5.1  Photoluminescence  spectra  of  a  selectively  intermixed  GaAs/AlAs 
superlattice  via  SiN,  encapsulation.  The  spectra  are  for  a  GaAs/AlAs  superlamce^ 
grown  with  79-A  wells  and  40-A  barriers,  both  before  and  after  ann^ing  at  850  C  tor 
3  hours.  The  major  peak  is  due  to  the  lowest  state  of  the  first  well  of  the  ^perlatuce, 
which  was  asymmetric  since  it  was  bounded  by  AIq  sGa^ jAs  on  one  side.  The  lowest 

superlattice  state  transition  wavelength  shift  only  40  A  when  covered  with  a  SiN,  stripe, 
and  200  A  when  uncovered.  This  results  in  a  selectivity  ratio  of  the  interdiffusion 
coefficient  of  4.2:1. 
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These  results  indicate  that  this  method  of  selective  intermixing  is  a  viable 
technique  to  selectively  alter  the  electronic  properties  of  a  superlattice.  Notice  that 
since  both  regions  experienced  shifts  in  the  n=l  subband  energy  levels,  the  optical 
properties  would  be  changed  in  both  also.  This  particular  structure  was  designed  to 
result  in  the  largest  difference  in  the  weighted  average  index  of  refraction  of  the  AlAs 
and  GaAs,  and  the  index  of  the  corresponding  alloy,  AlojjGao  J7AS.  Since  the  wells 
were  neither  completely  preserved  nor  intermixed,  the  indices  of  both  the  "preserved" 
and  the  "intermixed"  regions  must  fall  somewhere  between  the  extremes.  The  problem 
of  determining  the  index  of  partially  intermixed  wells  is  examined  in  Chapter  6. 

Finally,  the  extent  of  intermixing  in  the  lateral  direction  at  the  boundaries  between 
the  open  and  covered  areas  needs  to  be  addressed  before  devices  such  as  optical 
waveguides  with  widths  of  only  S  to  10  %m  can  be  fabricated  using  this  technique. 
Methods  of  optimizing  this  technique  should  also  be  considered,  such  as  annealing  in  an 
arsenic-overpressure,  or  implanting  the  exposed  surface  with  arsenic  ions  to  enhance 
the  intermixing  process. 

5.3  Low-Temperature  Gallium  Arsenide 

Low-temperature  (LT)  GaAs  refers  to  a  GaAs  epitaxial  layer  grown  by  molecular 
beam  epitaxy  (MBE)  at  a  low  substrate  temperature  of  200 °C  (as  opposed  to  640 °C). 
This  growth  process  results  in  a  layer  which  is  arsenic-rich, poj  and  has  been  shown  to 
enhance  the  Al-Ga  interdiffusion  coefficient  when  used  as  a  buffer  layer.  In  this 
work,  an  LT-GaAs  buffer  layer  was  patterned  with  stripes  to  provide  local  arsenic 
overpressures. 

First,  a  0.2- %m  LT-GaAs  layer  was  grown  on  a  GaAs  substrate.  The  layer  was 
processed  using  photolithography  and  a  wet  chemical  etch  to  produce  lOO-^m-wide 
LT-GaAs  stripes  separated  by  1(X)  ?4m  on  top  of  the  original  wafer.  The  wafer  was  then 
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degreased  and  an  abbreviated  surface-preparation  etch  was  used  to  avoid  etching  away 
the  LT-GaAs  stripes.  The  wafer  was  loaded  back  into  the  MBE  for  a  second  growth. 
A  structure  was  grown  on  top  of  the  processed  wafer  which  consisted  of  a  0.5-%  m 
Alo.3sGao.65As  layer,  an  87-A  GaAs  quantum  well  layer,  and  another  0.5- %m 
Alfl  ssGao  osAs  layer.  Although  most  of  the  wafer  was  hazy,  indicating  poor  crystalline 
structure,  sufficient  parts  of  it  were  clear  enough  to  continue  an  annealing  experiment. 

The  quality  of  the  quantum  wells  was  verified  by  taking  the  PL  spectra  of  several 
stripes,  which  also  established  a  baseline  for  the  experiment.  The  sample  was  then 
annealed  in  an  ampoule  with  no  overpressure  for  4  hours  at  850°C.  The  PL  peaks 
measured  after  the  annealing  step  showed  a  wavelength  shift  of  about  210  A  from  the 
wells  grown  on  the  LT-GaAs  buffer  layer  stripes,  and  about  190  A  from  the  other 
wells.  The  data  is  presented  in  Figure  5.2. 

Although  the  differences  are  small,  they  are  consistent  between  adjacent  wells, 
which  suggests  that  the  Al-Ga  interdiffusion  coefficient  in  the  structures  grown  on  the 
LT-GaAs  stripes  was  enhanced  over  that  of  the  adjacent  stripes.  The  small  difference 
in  peak  shifts  could  be  a  result  of  the  LT-GaAs  stripes  being  thinned  too  much  during 
the  surface  etch  in  the  sample  preparation  step  before  the  second  growth,  as  the  etch 
rate  was  not  calibrated.  An  attempt  to  grow  additional  layers  on  a  sample  with 
LT-GaAs  stripes  without  using  a  surface  preparation  etch  resulted  in  a  completely 
unacceptable  surface.  Clearly,  this  technique  needs  further  refinement  in  order  to 
produce  samples  with  enough  difference  in  the  interdiffusion  rates  to  be  a  useful 
method  of  selective  intermixing. 
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Figure  5.2  Photoluminescence  peaks  of  a  selectively  intermixed  quantum  well  using 
low-temperature  GaAs.  The  structure  was  an  87-A  GaAs/ALjjGaojjAs  quantum  well 
grown  on  stripes  of  LT-GaAs  and  annealed  for  4  hours  at  850*C.  The  wells  over  the 
LT-GaAs  shifted  to  a  lower  wavelength  than  the  wells  grown  over  a  conventional  GaAs 
buffer.  The  overall  trend  in  the  wavelength  across  the  wafer  is  probably  due  to  non¬ 
uniformities  in  the  growth. 
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6  Optical  Properties  of  Intermixed  Quantum  Wells 


All  of  the  optical  properties  of  an  intermixed  quantum  well,  as  with  any 
semiconductor,  are  contained  in  the  complex  dielectric  constant 

e(co)  =  e'(tA))+ i6"(a))  ,  (6.1) 

and,  in  general,  are  a  function  of  frequency  (see  Appendix  A). 

Absorption  and  reflection  properties  can  be  found  by  relating  £  ( oo )  to  the 
complex  index  of  refraction  by 

€(0U)  o  , 

(6.2) 

where  e  o  is  the  free-space  permitivity ,  n  (  oo  )  is  the  real  part  of  the  index  and  ic  ( u) )  is 

the  extinction  coefficient.  The  absorption  coefficient  a  (oo)  is  related  to  the  extinction 
coefficient  by 

"  (6.3) 

where  c  is  the  speed  of  light.  Using  G"/£o  =  2nic,  a((A})is  related  to  the  imaginary 

part  of  e(a))  by 


a(cju)  = - e"(tx>)  . 

G„cn(oo) 


Finally,  the  reflectivity  at  normal  incidence  is  given  by 


/?(uj)  = 


(n(U))-  1)"-^(K(U3))  = 
(n((Aj)+  l)2  +  (K(a)))  = 


In  order  to  determine  the  actual  values  of  the  optical  properties  of  an  intermixed 
quantum  well,  therefore,  g  ( u) )  must  be  calculated.  In  the  oscillator  approach  to 
determining  6(00 ),  tlie  complex  dielectric  consta  it  is 
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e(a))-e„-HA/X 


/y 


y  jix>  (6  5) 

The  dielectric  function  is,  then,  the  sum  of  the  dielectric  responses  of  all  the  oscillators 
in  the  system,  each  with  a  resonant  frequency  ,  oscillator  strength  /  y ,  and  damping 
coefficient  y  y.  The  oscillator  strength  is 

^  — |M|^ 

(6.7) 


fr 


mgfuxi  y 

where  |  M  I  ^  is  the  momentum  matrix  element  of  the  oscillator. 


The  frequency  is  calculated  by  solving  for  the  energy  levels  in  the  quantum  well 
as  described  earlier.  This  also  yields  the  envelope  wavefunctions  of  the  electron  and 
hole  states,  which  is  necessary  for  calculating  the  matrix  element. 

In  order  to  isolate  the  dielectric  constant  of  the  intermixed  quantum  well,  the  sum 
of  all  the  oscillators  in  Equation  (6.6)  is  restricted  to  the  energy  levels,  or  oscillators  in 
the  quantum  well.  The  remaining  oscillators  of  the  system  are  combined  with  to 
form  a  background  dielectric  constant,  This  assumes  that  the  changes  in  the 
dielectric  constant  are  due  solely  to  the  modifications  of  the  bound  states  in  the  well, 
since  the  electronic  states  of  the  bulk  barrier  material  remain  unaffected  by  the 
annealing  process. 


6.1  The  Matrix  Element  of  Intermixed  Quantum  Wells 

Except  for  the  bound  state  energies,  the  matrix  element  is  the  most  important 
factor  in  determining  the  optical  properties  of  a  quantum  well.  It  describes  the 
interaction  between  states  in  the  conduction  and  valence  bands,  and  directly  influences 
the  oscillator  strength  and  the  radiative  lifetime  of  excited  states.  Following  is  the 
theoretical  development  of  the  matrix  element  as  it  relates  to  intermixed  quantum  wells, 
as  well  as  the  description  of  an  experiment  to  observe  the  calculated  changes  to  it. 
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6.1.1  Theory 


The  matrix  element  for  a  single  oscillator  between  two  states,  m  and  n,  with  a 
single  polarization  is 

where  <  n  |  m  >is  the  overlap  of  the  electron  wavefunction  in  state  n  and  the  hole 

wavefiinction  in  state  m,  Pmn  is  the  interband  matrix  element  including  the  selection 
rules,  and  the  delta  function  conserves  momentum  in  the  directions  parallel  to  the  layers 
in  the  quantum  well,  described  by  the  wavevectors  it  „ , ,  and  A:  m  1 1  (see  Appendix  B).  The 
interband  matrix  element  in  III-V  semiconductors  for  electron-to-heavy-hole  (TE) 
transitions  at  Ac  ^  -  0  is 


\P 


mn 


(6.9) 


where  the  factor  of  ;  comes  from  ^  for  the  average  over  all  directions  ^  for  bulk)  and  I 


for  the  conservation  of  spin.  P„is  found  using  the  Kane  model!*’! 

m. 


2  0  ^  y  I 

‘  0 


1£_  1  I 


m. 


(6.10) 


where  mo  is  the  free-electron  mass,  m^  is  the  conduction-band  effective  mass,  is  the 


bulk  material  bandgap,  and  A  is  the  split-off  energy  for  the  split-off  valence  band. 

For  GaAs/AlGaAs  quantum  wells  with  thicknesses  greater  than  about  40  A,  the 
matrix  element  can  be  considered  a  constant,  taking  the  value  of  GaAs,  which  is  about 
|F„,„|^=  1 .87  mo  Fg  for  a  quantum  well  (1 .2  SmpFj,  for  the  bulk  material).  For 
thinner  wells,  a  significant  portion  of  the  envelope  wavefunctions  extends  into  the 
AlGaAs  barriers  where  the  interband  matrix  element  is  reduced.  For  instance,  the 
value  of  Fo  is  1.15  times  greater  in  GaAs  than  it  is  in  Alo  jGao^As.  To  include  this 
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effect  in  the  calculation  of  the  matrix  element  for  thin  wells,  the  interband  matrix 
element  for  each  layer  must  be  calculated  and  then  weighted  by  the  amount  of  overlap 
of  the  electron  and  hole  wavefunctions  for  each  region.  Equation  (6.8)  must  then  be 
modified  to  give 

-ti 

~  T 

2 

where  Y„,  and  Y„  are  the  envelope  wavefunctions  of  the  carriers,  and  momentum 

conservation  is  implied.  A  schematic  for  calculating  M  „„  is  given  in  Figure  6.3.  The 
barriers  are  denoted  as  regions  I  and  m  while  the  well  is  region  U.  The  matrix  element 
normalized  to  Fo  for  GaAs  for  several  cases  is  shown  in  Figure  6.4.  The  solid  line 
depicts  the  theory  for  a  GaAs  well  between  Alo3Gao7As  barriers.  For  thin  wells,  the 
overlap  of  the  electron  and  hole  wavefunctions  is  reduced  because  of  the  differences  in 
the  effective  masses  of  the  carriers,  which  causes  the  curve  to  diverge  from  the  ideal 
2D  case  of  an  envelope  overlap  of  unity.  Including  the  effect  of  the  barrier  interband 
matrix  element  results  in  a  further  reduction  of  the  matrix  element  for  ultrathin  wells 
compared  to  the  use  of  only  the  GaAs  interband  matrix  element.  The  matrix  elements 
for  both  the  electron-to-heavy-hole  (e-hh)  and  electron-to-light-hole  (e-lh)  transitions 
with  the  barrier  effects  as  well  as  the  GaAs-only  approximation  are  shown  in 
Figure  6.5. 

In  recent  work  performed  by  Benjamin,**®^  this  model  of  the  optical  matrix 
element  was  used  in  the  calculation  of  the  dynamics  of  stimulated  photon  emission  from 
the  n=l  e-hh  transitions  in  54-A  and  17-A  single  quantum  well  heterostructure  lasers. 
When  the  matrix  element  was  used  to  calculate  the  gain  coefficient,  a  reduced  value 
was  obtained  compared  to  a  model  which  is  valid  for  thick  wells.  Although  both 
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Figure  6.1  Schematic  for  calculating  the  matrix  element  of  a  thin  quantum  well.  The 
electron  and  hole  envelope  wavefunctions  are  shown  with  different  shading  for  each 
region  of  the  structure  to  emphasize  that  each  region  is  considered  separately.  The 
overlap  of  the  wavefunctions  are  calculated  for  each  region  and  used  as  a  weight  to  the 
interband  matrix  element  of  the  corresponding  region.  The  total  matrix  element  is  then 
the  sum  of  the  matrix  elements  for  each  region. 


Fisure  6.2  Matrix  element  for  e,-hh,  transitions  in  a  GaAs/AlojGao  yAs  quantum 
well  as  a  function  of  well  thickness.  The  ej-hh,  matrix  element  (heavy  solid  line) 
is  less  than  the  ideal  2D  case  for  thin  wells  because  of  the  reduced  overlap  of  the 
electron  and  hole  wavefunctions.  The  consideration  of  the  lower  bamer 
interband  matrix  element  causes  a  further  reduction  in  the  overall  matnx 
element  for  ultrathin  wells  since  the  wavefunctions  are  almost  enurely  in  the 

barriers. 
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Figure  6.3  Matrix  elements  as  a  function  of  thickness  for  ej-hh,  and  e,-lh, 
transitions  in  a  GaAs/Alg  jGag  7AS  quantum  well.  Notice  that  the  matrix 
elements  for  thin  wells  get  smjdler  due  to  the  overlap  between  electron  and  hole 
wavefiinctions.  The  dotted  lines  are  for  the  case  where  the  difference  in  the 
interband  matrix  elements  between  the  well  and  barrier  materials  are  neglected. 
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models  should  agree  for  thick  wells,  these  thin  wells  have  carrier  wavefunctions  which 
extend  into  the  barriers,  and  the  overlap  between  electron  and  hole  wavefunctions  is 
reduced,  which  is  not  taken  into  account  in  the  thick-well  model.  For  the  54-A  well, 
the  thick-well  gain  coefficient  calculated  using  the  thick  well  model  was  8.52  times  that 
of  the  model  using  this  matrix  element.  (Benjamin  reports  a  factor  of  4.26,  but  he 
neglected  spin  conservation  in  his  calculation.)  To  fit  the  experimentally  measured 
lasing  delay  to  the  calculated  delay,  the  gain  coefficient  had  to  be  scaled  down  by  an 
additional  factor  of  1.15.  For  the  17-A  well,  the  coefficient  was  scaled  down  by  1.60 
to  fit  the  data.  The  agreement  would  be  closer,  however,  if  the  electron-hole 
wavefunction  overlap  was  used  instead  of  only  the  electron  wavefunction  in  calculating 
the  overlap  between  the  carrier  wavefunction  and  the  optical  field. 

The  theory  is  easily  extended  to  intermixed  wells  by  dividing  the  well  into  N 
regions  such  that 


M 


mn 


imO  •'I 


(6.12) 


The  results  for  several  well  thicknesses  as  a  function  of  annealing  time  are  shown  in 
Figure  6.6.  Notice  that  the  matrix  element  gets  larger  in  time  for  the  e-hh  transitions  in 
the  thinnest  wells,  is  relatively  constant  for  e-lh  transitions  in  the  thinnest  wells  and 
e-hh  transitions  in  the  thicker  wells,  and  decreases  for  e-lh  transitions  in  the  thicker 
wells.  This  can  be  explained  by  considering  the  aluminum  composition  dependence  of 
the  interband  matrix  element,  and  the  overlap  of  the  electron  and  hole  envelope 
wavefunctions,  which  are  competing  mechanisms.  Each  case  is  explained  separately 
below.  Note  that,  in  general,  while  an  8-monolayer  (ML)  well  is  considered  thin,  in 
this  context,  it  is  considered  "thick"  when  compared  to  wells  of  only  a  few  monolayers 
or  less. 
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Figure  6.4  Matrix  elements  for  intermixed  quantum  wells.  The  e,-hh,  TE  matrix 
elements  are  shown  in  (a)  and  the  e,-lh,  TM  matrix  elements  are  shown  in  (b)  for 
various  widths  of  a  GaAs/Alg jGao  7AS  well. 
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(a)  Thin-well  e-hh  transitions:  In  thin  wells,  the  wavefunctions  of  both  carriers 
are  almost  entirely  in  the  barriers,  so  the  interband  matrix  elements  of  the  well  regions 
are  nearly  insignificant.  Instead,  the  overall  matrix  element  is  dominated  by  the 
overlap  of  the  carrier  wavefunctions,  which  both  get  larger  as  the  well  intermixed  due 
to  the  heavy-hole  wavefunction  spreading  into  the  barrier  faster  than  that  of  the 
electron.  The  result  is  that  the  matrix  element  gets  larger  as  the  well  intermixes. 

(b)  Thick-well  e-hh  transitions:  For  thicker  wells,  significant  parts  of  the  carrier 
wavefunctions  are  in  the  well  region,  so  the  effect  of  the  interband  matrix  elements 
becomes  important.  As  the  well  intermixes,  aluminum  diffuses  into  the  well,  and  the 
interband  matrix  elements  in  the  well  region  decrease  in  value.  At  the  same  time,  the 
overlap  improves,  but  not  as  quickly  as  for  thin  wells,  since  the  carriers  are  more 
confined  in  the  thicker  wells.  These  competing  mechanisms  result  in  an  overall  matrix 
element  which  is  relatively  constant. 

(c)  Thin-well  e-lh  transitions:  As  with  the  thin-well  e-hh  transition,  the  carrier 
wavefunctions  are  almost  entirely  in  the  barrier,  even  more  so  since  the  light-hole 
effective  mass  is  less  than  that  of  the  heavy  hole.  Therefore,  the  interband  matrix 
element  is  dominated  by  the  barrier  and  composition  changes  in  the  well  are  inconse¬ 
quential.  The  overlap  in  this  case  is  also  nearly  constant  since  the  effective  masses  of 
the  electrons  and  light  holes  are  closely  matched,  producing  envelope  wavefunctions  of 
nearly  the  same  size  and  shape.  Since  both  mechanisms  result  in  very  small  changes, 
the  overall  matrix  element  rem<uns  nearly  constant  as  well. 

(d)  Thick-well  e-lh  transitions:  The  overlap  in  this  case  changes  only  slightly 
because  of  the  nearly  matched  electron  and  light  hole  masses  as  in  the  thinner  wells. 
But  since  significant  portions  of  the  wavefunctions  are  confined  in  the  well  region,  the 
composition  dependence  of  the  interband  matrix  elements  is  important,  as  in  the 
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thick-well  e-hh  transitions.  In  this  case,  however,  the  overlap  does  not  compete,  so  the 
interband  matrix  elements  dominate  the  overall  matrix  element.  As  the  well  intermixes, 
the  aluminum  composition  increases,  and  the  overall  matrix  element  decreases. 

6.1.2  Lifetime  Experiment 


In  order  to  test  this  theory,  an  experiment  was  devised  in  which  ultrathin  quantum 
wells  of  various  thicknesses  were  annealed,  and  the  spontaneous  emission  spectra  and 
lifetimes  were  measured.  Th-*  radiative  lifetime  is  related  to  the  Einstein  relations  as 


1  _  ^  _  8nn^n(/iv)^  ^ 

“  *  ^  mn  “  ^  B 

X  (/ic)^ 


(6.13) 


where  is  the  probability  for  spontaneous  emission,  Bmn  is  the  probability  for 


stimulated  emission,  n  is  the  index  of  refraction  (which  must  be  determined  iteratively), 
n  is  the  group  refractive  index  which  accounts  for  dispersion,  A  is  Planck's  constant,  c 
is  the  speed  of  light,  and  v  is  the  photon  energy.  Equation  (6. 13)  assumes  that  the 
probabilities  that  the  upper  electronic  state  is  empty  and  that  the  lower  state  is  filled  are 
both  unity.  The  group  index  is  given  by 


where  e  is  the  electronic  charge.  Finally,  is  given  as 


2mo€o/z^(/iv) 


(6.15) 


Several  samples  were  grown  by  molecular  beam  epitaxy  (MBE)  with  GaAs 
quantum  wells  of  varying  thickness  ranging  from  1  ML  to  8  ML.  The  wells  were  in 
the  center  of  a  1(X)0-A  AIq ^GaojAs  separate  confinement  heterostructure  bound  by 
Alo4Gao  jAs  barriers.  The  basic  structure  of  these  samples  are  shown  in  Figure  6.7. 
The  samples  were  then  annealed  for  various  times.  All  samples  were  annealed  at 
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Reure  6.5  Separate  confinement  heterostructure  for  ultrathin  quantum  wells.  Ae 
outer  barriers  are  AL^Gao  ^As,  the  inner  barriers  are  Al^  3Gao7As,  mo  the  we  i 
GaAs.  The  well  widths,  L  „  ranged  from  1  ML  to  8  Ml  in  the  lifetime  expenments 

performed  here. 
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850°C  with  one  atmosphere  AS4  overpressure,  which  resulted  in  an  interdiffiision 
coefficient  of  8  x  1 cm^s-i.  The  lifetime  of  each  sample  is  inversely  proportional  to 
the  emission  energy  and  the  matrix  element,  and  should  show  a  downward  trend 
according  to  the  theory.  This  trend  was  observed,  but  only  after  a  sharp  increase  in  the 
lifetime  in  the  first  30  minutes  of  annealing  as  seen  in  Figure  6.8.  Also,  the  calculated 
lifetimes  were  several  times  higher  than  experiment,  which  is  not  surprising  since  the 
effects  of  scattering  are  neglected  in  the  theoretical  calculation. 

As  for  the  sharp  increase  early  in  the  annealing  process,  defects  at  the  layer 
interfaces  were  most  probably  responsible.  During  the  growth  of  the  layers,  a 
10-second  stop-growth  was  used  in  order  to  smooth  out  the  layers. *21  In  order  to  stop 
the  growth,  the  shutters  to  the  source  furnaces  in  the  MBE  were  closed.  The  molecular 
beams  then  heated  the  shutters,  thereby  outgassing  impurities  which  had  adsorbed  on 
the  surface  of  the  shutters  while  they  were  cold.  These  impurities  were  then  deposited 
on  the  epitaxial  layer  where  they  generated  defects  at  the  layer  interfaces.  During  the 
first  30  minutes,  the  annealing  process  removed  some  of  the  defects  while  at  the  same 
time,  smoothing  out  much  of  the  interface  roughness  remaining  despite  the  stop-growth. 

For  one  sample  in  this  set,  all  lifetimes  were  shifted  up  by  about  200  ps.  This 
was  the  S-ML  sample,  and  the  growth  record  indicated  that  a  gallium  shutter  was  stuck 
open  during  one  of  the  stop-growths.  The  growth  conditions  were  otherwise  identical 
to  the  rest  of  the  samples  in  this  set. 

To  determine  if  the  lifetimes  of  the  wells  could  be  improved,  a  second  set  of 
samples  was  prepared  using  the  same  basic  structure  as  the  first  set.  The  only 
differences  were  that  the  structures  were  grown  without  stopping  between  the  layers, 
and  the  substrates  were  silicon-doped,  as  opposed  to  the  previous  set  which  were  semi- 
insulating.  Both  changes  were  expected  to  reduce  the  number  of  defects  at  the  layer 
interfaces  and  in  the  epitaxial  layers  as  well. 
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Rgure  6.6  Lifetimes  of  ultrathin  GaAs/Alo  jCa^ jAs  wells  as  a  function  of  annealing 
time.  The  samples  were  annealed  at  850“C  m  1  atm  AS4.  The  data  for  the  samples 
grown  with  a  stop-growth  are  presented  in  (a),  while  the  continuous-growth  sample  data 
are  plotted  in  (b).  Both  have  die  theoretical  lifetimes  for  comparison.  The  well  widths 
correspond  to  the  data  as  follows:  plus  =  1  ML,  circle  =  2  ML,  triangle  =  3  ML, 
square  =  4  ML,  star  =  5  ML,  and  asterisk  =  8  ML. 
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The  data  showed  a  slight  increase  in  the  as-grown  lifetime,  and  a  large  increase  in 
the  lifetimes  for  all  samples  annealed  for  30  minutes.  Assuming  that  the  number  of 
impurities  deposited  at  the  layer  interfaces  was  reduced  without  a  stop-growth  (some 
shutters  still  open  and  close  at  interfaces  in  order  to  change  material  composition, 
however),  the  as-grown  lifetimes  would  be  higher.  However,  the  interfaces  would  be 
rougher  in  the  second  set  of  samples,  thereby  limiting  the  effect  of  the  continuous 
growth.  After  annealing  for  30  minutes,  however,  a  dramatic  increase  in  the  lifetime 
was  observed  as  the  interface  roughness  was  smoothed,  which  supports  the  assumption 
that  there  are  fewer  impurities  deposited  at  the  interfaces.  The  lifetime  values  at  30 
minutes  are  also  in  the  vicinity  of  the  theoretical  values,  lending  further  support  to  this 
assumption. 

The  next  30  minutes  of  annealing  yielded  an  interesting  trend.  Rather  than 
following  the  curve  for  theoretical  radiative  lifetime,  the  observed  lifetimes  dropped 
sharply  and  approached  the  same  values  as  the  previous  set  of  samples. 

Several  explanations  of  this  trend  are  possible  when  it  is  examined  in  terms  of 
non-radiauve  scattering,  f^irst,  as  the  well  intermixes,  aluminum  atoms  move  into  the 
center  of  the  well  where  the  largest  part  of  the  carrier  wavefunction  can  interact  with 
the  randomly  fluctuating  atomic  potentials  characteristic  of  the  alloy.  In  other  words, 
alloy  scattering  increases,  which  in  turn  decreases  the  overall  lifetime  of  the  carriers. 

While  alloy  scattering  is  certainly  a  contributor  to  the  trend  seen  in  Figure  6.6,  it 
is  probably  not  the  dominant  one.  If  it  were,  the  alloy  scattering  for  the  thinner  wells 
would  saturate  sooner  than  in  the  thicker  wells  since  most  of  their  wavefunctions  are  in 
barriers,  and  the  wells  would  fill  with  aluminum  atoms  faster.  Such  a  trend  is  not 
evident  in  the  data. 
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Another  explanation  based  on  non-radiative  scattering  that  may  explain  the  data  is 
that,  after  about  an  hour,  other  types  of  defects  diffuse  into  the  area  of  the  quantum 
wells.  These  defects  could  be  gallium  vacancies  (VqJ  or  arsenic  antisite  defects  (ASqJ 
coming  from  the  surface  of  the  wafer,  due  to  the  AS4  overpressure  during  the  annealing 
process.  Guido  et  al.  showed  that  the  quantum  well  intermixing  rate  was  dependent  on 
the  depth  of  the  well  during  an  annealing  step  in  an  arsenic  overpressure.  This 
indicated  that  the  transport  time  of  the  defects  from  the  surface  had  a  role  in 
determining  the  Al-Ga  interdiffusion  coefficient.  They  determined  that  the  diffusion 
coefficients  for  the  defects,  D and  D  .  .  ,  were  5  x  10"'^  cm^s  *  at  825 “C  with  a 
2-atm  arsenic  overpressure.  The  depth  to  which  the  defects  diffused  in  30  minutes  is 
therefore  on  the  order  of  1  ^m  (10,000  A),  which  is  within  a  factor  of  two  of  the  depth 
of  the  structures  in  this  study.  Similarly,  they  could  be  defects  diffusing  upward  from 
the  substrate.  Further  experiments  are  required  to  support  this  explanation.  For 
instance,  the  well  could  be  buried  at  various  depths  from  the  surface  to  determine  when 
the  sudden  drop  in  lifetime  occurs.  If  the  drop  occurs  after  longer  annealing  for  deeper 
wells,  then  the  explanation  of  defects  coming  from  the  surface  would  be  supported. 

A  third  explanation  of  the  drop  is  possible  by  looking  at  the  carrier  dynamics 
between  the  bound  states  and  the  barrier,  and  then  comparing  their  respective  lifetimes. 
If,  after  an  hour  of  annealing,  the  barrier  became  a  sink  for  carriers,  then  the  sudden 
drop  in  lifetime  would  be  expected.  In  order  for  this  to  happen,  the  barrier  lifetime 
would  have  to  be  much  shorter  than  that  of  the  well,  and  there  would  have  to  be  a 
scattering  mechanism  present  to  kick  the  carriers  to  the  barrier's  higher  energy  level. 
This  explanation  of  thermal  scattering  is  plausible  for  1-  or  2-ML  (and  possibly  3-ML) 
wells,  but  for  larger  wells,  the  separation  between  the  barrier  and  bound-state  energies 
is  several  times  the  thermal  energy  at  77K.  Therefore,  another  scattering  mechanism 
must  be  considered. 
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Another  problem  with  this  third  explanation  is  that  if  the  short  lifetime  in  the 
barriers  was  radiative,  the  barrier  photoluminescence  peak  would  be  dominant,  yet  the 
well  peak  was  always  much  stronger  than  that  of  the  barrier.  If  the  short  barrier 
lifetime  was  due  to  non-radiative  mechanisms,  then  the  carriers  would  never  get  to  the 
bound  state  and  the  photoluminescence  peak  at  the  bound  state  energy  would  be  weak, 
when  in  fact  ^he  peak  remained  strong  throughout  the  set  of  annealing  points. 

The  most  probable  explanation  of  the  sudden  drop  in  the  lifetime  at  about  one 
hour  of  annealing  is  the  second  one  given  above.  The  transport  of  defects  to  the 
quantum  well  region  probably  yielded  the  data  trend  observed.  This  theory  can  be 
tested  by  performing  a  series  of  lifetime  measurements  on  annealed  samples  with 
varying  cladding  layer  thickness. 


6.2  Damping  Coefficient 


The  last  variable  needed  to  describe  the  dielectric  function  is  the  damping 
coefficient,  which  is  defined  by  the  spectral  line  width  of  the  oscillator  as2Y  =  Acjo, 
where  Acxiis  the  full-width-at-half-maximum  (FWHM)  of  the  peak  at  the  resonant 
frequency. 

The  line-width  broadening  is  due  to  the  finite  lifetimes  of  the  electronic  states  and 
is  related  to  the  lifetime  through  the  Heisenberg  Uncertainty  Principle  ft  A  co  t  «  ft.f’^i  or 
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If  the  lifetime  is  radiative  only,  then  it  is  referred  to  as  a  natural  broadening.  This 
treatment  neglects  the  effects  of  carrier  scattering  which  may  be  important  since  the 
interaction  time  of  collisions  can  be  several  orders  of  magnitude  lower  than  the  mean 
radiative  lifetimes.  Using  Equations  (6. 13)  to  (6. 16),  the  damping  coefficient  due  to 
natural  broadening  is  then 
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(6.17) 


The  oscillators  over  which  the  sum  is  taken  in  Equation  (6.6)  are  transitions 
between  the  electronic  states  in  the  system.  In  a  quantum  well,  each  energy  level  has  a 
corresponding  density  of  electronic  states,  so  if  the  dielectric  response  of  each  level  is 
multiplied  by  the  density  of  states  (DOS)  function  and  integrated  over  the  energy,  the 
dielectric  constant  as  a  function  of  energy  become  s 
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where  p  ( £ '  -  Eg)  is  the  reduced  DOS,  r„„  =  ft  Y  nm.  and  the  sum  is  now  over  the 


pairs  of  energy  levels  n  and  m  between  the  conduction  and  valence  bands, 
respectively.  A  distinction  is  made  here  between  energy  states  and  energy  levels, 
where  each  level  is  a  solution  to  the  Schrodinger  Wave  Equation  in  the  quantum  well, 
and  contains  a  number  of  states  described  by  the  DOS.  Normally,  the  DOS  for  a 
quantum  well  is  taken  as  a  constant  in  energy,  but  the  three-dimensional  (3D)  DOS 
can  be  used  with  the  value  of  the  bound  energy  level  of  the  well  used  in  the  energy 
dependence  (see  Appendix  C).  To  simplify  the  problem,  the  integral  is  taken  over  the 
DOS.  By  doing  this,  the  number  of  oscillators  contributing  to  each  pair  of  energy 
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levels  is  determined,  which  allows  the  equation  to  retain  the  general  form  of 
Equation  (6.6).  The  lineshape  is  broadened  by  r„„„  which  means  the  energy  levels 
are  actually  spread  into  bands,  but  this  approximation  keeps  the  energy  dependence  in 
the  dielectric  response  while  greatly  simplifying  the  calculations.  The  dielectric 
constant  is  then  approximated  by 
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where  N  „„  is  the  number  of  oscillators  with  the  resonant  energy  E  The  function 


p  is  the  2D  DOS  function  derived  in  Appendix  C  without  spin-degeneracy,  and  is 
given  by 
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where  p  is  the  reduced  mass  of  the  combined  carriers  in  the  conduction  and  valence 


bands,  which  will  be  discussed  shortly,  and  A:„m  is  the  wavevector  for  the  energy  levels 
given  by 
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The  energy  gap,  £  is  between  the  bottom  of  the  conduction  band  potential  well  and 


the  bottom  of  the  valence  band  well  (highest  absolute  energy). 
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Since  the  material  composition  varies  across  the  intermixed  well,  the  effective 
mass  should  be  the  average  value  of  the  material  at  the  bottom  of  the  well  and  the 
barriers.  This  can  be  handled  in  the  same  manner  as  the  interband  matrix  element, 
using  the  carrier  wavefunction  to  weight  the  values  of  the  mass  for  each  region.  In  this 
case,  however,  only  the  individual  wavefunctions  are  used  instead  of  the  overlap.  For 
example, 

‘"0  •'*  (6.23) 

and 


11^1 

^  (6.24) 

Finally,  for  a  superlattice,  the  DOS  is  also  given  by  Equation  (6.21),  exceptA:„m 

is  replaced  by  the  superlattice  wavevector  q„m,  which  is  given  by 


<7«m  =  ^C0S‘‘(F(£)) 
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(6.25) 


since  the  superlattice  dispersion  relation  is  described  by 

cos(ga)  =  f(£)  ,  (6.26) 

where  a  is  the  period  of  the  superlattice,  and  £(£)  is  a  function  of  energy.  For  the 
intermixed  superlattice,  £(£)  takes  the  form,  previously  derived  in  Section  4.3,  of 


^(f)  =  ^(5n-^S22) 


(6.27) 


where  S 1 1  and  S  22  are  elements  of  a  transfer  matrix  that  describes  the  superlattice 


period.  The  DOS  for  an  intermixed  superlattice  is  then 
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6.4  Index  of  Refraction 


The  index  of  refraction  of  an  intermixed  quantum  well  as  function  of  photon 
energy,  E  =  huj,  is 


=  Re 


(6.29) 


where  g  ( f )  is  given  by  Equation  (5.20).  The  background  dielectric  function  defined 

earlier  includes  the  effects  of  the  barrier  oscillators  which  do  not  change  as  the  quantum 
well  is  annealed.  If  the  frequency  range  of  interest  is  far  enough  away  from  the 
spectral  features  of  the  background  diel^tiic  constant,  then  the  background  constant 
will  be  real.  In  the  case  of  a  quantum  well,  the  background  dielectric  constant  can  be 
taken  from  the  barrier  material  using  an  empirically  determined  index  of  refraction 
given  byf^^i 


(6.30) 


where  x  is  the  aluminum  fraction  in  the  Al,jGa,.,As  alloy,  A.  „  is  the  free-space 
wavelength  in  microns,  and 

A=  10. 9060-02. 92x 
5  =  0.97501 

C  =  (0.52886 -0.735x)^  ;  x<0.36 
C  =  (0.30386 -0.105x)^  ;  x>0.36 
D  =  0.002467(l  +  1.41x)  . 


The  background  dielectric  constant  is  then  given  by  Gj,  =  and  is  a  function  of 


or  5. 
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A  potential  problem  with  this  method  occurs  when  the  bound  states  are  near  the 
top  of  the  well,  as  with  a  1-  or  2-ML  well,  where  £5  may  be  complex  and  can  no  longer 
be  described  by  the  index  of  refraction.  However,  if  the  spectral  range  of  interest  for 
passive  waveguides  is  far  enough  away  from  the  absorption  edge  of  the  barrier,  this 
method  should  still  be  valid.  For  a  passive  waveguide  to  have  acceptable  losses,  the 
frequency  of  the  guided  light  must  be  far  enough  below  the  absorption  edge  of  the 
bound  state,  which  therefore  places  the  frequency  below  the  absorption  edge  of  the 
barrier.  As  long  as  the  dielectric  properties  in  the  barrier  remain  unchanged  with 
annealing,  all  changes  at  the  desired  frequency  will  be  due  to  changes  in  the  bound 
state. 

The  refractive  index  of  a  5-ML  GaAs/Al<j  jCa^  7AS  quantum  well  as  a  function  of 
energy  for  several  annealing  times  was  calculated  and  the  results  are  presented  in 
Figure  6.9.  Only  the  effects  of  n=m=l  e-hh  transition  are  shown  since  their  influence 
is  the  greatest  for  the  energies  of  interest. 

The  abscissa  covers  the  energies  between  the  bottom  of  the  as-grown  well  and  the 
barrier  with  the  lowest-energy  e-hh  transition  of  the  as-grown  well  labeled.  As  the  well 
intermixes,  the  transition  energy  increases  since  the  bottom  of  the  well  fills  in  as  the 
aluminum  composition  increases.  The  widths  of  the  peaks  at  the  resonance  points  are  a 
function  of  the  damping  coefficient  for  each  discrete  energy  level  as  calculated  above. 
This  is  for  the  ideal  case  where  the  energy  level  is  not  spread  out  at  all  and  the  damping 
coefficient  is  due  to  natural  linewidth  broadening  only.  The  effects  of  thermal 
broadening  is  shown  in  Figure  6.7  (b)  where  the  damping  coefficient,  f  u.  is  replaced 
with  the  thermal  energy  at  77K  (6  meV).  The  effects  this  has  on  the  index  are  isolated 
near  the  resonance,  so  as  long  as  this  region  is  avoided,  the  actual  index  will  be  well 
described  by  the  ideal  theoretical  case. 
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Figure  6.7  Refractive  index  for  an  intermixed  GaAs/Alo  3Ga(,  7As  quantum  well.  The 
index  for  a  S  ML  well  as  grown  (solid  curve)  and  annealeo  for  30  minutes  (dot-dash) 
and  120  minutes  (short  dashed)  is  shown  in  (a)  for  the  energy  range  between  the  GaAs 
and  barrier  band  edges.  The  effects  of  thermal  lineshape  broadening  are  included  in  (b) 
to  point  out  that  the  index  does  not  change  except  in  the  immediate  vicinity  of  the 
oscillator. 
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Another  trend  to  note  is  the  decrease  in  the  strength  of  the  resonance  as  it 
approaches  the  barrier  energy.  While  this  is  difficult  to  quantify  in  Figure  6.7  (a)  due 
to  the  selection  of  plotting  points  near  the  sharp  resonances,  the  thermally-broadened 
curves  in  Figure  6.7  (b)  show  this  well.  This  makes  intuitive  sense  because  the  index 
will  eventually  become  that  of  the  barrier  as  the  well  becomes  completely  intermixed. 
This  is  probably  not  a  function  of  oscillator  strength  because  the  oscillator  in  the 
denominator  increases  only  slightly  with  annealing  time,  which  is  offset  by  the 
increased  matrix  element  in  the  numerator.  Instead,  the  trend  is  due  to  the  decrease  in 
the  density  of  oscillators  at  the  resonant  frequency  as  described  in  Equation  (6.20). 

This  number  is  proportional  to  the  difference  between  the  energy  level  and  the  bottom 
of  the  well  which  decreases  as  the  well  fills  in. 

In  this  particular  example,  if  the  emission  energy  of  GaAs  (1.508  eV)  is  assumed, 
then  the  index  of  refraction  for  a  30-minute  annealing  step  is  3.4231,  and  a  120-minute 
annealing  step  results  in  an  index  of  3.4217.  While  these  values  represent  different 
times  t  with  the  same  interdiffusion  coefficient  A  they  could  also  represent  the  same 
annealing  times  with  different  interdiffusion  coefficients  since  the  amount  of 
intermixing  depends  on  the  product  Dt.  In  that  case,  this  example  could  be  used  to 
represent  a  selective  intermixing  process  with  a  4: 1  selectivity  ratio.  The  index 
difference  between  the  "preserved"  and  "intermixed"  regions  would  then  be  0.0014 
(0.04%),  which  is  less  than  half  that  of  the  ideal  case  of  0.0034  (0. 1  %)  where  the  well 
is  either  completely  preserved  or  completely  intermixed. 

By  selectively  intermixing  a  5-ML  quantum  well  with  a  4: 1  selective 
interdiffusion  coefficient  ratio,  an  optical  waveguide  can  be  formed.  Using  the  data 
above,  the  electric  field  distribution  for  a  5-?4m-wide  waveguide  would  look  like  that 
shown  in  Figure  6. 10.  Such  a  waveguide  could  be  useful  in  a  directional  coupler  where 
two  waveguides  are  placed  in  close  proximity  so  that  the  electric  field  modes  overlap 
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Figure  6.8  Electric  field  distribution  for  an  optical  waveguide  defined  by  selective 
intermixing.  The  indices  of  refraction  are  based  on  a  4: 1  selectivity  ratio  of  the 
interdiffusion  coefficients  in  the  cladding  and  guide  regions.  The  as-grown  quantum 
well  is  5  ML  wide  and  the  waveguide  is  5  #tm  wide. 


and  power  is  transferred  between  guides.  The  coupling  length  for  such  a  device,  which 
is  the  distance  that  the  optical  field  must  propagate  in  order  to  transfer  all  power  from 
one  guide  to  the  other,  is  given  by 


n  ^  mn 
2k,  k. 


(6.31) 


where  k,  is  the  coupling  coefficient,  and  mis  the  mode  number  for  the  transverse 


electric  (TEJ  waveguide  mode.  The  coupling  coefficient  is  defined  by  the  overlap  of 
one  waveguide  mode  with  the  evanescent  tail  of  the  adjacent  waveguide  mode.  1*^1  If  the 
overlap  is  small  enough  to  be  treated  as  a  perturbation,  then  the  coupling  coefficient  for 
a  coupler  with  propagation  in  the  x-direction  and  guide  separation  in  the  y-direction  is 
given  byt»^i 

K  rt 

+  ’  (6.32) 

where  Py  is  the  mode  propagation  coefficient  in  the  y-direction,  p ,,  is  the  propagation 

coefficient  in  the  x-direction,  y  „  is  the  1  /  e  point  of  the  evanescent  field,  W  is  the  width 
of  the  waveguides,  and  s  is  the  separation  between  them.  The  length  of  a  3  dB  coupler, 
which  splits  the  power  evenly  between  the  waveguides,  for  a  TE„  mode  is  shown  in 
Figure  6.11. 

For  a  reasonable  5-%m  separation,  the  splitter  would  need  to  be  310  %m  long  to 
transfer  50%  of  the  optical  power  between  waveguides,  which  is  also  reasonable  for 
integrated  optics.  Shorter  lengths  can  be  easily  obtained,  but  the  control  over  the  length 
becomes  more  difficult,  and  variations  in  length  become  more  critical  to  the  coupling 
ratio.  This  is  just  a  single  example  of  how  selectively  intermixed  quantum  wells  can  be 
used  in  integrated  optics. 
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Figure  6.9  Interaction  length  of  a  3  d6  optical  waveguide  coupler  designed  with 
disorder-defined  waveguides.  The  waveguides  are  comprised  of  selectively 
intermixed  5  ML  GaAs/Alg  jGaQ  7AS  quantum  wells.  The  waveguide  widths  are 
indicated  in  the  legend,  and  the  interaction  length  for  50%  power  transfer 
between  the  guides  is  given  as  a  function  of  guide  separation  for  each  case.  The 
inset  shows  a  schematic  of  the  device. 
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7  The  Trek 


This  final  chapter  is  a  review  of  this  research  endeavor.  It  begins  with  a  summary 
of  the  work  performed.  The  next  section  outlines  the  specific  conclusions  drawn  from 
the  various  experiments  and  calculations  presented  earlier.  Finally,  recommendations 
are  made  to  the  enterprising  engineer  who  would  use  this  thesis  as  a  starting  point  for 
designing  integrated  optical  circuits  based  on  thermally  intermixed  quantum  wells  and 
superlattices. 

7.1  Voyages 

This  research  has  focused  on  quantifying  the  changes  of  the  optical  properties  of 
quantum  wells  and  superlattices  due  to  the  thermal  interdiffusion  of  the  heterostructure 
interfaces.  Various  interdiffusion  experiments  were  performed  and  characterized  using 
photoluminescence  (PL)  techniques  and  theoretical  models  in  order  to  develop  a  new 
model  that  describes  the  changes.  A  summary  of  this  study  is  presented  here, 
emphasizing  the  original  research  performed. 

The  potential  well  shape  resulting  from  heterostructure  intermixing  was  examined 
in  Chapter  4.  It  started  with  the  refinement  of  a  computer  model  to  include  an  error 
function  shape  in  studying  intermixed  heterostructure  interfaces  (Section  4.1).  Further 
improvements  in  the  model's  speed  were  achieved  with  a  thoughtful  scheme  of  dividing 
the  intermixed  well.  The  scheme  presented  yielded  greater  accuracy  for  the  same 
number  of  divisions.  The  validity  of  the  error  function  composition  profile  was  then 
questioned,  and  an  experiment  was  performed  to  determine  if  the  Al-Ga  interdiffusion 
coefficient  is  dependent  on  aluminum  composition  (Section  4.2).  A  computer  model 
was  developed  which  used  a  composition-dependent  interdiffusion  coefficient,  and 
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simulations  were  compared  to  the  experimental  results.  While  a  small  compositionally 
dependent  trend  was  observed,  the  differences  were  within  the  experimental  error  and 
deemed  insignificant  within  the  conditions  of  the  experiment. 

The  next  section  examined  models  to  calculate  the  energy  dispersion  relations  of 
superlattices  (Section  4.3).  A  one-dimensional  tight-binding  calculation  in  the 
envelope-function  approximation  was  studied  first,  but  this  method  was  limited  to 
weakly  coupled  wells.  A  Kronig-Penney  type  calculation  was  examined  next.  This 
model  was  able  to  handle  strongly  coupled  wells,  and  was  successfully  adopted  to 
thermally  intermixed  superlattices  using  error-function-shaped  interfaces. 

Chapter  5  addressed  several  potential  techniques  to  selectively  control  the 
interdiffusion  process.  First,  the  use  of  gallium  ions  implanted  near  the  surface  to 
create  a  local  gallium  overpressure  was  studied  (Section  S.l).  A  gallium  focused  ion 
beam  (FIB)  was  used  to  selectively  implant  ions  into  the  surface  near  (not  in)  a 
GaAs/AlGaAs  superlattice,  and  the  samples  were  annealed.  PL  was  used  to  measure 
the  changes  in  energy  with  and  without  the  implanted  ions.  The  results  suggested  that 
the  implant  generated  mobile  defects  which  served  to  preserve  the  superlattice  in  and 
around  the  implanted  areas.  This  led  to  a  side  experiment  to  ascertain  the  effects  of 
different  ambient  overpressures  on  the  mobility  of  electrons  in  a  high-electron-mobility 
transistor  (HEMT)  structure.  The  nature  of  the  structure,  however,  prevented  the 
conclusive  determination  of  the  effects  of  the  overpressure.  A  more  appropriate 
experiment  was  proposed  based  on  the  availability  of  an  FIB  apparatus. 

The  next  method  of  selective  intermixing  examined  was  the  used  of  SiN,  as  an 
encapsulant  to  preserve  the  morphology  of  the  surface  (Section  5.2).  A  GaAs/AlAs 
superlattice  was  prepared,  and  SiN,  stripes  were  patterned  on  the  surface.  The  sample 
was  annealed  and  the  PL  peak  shifts  were  measured.  This  method  provided  a  suitable 
selectivity  between  the  covered  and  uncovered  stripes. 
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Finally,  arsenic-rich  GaAs  grown  at  low  substrate  temperature  (LT-GaAs)  was 
considered  as  a  source  of  local  arsenic  overpressure  (Section  5.3).  An  LT-GaAs  layer 
was  grown  and  patterned  on  a  GaAs  substrate,  then  reloaded  into  the  molecular  beam 
epitaxy  machine  for  subsequent  layer  growth.  The  sample  was  then  annealed  and  the 
PL  peak  shifts  showed  a  trend  verifying  that  the  LT-GaAs  stripes  enhanced  the  Al-Ga 
interdiffusion  coefficient.  Although  the  processing  techniques  need  refinement,  this 
proof-of-concept  experiment  demonstrated  that  this  is  a  possible  method  of  selective 
intermixing. 

The  primary  emphasis  of  this  thesis  was  chronicled  in  Chapter  6.  After  a  brief 
overview  of  the  relationships  between  the  various  optical  properties  and  origin  of  the 
dielectric  constant,  a  model  for  the  optical  matrix  element  was  developed  (Section  6.1). 
The  matrix  element  based  on  the  Kane  model  was  discussed,  and  the  theory  was 
extended  to  include  the  effects  of  the  electron-hole  envelope  wavefunctions  in  thin 
wells,  as  well  as  the  effects  of  varying  interband  matrix  elements  between  the  well  and 
barriers.  The  new  model  was  then  applied  to  thermally  intermixed  quantum  wells. 
Simulations  based  on  the  models  were  performed  and  the  radiative  electron-hole 
recombination  lifetimes  were  calculated.  These  simulations  were  compared  to 
time-resolved  PL  of  thermally  annealed  thin  wells.  The  data,  however,  was  masked  by 
what  was  probably  non-radiative  defects  that  were  transported  into  the  structure  from 
the  surface  during  the  annealing  process. 

The  damping  coefficient  was  calculated  next  using  a  simple  model  describing  the 
natural  linewidth  broadening  (Section  6.2).  The  derivation  of  the  dielectric  constant  for 
intermixed  wells  was  then  completed  using  a  modified  density  of  states  function 
(Section  6.3).  The  density  of  states  function  for  a  superlattice  was  then  formulated 
based  on  a  previous  calculation  of  the  energy  dispersion  relation. 
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The  chapter  culminated  with  the  first-ever  calculation  of  the  index  of  refraction  of 
a  partially  intermixed  quantum  well  (Section  6.4).  A  5-monolayer  (ML)  well  was  used 
as  an  example,  and  the  index  difference  was  calculated  for  a  selective  intermixing  case. 
The  example  was  extended  by  including  the  intermixed  wells  in  an  optical  waveguide, 
calculating  the  electric  field  distribution,  and  modeling  a  3-dB  power  splitter.  The 
calculation  was  presented  as  an  example  of  using  realistic  values  of  indices  of  refraction 
of  partially  intermixed  wells  in  a  useful  integrated  optical  device. 

7.2  Where  No  One  Has  Gone  Before 

The  following  is  list  of  conclusions  which  can  be  drawn  from  the  work  performed 
in  this  thesis. 

(a)  Compositional  dependence  of  the  Al-Ga  interdiffusion  coefficient  is  negligible  at 

850®C  with  1  atm  AS4  overpressure.  While  a  slight  dependence  is  evident,  it  is 
within  the  experimental  error.  Higher  temperatures,  however,  may  exhibit  a 
stronger  dependence. 

(b)  The  accuracy  of  a  numerical  solution  of  the  Schrodinger  Wave  Equation  in  an 

error-function-shaped  potential  well  is  improved  if  the  finite  regions  of  the  well 
are  chosen  such  that  there  are  more  regions  in  areas  of  rapidly  varying  potential. 

(c)  Energy  levels  of  quantum  wells  with  widths  between  8  and  10  ML  (20  and  30  A) 

are  the  most  sensitive  to  thermal  annealing.  The  energy  levels  of  thinner  wells 
tend  to  saturate  near  the  barrier,  and  levels  in  thicker  wells  rely  on  the  well  to  fill 
in  before  they  shift  in  energy. 

(d)  SiN„  can  be  used  to  selectively  mask  areas  of  a  quantum  well  or  superlattice  from 

overpressure-enhanced  intermixing.  The  minimum  feature  sized,  however,  has 
yet  to  be  determined. 
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(e)  Patterned  GaAs  grown  at  low  substrate  temperature  (LT-GaAs)  can  be  used  as  a 

localized  "arsenic  oveipressure"  to  enhance  intermixing.  Additional  work  is 
needed  to  improve  the  fabrication  process  in  order  to  determine  how  much 
interdiffusion  enhancement  is  possible. 

(f)  The  electron-to-heavy-hole  (e-hh)  envelope  wavefunction  overlap  is  critical  to  the 

optical  matrix  element  of  thin  quantum  wells.  The  interband  matrix  elements  of 
both  the  well  and  the  barriers  are  also  important  when  significant  portions  of  the 
wavefunctions  are  in  each. 

(g)  The  optical  matrix  element  for  e-hh  transitions  in  ultrathin  quantum  wells  gets 

larger  as  the  wells  are  thermally  intermixed.  This  is  primarily  due  to  the 
improvement  in  the  carrier  wavefunction  overlap.  This  effect  also  translates  to  a 
theoretical  reduction  in  the  radiative  lifetime  of  the  wells. 

(h)  The  refractive  indices  of  a  selectively  intermixed  quantum  well  are  somewhere 

between  that  of  a  completely  intermixed  well  and  a  completely  preserved  well. 

The  index  difference  between  the  preserved  and  intermixed  areas  is  then  much 
less  than  between  the  two  extremes,  which  is  important  when  applied  to  the  design 
of  a  device,  such  as  an  optical  waveguide. 

The  overall  conclusion  of  this  thesis  is  that  the  optical  properties  of  quantum  wells 
and  superlattices  can  be  altered  through  selective  intermixing,  and  the  resulting 
differences  are  large  enough  to  be  used  in  integrated  optical  applications.  This  work 
quantifies,  for  the  first  time,  these  Cuanges  in  relation  to  the  extent  of  intermixing  the 
quantum  well  experiences.  While  there  are  several  more  elaborate  theories  that 
describe  the  optical  properties  of  both  quantum  wells  and  bulk  material,  this  work 
bridges  the  gap  between  them,  which  is  essential  if  selective  intermixing  is  to  become  a 
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useful  technology.  From  here,  then,  design  rules  can  begin  to  emerge  which  govern 
the  fabrication  of  integrated  optical  components  that  take  advantage  of  the  unique 
properties  of  quantum  wells  and  superlatdces. 

7.3  The  Continuing  Mission 

While  this  thesis  is  a  significant  advancement  in  the  science  of  the  optical 
properties  of  semiconductor  quantum  structures,  more  work  is  required  to  make  it  a 
complete  technology.  This  section  contains  recommendations  which  well  build  on  this 
thesis  and  approach  such  a  level.  Explicit  suggestions  for  individual  experiments  to 
improve  the  work  can  be  found  in  the  body  of  this  document,  while  this  section  is 
reserved  for  more  general  recommendations. 

Two  recommendations  are  made  from  an  experimental  standpoint.  The  first  is  to 
find  an  optimum  method  of  selectively  controlling  the  intermixing  process.  The 
techniques  tested  here  are  just  a  few  of  the  possible  candidates.  While  SiN, 
encapsulation  shows  the  most  promise,  its  limits  need  to  be  tested,  and  other  methods 
should  be  considered,  or  perhaps,  combined  for  the  best  results. 

The  second  experimental  recommendation  is  to  fabricate  some  integrated  optical 
devices,  such  as  the  one  described  in  Section  5.4,  in  order  to  test  the  theory.  This,  of 
course,  applies  to  any  extensions  of  this  theory  as  well. 

On  the  theoretical  side,  other  effects  can  be  factored  into  the  model.  These  may 
include  the  effects  of  X  and  L  band  transitions  (especially  at  higher  aluminum 
composition),  the  inclusion  of  non-radiative  scattering  processes  (such  as  phonons)  in 
the  calculation  of  the  lifetime  and  damping  coefficient,  or  excitons  for  low-temperature, 
low-carrier-excitation  applications. 
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Another  possibility  is  to  extend  the  basic  theory  to  include  electro-optic  effects. 
These  effects  are  useful  in  such  devices  as  integrated  optical  interferometers,  bistable 
devices,  spatial  light  modulators,  and  self-electro-optic  devices  (SEED’S).  This 
particular  recommendation  would  be  especially  well-complemented  by  appropriate 
experimental  investigations. 

These  are  just  a  few  .suggestions  of  research  that  builds  on  this  study.  This 
dissertation  should  serve  as  a  good  starting  point  for  an  investigation  of  any  of  these,  or 
any  other  studies  which  are  concerned  with  the  optical  properties  of  intermixed 
quantum  wells  and  superlattices. 
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9.1  Appendix  A 
The  Dielectric  Constant 


The  optical  properties  of  a  semiconductor  are  contained  in  the  complex  dielectric 
constant 

e(oo)  =  e'(co)  + 

Absorption  and  reflection  properties  can  be  found  by  relating  e( go)  to  the  complex 
index  of  refraction  by 


^^^  =  n£(GO)  =  (^(GO)  +  iK(GO))^  . 


(9.2) 


where  n  (  go  )  is  the  real  index  and  k  (  go  )  is  the  extinction  coefficient.  The  absorption 
coefficient  a  (  go  )  is  related  to  the  extinction  coefficient  by 


,  ^  2ic(a>) 

a(Go)  = - -uo 


(9.3) 


where  c  is  the  speed  of  light.  Using  e"/€o  =  2nK,  a(Go)is  related  to  the  imaginary 


part  of  e( go)  by 


a(cw) 


GO 

eoCn(Go) 


e-(uj)  . 


Finally,  the  reflectivity  at  normal  incidence  is  given  by 

'^“(n(a))^l)2  +  (K(u)))" 


(9.4) 


(9.5) 


The  usual  method  of  experimentally  determining  the  optical  constants  is  to 
measure  the  absorption  spectrum  (and  hence  the  extinction  coefficient,  ic(go))  or  the 
reflectivity  (which  gives  n(  go  )  for  small  k(go  ))  for  a  wide  frequency  range,  and  then 
use  the  Kramers-Kronig  relations  to  transform  the  information  back  to  the  index  of 
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refraction  or  extinction  coefficient,  respectively.  Similarly,  if  the  imaginary  part  of  the 
dielectric  constant  is  calculated,  the  real  part  can  be  found  via  the  Kramers-Kronig 
relations,  and  vice  versa.  The  transformations  are  given  by 


‘  (A) '  e  "  ( CO ' )  d  GO ' 


nJo 

2<o  /‘-[e'(co')-eo]d(o' 


(o'^-to^ 


e"((o)=-^  r 


'2  2 
O)  -  (A) 


9.1.1  The  Oscillator  Model 


(9.6) 


More  fundamental  than  the  dielectric  constant  is  the  polarizability,  which  is 
described  by  the  atomic  dipole,  p.  The  classical  equation  of  motion  for  an  electron 
oscillating  against  a  fixed  nucleus  is 

mgX  +  2ymgX-*-  Kx^-e£cos(X)t  (9.7) 

where  mg  is  the  mass  of  an  electron  oscillating  against  the  nucleus,  x  is  the 

displacement,  y  is  a  frictional  coefficient  due  to  dissipative  forces  such  as  phonons,  K  is 
the  spring  constant,  eis  the  electronic  charge,  and  Ceos  cot  is  the  electric  field  from 
the  photons.  The  solution  to  the  displacement  is 

eS/m„ 

X  - - ; - ; - 

(u):-(A)^)-2iu)Y  (9g) 

where  oog  -  yjK/rUg  is  the  resonant  frequency  of  the  oscillator.  The  polarization  is 


then 


e^/m. 


p  =  -ex  = 


•£  . 


(ou‘-u>^)-2icoY  (9.9) 

The  macroscopic  polarization  P  is  the  sum  of  the  atomic  dipoles  per  unit  volume,  so  for 


the  single  oscillator,  the  dielectric  constant  is 
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e(u))-e„  +  = 
£' 


TTioV  ujf-u)^-2/a)Y 

_ 2yoo _ 

*  (o),  -  t«^)^  +  4Y^t*)^  mpK  ((A)o  -  00^)*  +  4Y^tw^ 

=  e'(u))+ie"(tA))  .  (9.10) 


For  CO  »  Y,  the  imaginary  part,  g  "  ( lo  ),  is  non-zero  only  for  co  near  co  In  that  case, 
it  can  be  approximated  by  a  Lorentzian  curve  with  full-width-half-max  of2Y, 

Y 


G"(CO)* 


2m„uj„V  (oOo-(o)^+Y' 


(9.11) 


If  Y  approaches  zero,  g  "  ( co  )  can  be  approximated  by  a  delta  function 


G"((0) 


e^n 


2moCOpK 


6(to  -CO) 


2mo(OoV' 


6(/i(Oo-^co) 


(9.12) 


Note  that  the  Kramers-Kronig  transformation  of  Equation  (9. 12)  gives  the  real  part  of 
Equation  (9. 10)  with  y  =  0.  This  transformation  is  unnecessary  in  this  case,  however, 
since  all  the  parameters  used  to  calculate  g  "  ( co )  can  be  used  to  calculate  £ '  ( (^ ) 
directly. 


9.1,2  Quantum  Mechanical  Model 

The  result  above  is  only  for  a  single  classical  oscillator,  and  has  been  simplified 
by  removing  friction  losses,  which  is  analogous  to  the  conservation  of  energy.  In  a 
semiconductor,  however,  many  oscillators  with  various  frequencies,  damping 
coefficients  and  strengths  are  present,  and  the  classical  case  breaks  down.  In  order  to 
properly  determine  the  dielectric  constant,  a  quantum  mechanical  description  is  needed. 
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The  derivation  begins  with  the  calculation  of  the  absorption  coefficient,  which  is  then 
related  to  €  "  ( ^  )  by  Equation  (9.4).  This  begins  with  the  absorption  and  emission 
rates  and  respectively,  which  are 


1/ 


nmPX  2 


mlvo^V 


(9.13) 


The  subscripts  denote  the  following:  mis  the  upper  conduction  band,  n  is  the  lower 
valence  band,  (3  is  the  photon  wavevector,  \  is  one  of  two  possible  polarizations,  and  k 
is  the  electronic  wavevector  whic  defines  the  particular  state  in  the  band  mor  n. 

The  factor  nxp  is  the  photon  occupation  number  given  by  the  Bose-Einstein 


distribution  function.  The  absorption  rate  and  the  first  term  in  the  emission  rate  are  the 
stimulated  terms  which  depend  on  the  number  of  photons  n  ^6,  while  the  other  emission 
term  is  independent  of  then^p,  and  hence  is  the  spontaneous  emission  term.  The  first 
delta  function  maintains  the  conservation  of  momentum,  assuming  negligible  photon 
momentum,  and  the  second  delta  function  is  the  conservation  of  energy  from  Fermi's 
Golden  Rule.  Finally,  P„mx  is  the  momentum  matrix  element,  which  is  assumed  to  be 
independent  of  fc,  given  by 


nmX 


=  -m  J  u.\ 


(r)<5x*V,u„,(r)d  r 


(9.14) 


In  Equation  (9.14)  the  dipole  approximation  is  assumed,  and  is  the  unit  vector 


in  the  direction  of  the  electric  field  of  the?v-polarized  photons.  The  functions  u(  r )  are 
the  periodic  parts  of  the  wavefunctions,  which  are  usually  taken  near  the  zone  center 
with  a  small  range  of  k. 
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The  momentum  matrix  elements  have  been  worked  outi^i  for  two  conduction  band 
states  (spin-up  and  spin-down)  and  six  valence  band  states  (hh,  Ih,  so;  up  and  down) 
including  spin-orbit  interaction.  The  values  are  given  in  Table  9.1. 


Table  9.1.  Momentum  Matrix  Elements  between  Conduction-  and 
Valence-Band  States 


Cl 

C2 

hhl 

0 

hh2 

0 

i^(X,-iXy) 

Ihl 

lh2 

P  0  ..  ..  . 
^(Xjj  -  iXy) 

sol 

^/3  ' 

^(X,-iXy) 

so2 

.P  0  ..  . 

-l^(X^  iXy) 

'73  ' 

The  conduction  band  dispersion  relation  was  worked  out  by  Kane  as 

f  Po^Yk^f  2  1 


2m 


fi^k^ 

2m\ 


E  =  E  +  ^  + 1  —2—  I  —I  — 

^  I  ^  I  o  I  ir 


3  yEg  Eg*  /S 


(9.15) 
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where  A  is  the  split-off  energy  between  the  light-hole-  and  split-off-band  maximas.  P  < 
can  then  be  solved  for  to  find  its  dependence  on  the  known  conduction  band  effective 
mass,  and  is  given  by 


2  m^Egf  m„ 

“  2  [m]  J3£p+2A 


(9.16) 


The  rate  of  change  in  the  number  of  photons  is  given  by  the  summation  over  pairs 
of  states  of  the  emission  minus  the  absorption  rates,  each  multiplied  by  the  their 
respective  probabilities  governed  by  Fermi-Dirac  statistics.  The  rate  is  given  by 


(9.17) 


where  / ,»  /  „  are  the  Fermi-Dirac  distribution  functions  of  the  upper  and  lower  states, 


respectively.  Expanding, 

d  Y-  ne^ 

2-  2  7/ 


^6[^„-£„]6[£„-£„-fico]/„(l-/„) 


(9.18) 


where  the  first  term  in  the  change  in  rix^  is  due  to  stimulated  emission  and  absorption, 

and  the  second  term  is  due  to  spontaneous  emission.  The  absorption  coefficient  is 
given  by 


Xxpcnfco) 


(9.19) 


which  is  related  to  the  imaginary  part  of  the  dielectric  constant  via  Equation  (4). 
Solving  for  e  "  ( CO  )  gives 
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(9.20) 


e"(oo) 


gp 


mluo^V 


n  m 


9.1.3  Oscillator  Strength  and  Lineshape  Broadening 


To  link  Equation  (9.20)  to  the  oscillator  model  of  Equation  (9. 12),  examine  the 
response  between  two  bands  n  and  mat  wavevector  k.  Let  £„-£„  =  nu)  1, 

/  „  =  0,  and  consider  only  one  polarization.  Equation  (9.20)  then  becomes 


rrioUi^V 


nhQ^  2\P„nx\‘ 


2mn(x>V  m«ha) 


6[hcOo-/icxj] 


nhe" 


ZrUoiXiV 


(9.21) 


The  factor  /  „mxis  called  the  oscillator  strength  which  describes  the  relative  strength  of 


the  oscillator  between  band  mand  band  n  for  the  polarization  K  compared  to  the  other 
oscillators  in  the  entire  system.  The  oscillator  strength  obeys  the  sum  rulei^) 


I ■  (’•“) 

n,  m 

If  the  oscillator  in  Equation  (9.21)  is  allowed  to  be  damped  as  in  Equation  (9. 1 1), 
the  delta  function  can  be  replaced  by  a  Lorentzian  lineshape, 

g  -  Ti(x>)  =  nb{uo g  -  uo) 


(a>„-cu)^  + Y' 


so  Equation  (9.21)  becomes 


6"(00)  = 


eVn 


nm\ 


2mgOoV  (a)„-ou)^  +  Y 


nm\ 


(9.23) 


(9.24) 
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Re-introducing  the  sum  over  discrete  states, 


e"(oo) 


2m„iK>V  „  „ 


ZII 


^ninxYn/TiX  7*  i 
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r  (C«o-lA))^  + Y„m).  (9.25) 

In  a  real  system,  the  lifetimes  of  the  states  are  finite,  so  the  energy  (and  hence  the 
frequency)  is  broadened  due  to  the  uncertainty  principlei**i/iAa)T  ~Ti,  or 


2Y„m„  =  Ao)  «-  =  —+  — 


1_  1 

t  T, 


(9.26) 


This  lifetime-broadening  is  referred  to  as  natural  broadening  if  x  is  radiative  only.  The 
effects  of  collision  broadening  have  been  neglected  in  Equation  (9.26),  although  they 
might  be  important  (even  dominant)  since  the  interaction  time  of  the  collision  can  be 
several  orders  of  magnitude  lower  that  the  mean  radiative  lifetimes.  Fortunately,  y  can 
be  determined  experimentally  from  absorption  or  time-resolved  spectroscopy  data. 

Equation  (9.2S)  is  particularly  useful  for  discrete  features  such  as  excitons,  while 
Equation  (9.20)  can  be  used  to  describe  systems  with  a  continuous  density  of  states. 
Broadening  mechanisms  can  be  included  in  the  later  case,  however,  by  convolving  it 
with  the  broadened  lineshape.i*^ 


9.1.4  Dielectric  Constant 

The  dielectric  constant  can  be  calculated  for  a  bulk  semiconductor  using 
Equation  (9.20).  For  each  band,  with -  1 ,  / „  =  0, 
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(9.27) 


The  conservation  of  momentum  is  contained  in  the  explicit  value  of  A:  for  and 


£■  min  the  conservation  of  energy  delta  function.  The  reduced  3-dimensional  (3D) 
density  of  states,  p  3,  contains  the  reduced  mass  for  bands  m  and  n, 

1  _  1  ^  1 

li  rn„  ^^  28) 

and  only  one  spin  is  considered  for  each  band  pair. 

To  get  the  real  part  of  the  dielectric  constant,  the  Kramers-Kronig  transformation 
can  be  used,  as  long  as  the  complex  part  is  known  over  a  sufficiently  large  frequency 
spectrum.  If  the  spectrum  of  interest  covers  a  small  region,  but  is  sufficiently  sqiaiated 
from  other  features  in  the  spectrum  which  contribute  to  the  local  real  dielectric 
constant,  the  transformation  can  still  be  used  with  proper  precautions.!*^)  In  this  case, 
the  significant  contributions  from  other  parts  of  the  spectrum  must  be  included  in  the 
background  dielectric  constant,  Ct,(uo),  which  is  real  and  replaces  £„  in  Equation  (9.6). 

If  only  a  small  spectral  region  is  known,  as  in  the  case  of  a  band  edge  where  only 
the  parabolic  part  near  k  =  0  is  assumed,  the  Kramers-Kronig  analysis  will  have  large 
errors.  To  get  around  this  problem,  the  real  part  of  Equation  (9. 10)  can  be  used  with 
€0  =  66(00)  and  the  real  oscillator  term  multiplied  by  the  oscillator  strength,  or 
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(cOo-CA)^)  +4y^oo  (9.29) 

For  the  quantum  well  dielectric  constant,  we  start  with  the  second  step  of 
Equation  (9.27),  and  substitute fcd/c  =  C\i/1i^)dE,  and  proceed  as  follows  with  the 
assumption  A: , ,  =  0,  and  k  =  which  is  quantized  in  the  5 direction. 
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(9.30) 


where  0  j  is  a  step  function  for  each  quantized  level,  E ,,  and  p  ?,d  is  the  reduced  density 


of  states  for  a  finite  quantum  well  with  only  one  spin  considered  (see  Appendix  C).  In 
fact,  the  reduced  density  of  states  in  Equation  (9.30)  is  valid  for  all  systems  as  long  as 
the  proper  form  of  A  z  is  used.P’> 
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9.2  Appendix  B 

Optical  Matrix  Elements  in  Quantum  Wells 


Under  the  envelope-function  approximation,  the  wavefunction  of  an  electron  in 
the  conduction  band  of  a  quantum  well  can  be  expressed  asi^i 

‘  ‘  (9.31) 

where  ;Sis  the  direction  of  quantization,  n  denotes  the  band  with  subscript  i  indexing  the 
bound  state,  <l>  (  z  )  is  the  envelope  wavefunction  and  li  (  r  )  is  the  periodic  Bloch 
function  of  the  bulk  material.  The  momentum  matrix  element  between  the  conduction 
band  (n)  and  valence  band  (m)  with  a  photon  field  polarized  in  the  direction  is 


^  f  -in  Ju 


,-.1.  - 


■  mil 


.(r)d^-V,u„(r)d  rj6(A:„||-k, 

.  (9.32) 

The  first  term,  <  I  mj  >is  the  overlap  of  the  envelope  wavefunctions  which  was 


assumed  to  be  6  by  Asada  er  ah  However,  <  c ,  |  /i/i  i  >  *  0 . 8  for  a  one-monolayer 
GaAs/Alg  jGag  7AS  quantum  well,  so  the  approximation  by  Asada  et  ah  is  no  longer 
valid  for  ultra-thin  wells.  /’  xmn  is  related  to  the  bulk  matrix  element  except  for  its 
average  with  respect  to  the  direction  of  k,  and  the  delta  function  is  the  momentum 
conservation  term. 

Using  the /c  •  p  method,  Kane  calculated  the  wavefunctions u(r).  If/c  =  /c;^the 


normalized  wavefunctions  can  be  written  in  spinor  rotation  asP'i 
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(9.33) 


The  basis  function  |  S  >  is  a  symmetrical  s-Iike  wavefunction,  while  \  X>,  IV  >,  and 


I Z  >  are  atomic  p-like  functions  antisymmetric  in  x,  y,  and  z,  respectively.  The 
subscripts  c,h,l  and  s  are  for  conduction,  heavy-hole,  light-hole,  and  split-off  bands, 
respectively,  and  T  and  I  denote  spin-up  and  spin-down  functions.  Because  of  the 
symmetry  of  the  basis  functions,  the  only  allowed  integrals  are 


-in<s i~i A'>--in<5 i-^ix >=-in<s i^iz . 

dx  bx  bx 

Finally,  the  matrix  elements  between  the  Bloch  wavefiinctions  are 

If  k  9^  k  5,  then  the  basis  functions  are  denoted  by  a  prime  ( I  A"  >)  and  are 


(9.34) 

(9.35) 


transformed  to  the  original  functions  (see  Figure  9.1)  usingi^’J 
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(9.36) 

5 '  =  5‘  . 

The  matrix  element  is  a  measure  of  the  dipole  moment  between  the  hole  in  band  m  and 


the  electron  in  band  n,  which  forms  a  vector,  P,  perpendicular  to  k.  The  interaction  of 
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Figure  9.1  Coordinate  system  for  the  transformation  of  basis  functions.  In  a 
quantum  well,  k,  is  taken  as  the  quantized  direction,  and  k^  is  fixed  at  k„  for  each 

energy  level  E^,.  For  transitions  in  the  n*  level  with  \k  \  >  ,  the  angle  0  is 

fixed  also. 
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the  dipole  with  the  electric  field  is  described  by  |  p .  ^  ^  1 2,  which  is  averaged  over  all 
directions  of  A;  in  bulk.  In  the  case  of  f  £  x,  the  matrix  element  for  the  conduction- 
heavy  hole  band  transition  is  derived  from  the  x  component  of  Pci  hi. 
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^x'^cihi  = -p(cos0cos(l)-isin(ti) 

v  2 

The  average  matrix  element  for  the  bulk  is  then 
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(9.38) 


Similarly,  the  matrix  elements  for  the  other  bands  can  be  calculated,  and  are  given  in 
Table  9.2. 


In  a  quantum  well,  however,  is  fixed  for  each  state,  so  for  a  particular 

magnitude  I  k  1, 0  will  also  be  fixed.  The  average  matrix  element  is  then,  for 
TE-polarized  photons  (d^  •  2  =  0), 
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(9.40) 


In  general,  the  effective  mass  in  a  superlattice  is  anisotropic,  having  a  larger  mass 
in  the  direction  of  the  superlattice  due  to  the  modified  dispersion  relation.  If,  however, 
the  mass  is  nearly  isotropic,  the  following  approximation  can  be  made  using 
A:,  =  kcosQ. 


/c^cos^G 


=  cos^G 


E  k-  (9.41) 

where  £  „  is  the  bound  state  energy  (A:  *  ”  A:  y  -  0),  and  E  is  the  total  energy  for  the 
three-dimensional  (3D)  system.  In  this  case.  Equation  (9.39)  becomes 
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(9.42) 


The  calculations  for  the  other  bands  are  given  in  Table  9.2.  The  approximation  of 
Equation  (9.40)  can  be  used  for  cos  ^  G,  while  si  n  ^  G  can  be  replaced  by  ( 1  -  £  „  /  £ ). 
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Table  9.2.  Bulk  and  Quantum- Well  average  matrix  elements <  >for 
conduction-heavy  hole,  light  hole,  and  split-off  band  transitions.  Note  that 
matrix  elements  for  T  T  transitions  are  the  same  for  i  i  transitions,  and  the 
i  T  transitions  are  the  same  for  T  i . 


X 

Bulk 

Quantum  WeU 

Spins 

T  t.  i  1 

ti .  it 

t  t,  i  i 

ti .  it 

c-hh 

TE 

-P^ 

3  ^ 

0 

ip^Ci+cos^e) 

0 

TM 

3  “ 

0 

^F^sin^e 

0 

c-lh 

TE 

-P^ 

9  ° 

9  “ 

^F^sin^e 

^Ff(l-cos"e) 

TM 

-P^ 

9  “ 

i/>2 

9  ° 

|F^cos^e 

X^oSin'e 

O 

c-so 

TE 

-P^ 

9^0 

2/>2 

9  “ 

O 

iF^(l+cos^0) 

D 

TM 

m 

-P^ 

9  ^ 

^Ffcos^e 

^F^sin'e 

Note  that  if  all  matrix  elements  for  each  transition  are  added  together,  including 
all  combinations  of  spin  and  polarization,  the  total  would  be  4/^0.  The  total  includes 
the  familiar  spin  degeneracy  as  well  as  the  sum  over  both  polarizations.  The 
polarization  selection  rules  can  be  obtained  from  Table  9.2  by  comparing  the  matrix 
elements  for  each  polarization.  For  TE,  the  selection  rules  give 


^  ‘  ch'^TE  •  ^  cl  ^  Tf  ‘  CJ  T£  O  •  1  •  ^ 


and  for  transverse  magnetic  (TM), 
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(9.43) 

(9.44) 
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In  the  conduction  to  heavy-hole  transition,  only  transitions  between  states  with  the 
same  spin  are  allowed,  thus  conserving  angular  momentum.  The  result  of  this  spin 
selection  rule  is  that  an  additional  factor  of  ^  must  be  multiplied  to  <  Pc/t  since  the 

probability  of  the  two  states  involved  having  the  same  spin  is^. 

Due  to  the  relativistic  motion  of  electrons  near  the  atomic  nuclei,  an  effective 
magnetic  field  is  present,  which  gives  rise  to  spin-orbit  interaction.  The  result  is  the 
lifting  of  the  degeneracy  of  the  valence  bands  so  that  the  light-hole  and  split-off-hole 
bands  are  separated  from  the  heavy-hole  band.P‘1  Hence,  the  spin  selection  rule  is 
relaxed  for  these  two  bands,  and  transitions  are  possible  between  states  with  opposite 
spins. 

To  complete  this  discussion  of  the  matrix  element,  we  must  determine  Po-  In  the 


Kane  model,  the  conduction  band  dispersion  is  given  by 

,  ^  n^k^Plfz  I 
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(9.45) 


which  is  the  origin  of  the  effective  mass,  m*.  In  Equation  (9.44),  £  „  is  the  bandgap  at 

k  =  0,  A  is  the  spin-orbit  splitting  of  the  valence  band,  and  is  the  free-electron 
mass.  is  then  related  to  the  effective  mass  in  the  conduction  band  by 
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In  Al^Ga,.,As,  the  values  of  and  AP^i  are 
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Eg=  1.519+  1.247X- 


5.405X 
7  +  204 


;x<0.45 


=  1.519+  1.247X+  1.147(x-0.45)^-^^^^^^^—^  ;x>0.45 
A  =  0.341  -0.066X 


(9.47) 

(9.48) 

(9.49) 


—  =  0.0665 +  0.0835  jc 

TTX  Q 

where  Eg  is  for  the  T  band,  and  both  Eg  and  A  are  given  in  eV.  Using  these  numbers 


and  finding  7 o  for  x  =  0  and  0.3,  for  example,  we  get 

Pl(ix  =  0)  =  7A8Egm, 

=  10.28X  10‘^° 
7^(jc  =  0.3)  =  S.22£*„m„ 

=  8.94x  10'^° 


with  units  of  eV-kg.  The  matrix  element  for  the  GaAs  in  this  case  is  1. 15  times  that  of 
the  alloy. 

For  a  quantum  well,  the  envelope  wavefunction  generally  extends  into  the  barrier, 
especially  in  very  thin  wells  where  most  of  the  wavefunction  is  outside  the  well.  For 
such  a  structure,  neither  the  envelope-wavefunction  overlap,  nor  the  matrix  elements 
for  the  Bloch  wavefunctions  in  the  barriers  is  negligible.  To  include  these  effects,  the 
matrix  element  given  in  Equation  (9.31)  must  be  found  for  each  region  (the  well  and 
both  barriers),  then  each  part  must  be  added  to  yield  the  total  matrix  element  for  the 
entire  structure.  The  matrix  element  is  then 
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where  (including  the  spin  selection  rule),  assuming (TE polarization)  and 

k^  =  ky~0. 


<p2  >'=ip2 


and  Pol  is  given  by  Equation  (9.46)  using  the  appropriate  parameters  from  that  region. 
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9.3  Appendix  C 

The  Density  of  States  for  Quantum  Wells  and  Superlattices 

M.W.  Prairie  and  R.M.  Kolbas,  "A  General  Derivation  of  the  Density  of  States 
Function  for  Quantum  Wells  and  Superlattices,"  Superlattices  and  Microstructures, 
Vol.  7,  No.  4,  pp.  269-277,  (1990). 
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The  intent  of  this  paper  is  to  provide  the  reader  with  a  detailed  summary  of 
the  development  of  the  density  of  states  (DOS)  functions  for  two-dimensional  sys¬ 
tems.  Spttifically,  the  DOS  is  derive  for  an  infinite  quantum  well,  a  finite  well, 
and  a  periodic  array  of  coupled  wells  (a  superlattice).  Many  authors  state  that  the 
DOS  is  ’simply  ...’without  references,  yet  many  who  are  new  to  the  subject  of  two- 
dimensional  systems  may  not  see  the  ’simplicity,’  for  instance,  of  the  derivation  of 
the  DOS  for  a  superlattice.  We  also  show  the  relationships  between  the  expressions 
for  each  case  when  the  appropriate  limits  are  taken.  This  comparison  shows  the 
consistency  that  such  a  general  derivation  furnishes  to  each  expression. 


Introduction 

The  distribution  of  carriers  (electrons  or  holes)  per 
unit  energy  interval  in  a  semiconductor  is  determine 
by  the  density  of  states  (DOS)  and  the  Fermi-Ditac  dis¬ 
tribution  functions.  The  DOS  function  tells  us  how 
many  electron  or  hole  states  are  in  an  interval  of 
energy,  while  the  Fermi-Dirac  function  tells  us  the 
probability  of  finding  a  carrier  for  each  of  those  states. 
For  instance,  in  the  case  of  a  superlattice,  the  distribu¬ 
tion  of  electronic  sutes,  or  the  DOS,  is  somewhere 
between  that  of  a  quantum  well  and  the  bulk 
semiconductor,  and  the  absorption  and  emission  line- 
shapes  are  greatly  affected  due  to  the  form  of  the  DOS 
function.'  This  suggests  that  an  understanding  of  the 
distribution  of  states  is  crucial  in  the  development  of 
optoelectronic  devices  employing  two-dimensional  sys¬ 
tems. 

In  this  review,  the  DOS  function  is  first  derived  for 
the  simple  case  of  the  ideal  (infinite),  two-dimensional 
(2D)  square  well.  This  pure  2D  an^ysis  is  found  to  be 
inadequate  and  is  improved  with  the  incorporation  of 
the  third  dimension  (the  well  width).  The  next  section 
considers  the  DOS  function  for  thin  or  shallow  wells 
with  the  bound  states  of  interest  near  the  top  of  the 
well.  In  this  case,  the  ideal  analysis  breaks  down,  and 
the  bound  states  must  be  found  numerically  and 
inserted  into  the  DOS  function.  The  effects  of  cou¬ 
pling  between  an  array  of  finite  quantum  wells  is  con¬ 
sidered  next  in  the  derivation  for  superlattices.  Finally, 
the  equivalence  between  each  case  is  shown  as  the 
appropriate  limits  are  taken.  First  is  an  example  of  the 
equivalence  of  the  single  well  DOS  to  that  of  the  bulk 
semiconductor  as  the  well  width  is  allowed  to  go  to 
infinity.  Next,  a  conceptual  analysis  is  used  to  show 
that  the  superlattice  DOS  reduces  to  the  single  quan¬ 
tum  well  case  as  the  barrier  width  gets  large,  and  the 


three-dimensional  (3D)  DOS  is  achieved  when  the 
barrier  width  goes  to  zero,  or  the  well  width  goes  to 
infinity. 


The  Infinite  Potential  Well 

The  density  of  states  (DOS),  g  ( f ),  is  defined  such 
that  the  number  of  orbital  states  per  unit  volume  with 
energy  between  E  and  d  E  is*’^ 

where  i  is  the  number  of  dimensions,  dk  is  the  differ¬ 
ential  volume  (3D),  area  (2D)  or  length  (ID)  element 
for  a  surface  of  constant  energy,  and  the  summations 
are  taken  over  spin,  degenerate  band  minima,  and  any 
other  mechanism  resulting  in  a  degeneracy  of  elec¬ 
tronic  states.  For  energy  surfaces  with  a  high  degree  of 
symmetry  (spherical,  parabolic,  etc.),  the  differential 
volume  element  can  be  coverted  to  a  differential  func¬ 
tion  of  energy,  and  the  DOS  is  obtained  by  taking  the 
derivative  of  both  sides  of  the  above  equation  wiSi 
respect  to  E.  However,  since  such  symmetry  is  not 
always  the  case,  we  will  use  the  general  expression  for 
the  DOS,  which  is  given  by^'^ 

^  (I) 

where  d  5  is  the  differential  surface  area  of  the  con¬ 
stant  energy  surface  in  3D,  while  we  use  a  differential 
length  for  the  2D  case. 
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The  energy,  E,  of  an  electron  near  the  minima  of 
the  conduction  band  in  the  free  particle  approximation 
is  given  by 


where  Ti  is  Plank’s  constant  divided  by  2n  ,  A:  is  the 
wavevector  of  the  particle,  and  m '  is  the  effective 
mass. 

For  the  pure  2D  case  with  only  spin  degeneracy, 
spherical  energy  surface,  and  -  kl*  k^. 
Equation  (1)  becomes 


gaff)  2/(2n)^tVt£- 

=  r _ 

Jo 


'•2* _ 2kd^ _ 

0  + 


Zn^h^k/m" 


Equation  (3)  is  the  DOS  for  the  pure  2D  case 
where  there  is  no  &  dimension,  and  bound  states 
higher  than  n= 1  are  neglected.  In  reality,  there  is  a  £ 
dimension,  which  in  fact  is  responsible  for  the  quantum 
size  effects  in  a  real  device.  To  find  the  DOS  using  the 
more  general  approach,  we  must  solve  the  problem  in 
three  dimensions.  The  energy  is  then 


Figure  1.  Differential  surface  area  of  a  constant  energy 
surface  in  cvlindrical  coordinates.  The  x-y  cross  section 
is  circular,  but  the  x-z  (or  y-z)  cross  section  is  not  nec¬ 
essarily  circular. 


Substituting  this  expression  into  d  S,  we  get 


dS  - 


i(  k. 


(dk.)= 


jk^  * k\dk^d^ . 


The  3D  single  quantum  well  (SQW)  DOS  is  then 


•^.ik^*kl),  (4 

Zm 

where  the  electron  is  free  in  the  x  and  y  directions 
(k ,  and  Ac ,  are  continuous),  while  it  is  bound  in  the  i 
direction  ( Ac ,  is  discrete).  With  the  help  of  Figure  1, 
the  differential  surface  area  is 


dS~kd^yj(,dky*(,dk,f.  (5 

Since  the  differential  change  in  energy  along  a  constant 
energy  surface  is  zero, 

dE - :(,kdk*  k,dk,)~  0  . 

m 


/•*»  r 

J Jo 


2"  Z^k^*  kld^dk. 


‘.•^0  8n^|v,[^(Ac"*A:^)]| 
*•  J k^  *  kldk. 


Recall  from  the  solution  of  the  infinite  square  well 
problem  that  k ,  is  quantized,  and  is  given  by 
A:  j  ■=  nn/  L,  where  L ,  is  the  width  of  the  well.*  Insert¬ 
ing  this  into  Equation  (8)  gives 


dk  -  -  • 


Superlattices  and  Microstructures,  Vol.  7.  No.  4. 1990 


271 


Comparing  this  result  to  Equation  (3),  we  see  that 


The  Finite  Potential  Well 


n- 1.2.3....  (10) 


This  same  result  can  also  be  obtained  starting  with  the 
3D  DOS.  Using  ^nation  (1),  the  3D  DOS  for  a  bulk 
semiconductor  is  given  by  the  well-known  expression^ 


(11) 


The  solution  to  the  infinite  square  well  (ideal  case) 
results  in  quantized  energies,  and  for  both  A; ..  and  k  ^ 
equal  to  zero,  is  given  by* 


2m'L\ 


n-  1.2.3....  (12) 


Using  this  form  of  the  energy  in  Equation  (11),  the  3D 
DOS  becomes 


SpU' 


(£■) 


v^(m')^  nfin 


m'  n 


(13) 


which  is  identical  to  Equations  (9)  and  (10). 

Equations  (10)  and  (12)  are  shown  g^hit^ly  in 
Figure  2  below.  Note  for  each  step  in  n,  a  constant 
92(f)/Z,  is  added  to  the  3D  DOS  function,  which 
results  in  a  set  of  steps  that  follow  the  bulk  DOS  curve. 


Figure  2.  The  DOS  of  an  infinite  quantum  well  as  it 
corresponds  to  the  bound-sUte  energies.  In  this  case, 
the  bound-state  energy  is  inversely  proportional  to  the 
well  width,  while  the  DOS  is  inversely  proportion^  to 
the  square  of  the  width,  so  a  well  that  is  half  as  thick  as 
the  original  one  will  have  a  bound  state  twice  as  high, 
and  a  DOS  four  times  the  original  value. 


We  have  seen  so  far  that  the  DOS  of  an  infinite 
potential  well  is  a  step-like  function,  where  a  constant 
g2iE)/L,  is  added  to  gf*''^(f)  for  each  quantized 
step  in  energy.  For  the  finite  well,  however,  the  wave- 
functions  near  the  top  of  the  well  are  less  confined, 
resulting  in  lower  bound  state  energies  compared  to 
those  of  the  infinite  well,  resulting  in  different  size 
steps  in  the  DOS.  This  can  be  seen  by  solving 
Schrodinger's  Equation  and  numericallv  or  graphically 
finding  the  values  for  the  bound  states.' 

The  finite  potential  well  is  shown  in  Figure  3.  Tak¬ 
ing  the  bottom  of  the  well  as  zero  potential  and  the  tc^ 
of  the  well  as  K  o  .we  can  define  the  following 
variables; 


2  2m-(K.-f) 

a  - - — - 


^2  2m‘E 

2  2m‘V„ 

Y  n— 


For  even  vravefunctions,  the  solution  is 

• 

(14) 

while  the  odd  solution  is 

a  -  -kcot( k~ 

\  2, 

)• 

(15) 

Also,  notice  that  from  the  defined  variables, 

*  k^ , 

(16) 

1(2) 
I  p 


// 


m 


Figure  3.  Finite  potential  well.  The  width  is  I,  ,  and 
the  potential  depth  is  I '  o  • 
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which  is  the  equation  for  a  circle  with  radius  Y  ■ 
Graphing  equations  (14-16)  as  in  Figure  4,  we  can  find 
the  set  of  allowed  values  k  „ ,  which  give  us  the  values 
of  the  energy  eigen  states 


2m 


-Xn 


(17) 


Using  Equation  (17)  in  (11),  we  find  the  DOS  for 
the  finite  quantum  well  (FQW)  to  be 


m 


(18) 


0  i,  9(i) 


Figure  5.  C^omparison  of  energy  states  and  DOS  for  an 
infinite  single  quantum  well  (SQW)  and  a  finite  quan¬ 
tum  well  (FQVO.  Notice  that  the  energies  of  the  states 
in  the  finite  well  are  no  longer  inversely  proportional 
the  well  width,  and  the  DOS  is  not  inversely  prq>or- 
tional  to  the  square  of  the  width. 


Note  that  in  this  case,  the  DOS  for  a  thin  well  has 
nearly  the  same  DOS  as  a  well  twice  as  big  as  the  first. 
This  IS  contrary  to  the  infinite  well  case,  where  the 
DOS  of  the  first  state  of  the  thinner  well  is  four  times 
as  great  as  the  DOS  of  a  well  that  is  twice  as  wide! 

The  consequence  is  that  in  a  real  system  with  very  thin 
wells,  calculations  which  depend  on  the  density  of 
states,  such  as  the  capture  or  emission  of  charge  carri¬ 
ers,  b^me  less  depmdent  on  the  width  of  the  well. 

In  this  case,  other  factors  such  as  the  spatial  extent  of 
the  wavefunction  become  more  important.' 


Figure  4.  Graphical  solution  for  bound  states  in  a 
finite  potential  well.  As  the  well  gets  thinner,  the  cir¬ 
cular  curve  of  radius  Y  will  contract.  In  this  case  (sym¬ 
metric  well),  there  will  always  exist  at  least  one  bound 
state. 


Notice  that  Ak  ^  k„.,  -  k„  is  no  longer  a  con¬ 
stant.  In  fact,  A  k  gets  smaller  as  n  gets  bigger  for  the 
same  reason  that  each  bound  state  of  the  finite  well  is 
lower  than  the  corresponding  eigen  ene^y  of  the  infi¬ 
nite  well  by  an  increasing  amount  as  n  increases  as 
shown  in  Figure  3. 

Notice  that  the  highest  bound  state  of  the  finite 
quantum  well  is  very  near  the  top  of  the  well.  Refer¬ 
ring  back  to  Figure  4,  the  intersection  on  the  circle  for 
this  state  occurs  at  a  point  where  a  tangent  to  the  circle 
will  be  almost  vertical.  This  means  that  if  the  width  of 
the  well  were  changed,  the  shift  in  the  bound  state 
would  be  very  small.  If  the  well  was  very  thin  such  that 
n/ Lz  »  y,  we  would  only  have  one  bound  state  in  the 
well  and  the  energy  of  that  state  would  be  near  the  top 
to  the  well.  In  this  case,  a  large  change  in  L ,  would 
only  produce  a  small  shift  in  the  energy,  as  shown  in 
Figure  6. 


SuperlatUces 

When  two  finite  quantum  wells  are  placed  near 
each  other  such  that  the  barrier  between  them  is  small 
enough  to  allow  interaction  between  the  wells,  the 
degenerate  states  split  into  two  bound  states  of  nearly 
the  same  energy.  If  this  multiple  quantum  well 
arrangement  were  extended  to  many  wells  placed 
together,  then  the  energy  levels  broken  into  bands, 
much  like  the  ene^y  bands  of  the  host  crystal  which 
are  due  to  the  peri^icity  of  the  atomic  lattice.  The 
differences  are  that  the  quantum  well  lattice  is  one 
dimensional  with  a  period  several  times  greater  than 
that  of  the  crystal,  and  the  energy  bands  due  to  the 
coupled  wells  are  superimposed  onto  the  allowed 
energy  bands  of  the  host. 

Such  a  structure,  first  proposed  by  Esaki  and  Tsu  in 
1969,  is  called  a  superlatbce  due  to  the  periodic  poten¬ 
tial  superimposed  on  the  crystal’s  periodic  potential, 
which  gives  rise  to  a  new  hand  structure.’  The  new 
potential  is  due  to  an  array  of  square  wells  of  width  L , , 
separated  by  barriers  of  width  and  height  U  o . 
Applying  the  tight-binding  approximation  to  the  new 
periodic  potential,  the  superlattice  energy  bands  can  be 
calculated.  For  the  periodic  superlattice  potential 
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(a)  (b)  (c) 

Figure  6.  (a)  Graphical  solution  for  well  with 
2i,  “  i',,  <c  n/ L\,  (b)  Single  bound  state  for 
both  wells,  and  (c)  The  DOS  for  the  corresponding 
wells.  Notice  that  the  DOS  of  the  thin  well  is  not  four 
times  that  of  the  thick  well.  In  fact,  for  ultra-thin  wells, 
the  DOS  of  both  wells  may  be  nearly  equal. 


U  (2)  -  U(_z-  nd),  whered  is  the  period  of  the 
superlattice,  Lt,  and  n  is  an  integer,  the  energies 
can  be  found  usingb* 

<tu.  I  I  \l;„> 

<V, 


- 


I 

for  N  wells.  Solving  (19)  yields 

-  f^‘^'^-o-2Pcos(k,d)  (22) 


a-  £(H;"«'"(2))’(/(2)H;"’‘"(2)d2  .  (23) 

which  is  the  shift  from  E  due  to  the  coulombic 
interaction  of  the  coupled  wells,  and  P  is  the  exchange 
energy  which  determines  the  extent  of  the  band  broad¬ 
ening,  and  is  given  by 

P- J(,|;"<'‘^(2-d))'i/(2)t(;''‘'‘^(2)d2  .  (24) 


assuming  appreciable  interaction  from  adjacent  wells 
only.  The  important  thing  to  note  about  a  and  P  is 
that  as  the  barrier  width  decreases,  the  interaction 


between  wells  should  increase,  and  thus  a  and  P 
should  increase. 

The  energy  bands  of  a  superlattice,  as  shown  in 
Figure  7(a),  can  be  found  using  Equation(19),  with  the 
to^  energy  written  as 

(d  TTX 

f.2 

•  Eo-2^cosik,d)*—~.{k\*k\)  (25) 
cm 

where  Eo~  E  ~  a .  The  divergence  of  the  total 
energy  with  respect  to  the  wavevector  is 

Vj£  “  2dpsin(kjd)i  *  ^(A:,x  + /c,y )  ■ 


so  the  magnitude  is  then 


The  differential  surface  area  of  the  constant  energy 
surface  is  given  by  Equation  (S).  Recall  that  the  differ¬ 
ential  change  in  energy  of  a  cotuiaru  energy  surface  is 
zero,  so  we  can  relate  d  k  to  d  k  I .  lMk^-k\*k\, 
then  d  f  -  0  becomes 

dE  ~0~2^dsix\(,k,d)dk,*—kdk  . 

m 

so 

dk  ~ -^^^J—s\n(,k^d)dk,  .  (27) 

The  differential  surface  area  of  Equation  (5)  then 
becomes 
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Figure  7.  (a)  Dispersion  relation  for  a  superlattice, 
and  (b)  DOS  for  a  superlattice  versus  single  well  and 
bulk.  Because  the  interaction  between  wells  causes  the 
bound  states  to  broaden  into  bands,  the  DOS  is  no 
longer  constant  as  with  the  finite  well,  but  it  does 
retain  a  step-like  nature. 


dS‘kdi 


^Jsin\k.d}(dk,y*tdk,y  (28) 


which,  using  Equation  (26),  can  be  re-written  as 


dS-d^dk, 


f4d*p*sin*(k 


Putting  this  in  the  DOS  Equation  (1),  we  get 

_L  r_£LL- 

4n=’J  |V.f I 


k* 

(29) 


-  -^2n  f  '  dk, 

J-k. 


m 


2n^fi* 


’  f  'dk, 

Jo 


n^fi 


which  is  the  same  as  Equation  (8).  The  next  step  is  to 
solve  for  k„  but  normidize  the  problem  first  by  defin¬ 
ing  a  dimensionless  energy  that  is  zero  below  the  bot¬ 
tom  of  a  subband,  and  hu  a  value  of  one  for  energies 
above  the  subhand.  Using  Equation  (25),  and  adding 
and  subtracting  2p, 

h^k^  I  4B 

f-fo*— -i4Pcos(k.d)*^-2p.  (31) 

d  c. 

Rearranging, 

f-f„-^*2p  , 

- - -(1 -cosfc.d)-?  (32) 

where  ^  is  the  new  dimensionless  energy.  Now  solving 
for  k , ,  we  get 

k,-iarccos(l -2?)  (33) 

d 

The  DOS  for  a  superlattice  is  then 
m-  (  ° 

-2?):  0<?<1  (34) 

"  n  :  ?>1 

which  is  shown  in  Figure  7(b).  To  express  ( ^ ) as 
9  f  ‘  ( ? ),  we  must  note  that 


(30) 
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and 

So  Equation  (3S)  becomes 

gr(Ud^-gr(f)4pd5 


or, 

gf(?)“^^arccos(l-25):  0<5<1 


(35) 

(36) 

(37) 


Recall  that  this  analysis  was  for  each  subband,  so 
we  must  modify  'Jie  equation  to  represent  the  general 
case.  From  the  r.uperlattice  dispersion  relation  (Figure 
7(a))  we  see  that,  ^ong  with  Equation  (33),  that  the 
general  form  of  the  wavevector  is 

^.-^-^arccosd -25.)  (38) 


where  n  is  the  number  of  filled  subbands,  and  is  the 
dimensionless  energy  of  the  i"'  subband  where 
(  -  n  +  1  (recall  that  each  is  unique  since 

depends  on  £  .  a  and  P ,  all  of  which  change  for 

different  subbands).  The  total  superlattice  DOS  is 
then 

9f(f)~  z'?;  .['in-"arccos(l -2^.)]  (39) 

n  n  o. 


for  0  <?,<!.  For  energies  that  fall  between  the  n"* 
and  subband,  we  have 


m'  n 
nh^  d 


-g2(f)^  (40) 


which  looks  a  lot  like  Equation  (10)  (  gf‘’'^(£)),  except 
that  we  have  the  factor  of  I  /  d  instead  of  1  /  £ , .  This 
is  because  the  superlattice  DOS  depends  on  the  inter¬ 
action  of  coupled  wells  in  a  periodic  potential  of  period 
d,  while  the  single  quantum  well  DOS  depends  only  on 
L ,  since  there  is  no  interaction  from  other  wells. 

Equivalence  of  DOS  Expressions 

The  single  quantum  well  DOS  derived  earlier 
should  reduce  to  the  3D  DOS  if  the  well  width  is 
allowed  to  go  to  infinity.  As  an  example  of  this,  let  us 
find  the  number  of  electrons  per  unit  volume  in  a 
quantum  well  of  a  non-degenerate  semiconductor  (i.e., 
let  the  conduction  band  energy,  £c ,  be  greater  than  a 
few k,T  above  the  Fermi  energy,  E  f).  The  solution 
should  be  equivalent  to  the  electron  density  found 


when  the  3D  DOS  is  used  for  the  bulk  semiconductor 
case.  The  expression  for  the  election  density  is" 

(41) 

where  /  ^ o(  f  ) is  the  Fermi-Dirac  distribution  func¬ 
tion,  which  can  be  ^proximated  by  the  Maxwell- 
Boltzmann  distribution  function,  /  «« ( £* ),  for  a 
non-degenerate  semiconductor  where  £  c  -  £  f  is 
greater  than  a  few  kgT,  or 


-/«.(£)  (42) 


So  the  electron  density,  n„  is  then 


r  m  n 

n,  -  — r  —  e  d£ 

*  Jr,  nft"  £. 

•  Ll  £  -£_ 

-  e  *  I  ne  •  dE  (43) 

nfi  L,  Je, 


Notice  that  n  is  a  function  of  E  that  increases  by  one 
for  each  additional  bound  state.  For  the  bottom  of 
each  step  in  the  DOS  curve,  let  £-£,-£„,  and  let 
(£„,,  -  E c)  »  *,r,  and  take  the  contribution  of 
each  bound  state  separately,  as  in  Figure  8, 


n,  -  A*  B*  ... 


or 


•  IL  .±. 

'Tl  t  ‘t’’  I  f 

-TTT^^  c  '  c(£ 

nfi  L,  Ji, 

Ll  ..  .±. 

*  c**'  r  e  ‘‘"dE  *...  ] 

A^{\'^^dE  (44) 

^r\.h^L,  Jr, 


m'k.T 


Zni  Pig/  k.T 

- - —  e 

n  nfi^L^ 


fU  K  g  I  ^ 
- - -  )  ( 


(45) 
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where  N  f  is  the  effective  DOS  for  the  conduction 
band  electrons  in  a  quantized  system.  As  f.,  grows 
large,  the  bound  states  will  drop  toward  E,. ,  and  will 
eventually  form  an  energy  band,  as  in  the  bulk  mate¬ 
rial.  When  this  happens,  we  can  say  that  the  number  of 
bound  states  in  the  band  is  large  (note  in  Equation  (10) 
that  Z.i«n),  so  asZ.2  goes  to  infinity,  n  also  goes  to 
infinity,  and  we  should  get  the  same  N  c  as  for  the  3D 
bulk  case.  In  the  limit  of  the  sum  going  to  an  integral. 


n 


dn 


fn  1 2m‘  Llk,T 

V  2' V 


(46) 


Substituting  Equation  (46)  into  N  gives 


m‘ktT 

2nfi^ 


3 


(47) 


which  is  exactly  the  value  for  the  bulk  case,  where' ■ 


(48) 


Since  Equations  (47)  and  (48)  are  equivalent, 

-  N  in  the  limit  as  i,  goes  to  infinity.  This 
is  equivalent  to  saying  that  the  general  2D  DOS  is  the 
same  as  for  the  bulk  case  in  this  limit  since  the  only  dif¬ 
ference  between  these  expressions  is  the  DOS  func¬ 
tions  used  at  the  start  of  their  derivations. 

To  see  how  Equation  (40)  for  the  superlattice  DOS 
relates  to  the  single  quantum  well  case,  we 

must  move  the  wells  of  the  superlattice  apart 
(Z.6  -»  »),  which  removes  their  interaction.  This 
results  in  a  and  |3  both  going  to  zero,  so  Equation  (25) 
becomes 


E 


2m‘ 


Figure  8.  The  DOS  near  the  bottom  of  a  wide  quan¬ 
tum  well.  The  contribution  of  electrons  from  ea^ 
bound  state  is  found  when  each  element  of  the  DOS  is 
integrated  with  the  distribution  function,  giving  the 
total  number  of  electrons  as  the  sum  of  each  contribu¬ 
tion. 


fc")  . 


(49) 


which  means  that  k ,  *  n  it  /  /. ,  in  the  limit  of  large  bar¬ 
riers.  Using  this  wavevector  in  Equation  (31),  we  get 


lim  gf{E)  ■ 

I.— 


m 


nn 

17 


m’  n 
nh^ 


-g®'^^(f)  (50) 

Similarly,  to  show  how  the  superlattice  DOS  relates 
to  the  bulk  case,  gj^E)  ,  we  can  either  let  /.(,  vanish, 
or  let  Z.,  go  to  infinity.  If  Z.^  gets  small,  the  coupling 
between  the  wells  get  stronger,  and  the  subbands 
broaden  until  the  steps  in  the  E  vs  gf'(f )  curve 
merge  and  form  a  continuous  parabolic  curve.  Note 
that  the  symmetrical  shape  of  the  steps  in  Figure  7(b) 
are  arc  cosine  functions,  which  is  due  to  the  assump¬ 
tion  that  only  nearest-neighbor  interactions  contribute 
in  the  tight-binding  theory.  For  very  strong  coupling, 
however,  wells  beyond  the  nearest  neighbors  must  be 
consider^  since  the  spatial  extent  of  tte  electron 
wavefunctions  become  large.  This  af^rent  failure  of 
the  tight  binding  method'^  will  result  in  the  loss  of  the 
arc  cosine  symmetry  of  the  steps,  and  should  allow  the 
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shape  of  the  steps  to  conform  more  to  the  parabolic 
bulk  case. 

If,  however,  we  let  L ,  get  large,  the  carrier  wave- 
function  will  be  more  confined  in  the  wells,  which 
reduces  the  coupling  and  causes  the  subbands  to 
narrow.  But  increasing  Z.,  also  has  the  effect  of  mov¬ 
ing  the  subbands  down  in  energy  (remember  that  the 
subbands  are  nearly  centered  on  the  bound  states  of 
the  single  quantum  well),  which  would  cause  then  to 
stack  up  at  the  bottom  of  the  wells.  For  a  large  well 
width,  the  steps  in  the  E  vsgPiE)  curve  will 
become  so  small  and  numerous  that  they  will  form  a 
quasi-continuous  curve  that  will  approach  the  bulk 
DOS  case.  In  either  case,  the  energy  band  would 
approach  the  bulk  case,  where  F  ( k )  can  be  approxi¬ 
mated  by  a  parabolic  curve,  and  k ,  would  become 

Ac  j  *  yl2m’  E  /fi  .  Using  this  intuitive  analysis  in 
Equation  (38),  we  get 
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and  Equation  <’39)  becomes 
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which  is  the  3D  DOS  given  in  Equation  (11). 

Conclusions 


(52) 


The  derivation  of  the  DOS  function  for  quantum 
wells  and  superlattices  is  straight  forward  when 
approached  with  some  insight  into  the  physical  con¬ 
straints  of  these  structures.  The  pure  2D  infinite  well 
case  proved  to  be  inadequate  for  the  treatment  of 
actual  quantum  well  systems  because  it  ignores  the  well 
width  and  higher-order  bound  states  of  the  quantum 
well.  The  finite  well  case  was  obtained  when  the 
appropriate  boundary  conditions  were  applied.  Special 


care  must  be  taken  when  analyzing  very  thin  wells 
because  the  energy  levels  may  or  may  not  be  strongly 
dependent  on  the  well  width.  The  superlattice  DOS 
was  then  derived  by  considering  the  interaction  of  an 
array  of  finite  quantum  wells.  Finally,  we  showed  that 
each  of  the  above  cases  is  consistent  with  the  others 
when  the  appropriate  limits  of  well  or  barrier  width  are 
taken  into  account. 
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